Re: [R] MSBVAR Package
Hi Ahmad, I don't know of any, but this might help: http://maths-people.anu.edu.au/~johnm/courses/r/ASC2008/pdf/Rtimeseries-ohp.pdf Jim On Thu, Sep 28, 2017 at 2:12 PM, ah awrote: > Hi Jim > > Thank you very much for your reply and your help. > By the way, where can I find some tutorial videos? > Thanks again for your help. > > Best regards > Ahmad > __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] building random matrices from vectors of random parameters
I would try something like n = 5 a <- rnorm(n,0.8,0.1) so <- rnorm(n,0.5,0.1) m <- rnorm(n,1.2,0.1) mats = mapply(function(sa1, so1, m1) matrix(c(0,sa1*m1,so1,sa1),2,2,byrow=T), a, so, m, SIMPLIFY = FALSE) > mats [[1]] [,1] [,2] [1,] 0.000 0.9129962 [2,] 0.4963598 0.7067311 [[2]] [,1] [,2] [1,] 0.000 1.0150316 [2,] 0.5489887 0.8469046 [[3]] [,1] [,2] [1,] 0.000 0.9516137 [2,] 0.3724521 0.8306535 [[4]] [,1] [,2] [1,] 0.000 1.0525355 [2,] 0.8075108 0.8314638 [[5]] [,1] [,2] [1,] 0.000 0.9400074 [2,] 0.4803386 0.7901753 On Wed, Sep 27, 2017 at 5:47 PM, Evan Coochwrote: > Suppose I have interest in a matrix with the following symbolic structure > (specified by 3 parameters: sa, so, m): > > matrix(c(0,sa*m,so,sa),2,2,byrow=T) > > What I can't figure out is how to construct a series of matrices, where the > elements/parameters are rnorm values. I'd like to construct separate > matrices, with each matrix in the series using the 'next random parameter > value'. While the following works (for generating, say, 5 such random > matrices) > > replicate(5,matrix(c(0,rnorm(1,0.8,0.1)*rnorm(1,1.2,0.1),rnorm(1,0.5,0.1),rnorm(1,0.8,0.1)),2,2,byrow=T)) > > its inelegant, and a real pain if the matrix gets large (say, 20 x 20). > > I'm wondering if there is an easier way. I tried > >> sa <- rnorm(5,0.8,0.1) >> so <- rnorm(5,0.5,0.1) >> m <- rnorm(5,1.2,0.1) > > matrix(c(0,sa*m,so,sa),2,2,byrow=T) > > but that only returns a single matrix, not 5 matrices as I'd like. I also > tried several variants of the 'replicate' approach (above), but didn't > stumble across anything that seemed to work. > > So, is there a better way than something like: > > replicate(5,matrix(c(0,rnorm(1,0.8,0.1)*rnorm(1,1.2,0.1),rnorm(1,0.5,0.1),rnorm(1,0.8,0.1)),2,2,byrow=T)) > > Many thanks in advance... > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] building random matrices from vectors of random parameters
Suppose I have interest in a matrix with the following symbolic structure (specified by 3 parameters: sa, so, m): matrix(c(0,sa*m,so,sa),2,2,byrow=T) What I can't figure out is how to construct a series of matrices, where the elements/parameters are rnorm values. I'd like to construct separate matrices, with each matrix in the series using the 'next random parameter value'. While the following works (for generating, say, 5 such random matrices) replicate(5,matrix(c(0,rnorm(1,0.8,0.1)*rnorm(1,1.2,0.1),rnorm(1,0.5,0.1),rnorm(1,0.8,0.1)),2,2,byrow=T)) its inelegant, and a real pain if the matrix gets large (say, 20 x 20). I'm wondering if there is an easier way. I tried > sa <- rnorm(5,0.8,0.1) > so <- rnorm(5,0.5,0.1) > m <- rnorm(5,1.2,0.1) matrix(c(0,sa*m,so,sa),2,2,byrow=T) but that only returns a single matrix, not 5 matrices as I'd like. I also tried several variants of the 'replicate' approach (above), but didn't stumble across anything that seemed to work. So, is there a better way than something like: replicate(5,matrix(c(0,rnorm(1,0.8,0.1)*rnorm(1,1.2,0.1),rnorm(1,0.5,0.1),rnorm(1,0.8,0.1)),2,2,byrow=T)) Many thanks in advance... __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] build a SpatialLines object from a list
Have you tried following the example in ?'SpatialLines-class' You'll probably get better help from R-sig-geo And please don't send html email, it makes your email hard to read. -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 On 9/26/17, 2:36 AM, "R-help on behalf of Ashraf Afana via R-help"wrote: Hi Eric, Thanks for the help.But this will not solve the problem as it will generate a list and what I need is an object of class sp using SpatialLine function from sp package.So, I need to convert each matrix to coordinates and then to a line and then to a spatial line as figured in the code. My data structure is a list of 141 matrices.Each matrix represents coordinates of the river lines position. Ashraf, cheers On Monday, 25 September 2017, 16:56, Eric Berger wrote: Hi Ashraf,It is not obvious to me what your structures are but one problem in your function is the assignment tt1 <- SpatialLines(list(tt[[i]])). This will set tt1 to just have one item. Consider the following test.func <- function(x) {tt1 <- list()for ( i in ... ) { ... tt1[[i]] <- SpatialLines(tt[[i]])}return(tt1)} HTH,Eric On Mon, Sep 25, 2017 at 3:40 PM, Ashraf Afana via R-help wrote: Hi all,I'm trying to build a SpatialLines object from a list that contains 124 river segments. Each segment in the list contains the x,y coordinates. I'm using the following code to create the SpatialLines object, but it just retrieves one segment. Any suggestions? test.func = function(x){ for (i in 1:length(x)) {tt[[i]] <- x[i]; tt[[i]] = Line(tt[[i]]); tt[[i]] = Lines(list(tt[[i]] ), 'i')tt1 = SpatialLines(list(tt[[i]])) }return(tt1) } Ashraf, [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/ listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/ posting-guide.html and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] need held in r coding.
It's pretty clear from the error message that dat$Date is not an appropriate time-based object. Since dat$Date is created by the get_majorlandmarks() function, and your question provides no information about that, it's hard to be more specific. However, I would suggest looking at the output of class(dat$Date) and str(dat$Date) The results should indicate that dat$Date somehow indicates a time-based ordering. I wouldn't be surprised if dat$Date is a factor (the class() function will tell you if it is), in which case something is wrong with how dat is constructed. For example, there might be invalid values in an input file. The other thing is whether dates are acceptable to xts() for the order.by argument. I don't know. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 Lab cell 925-724-7509 On 9/27/17, 2:31 PM, "R-help on behalf of istiyak ahamad"wrote: Need Help in Debugging below script: dat <- get_majorlandmarks(dat,Dmin,Per) fit_xts <- xts(dat$fit,order.by = dat$Date,frequency = 365) close_xts <- xts(dat$Close, order.by = dat$Date, frequency = 365 ) majorlandmarks_xts <-xts(dat$Close[dat$majorlandmarks==TRUE], order.by = dat$Date[dat$majorlandmarks==TRUE], frequency = 365 ) minorlandmarks_xts <- xts(dat$Close[dat$minorlandmarks==TRUE], order.by = dat$Date[dat$minorlandmarks==TRUE], frequency = 365 ) landmarkschart_data <- cbind(close_xts,fit_xts,majorlandmarks_xts,minorlandmarks_xts) Error:- order.by requires an appropriate time-based object dygraphs::dygraph(landmarkschart_data, xlab = "Time",ylab = "Price") %>% dySeries("..1",label = "Close", drawPoints = FALSE, strokeWidth = 1, color = "black") %>% dySeries("..2",label = "Fit", drawPoints = FALSE, strokeWidth = 1, color = "red") %>% dySeries("..3",label = "Major LandMarks", drawPoints = TRUE, pointSize = 5, strokeWidth = 0, color = "darkorange") %>% dySeries("..4",label = "Minor LandMarks", drawPoints = TRUE, pointSize = 1, strokeWidth = 0, color = "darkblue") %>% dyRangeSelector() ###Estimating error with different parameters # RMSE <- vector() # for(i in seq(0.01,0.1,0.01)){ # dat <- get_majorlandmarks(dat,i,0.01) # RMSE[i] <- sqrt(sum((dat$fit-dat$Close)^2)/nrow(dat)) # } # plot(RMSE) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] MSBVAR Package
Hi Ahmed, You seem to know about the package, so I would suggest downloading: https://cran.r-project.org/web/packages/MSBVAR/MSBVAR.pdf getting a big cup of coffee and going through the examples to see how the package works. If you can't make any sense out of that, you will probably have to contact someone who does know how to use the package and is willing to provide a tutorial (that is not me). Jim On Wed, Sep 27, 2017 at 5:22 PM, ah a via R-helpwrote: > dear sirs or madam, > > As I'm interested to search about the monetary transmission channel in our > country by MSVAR model, I would be grateful if you help me and tell me how > can I run MSVAR in R or send me the related code to run this model . > Actually, I'm new user of R and I don't know how to run Markov Switching > Var Model in R. > > Thank you very much in advance for your help. > > Best Regards, > Ahmad Arabian > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] need held in r coding.
Need Help in Debugging below script: dat <- get_majorlandmarks(dat,Dmin,Per) fit_xts <- xts(dat$fit,order.by = dat$Date,frequency = 365) close_xts <- xts(dat$Close, order.by = dat$Date, frequency = 365 ) majorlandmarks_xts <-xts(dat$Close[dat$majorlandmarks==TRUE], order.by = dat$Date[dat$majorlandmarks==TRUE], frequency = 365 ) minorlandmarks_xts <- xts(dat$Close[dat$minorlandmarks==TRUE], order.by = dat$Date[dat$minorlandmarks==TRUE], frequency = 365 ) landmarkschart_data <- cbind(close_xts,fit_xts,majorlandmarks_xts,minorlandmarks_xts) Error:- order.by requires an appropriate time-based object dygraphs::dygraph(landmarkschart_data, xlab = "Time",ylab = "Price") %>% dySeries("..1",label = "Close", drawPoints = FALSE, strokeWidth = 1, color = "black") %>% dySeries("..2",label = "Fit", drawPoints = FALSE, strokeWidth = 1, color = "red") %>% dySeries("..3",label = "Major LandMarks", drawPoints = TRUE, pointSize = 5, strokeWidth = 0, color = "darkorange") %>% dySeries("..4",label = "Minor LandMarks", drawPoints = TRUE, pointSize = 1, strokeWidth = 0, color = "darkblue") %>% dyRangeSelector() ###Estimating error with different parameters # RMSE <- vector() # for(i in seq(0.01,0.1,0.01)){ # dat <- get_majorlandmarks(dat,i,0.01) # RMSE[i] <- sqrt(sum((dat$fit-dat$Close)^2)/nrow(dat)) # } # plot(RMSE) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Need Help in Debugging
I am getting following error when running this script :Error in length(runsum) : object 'runsum' not found ##error function dat$fit <- NULL dat$fit[dat$majorlandmarks] <- dat$Close[dat$majorlandmarks] run <- rle(dat$majorlandmarks) runvalue <- run$values runsum <- cumsum(run$lengths) run <- run$lengths for(i in 1:(length(runsum)-1)){ if (runvalue[i]==FALSE){ left <- runsum[i-1] right <- runsum[i+1]-(run[i+1]-1) slope <- (dat$Close[dat$seq==right]-dat$Close[dat$seq==left])/(run[i]+1) dat$fit[(left+1):(right-1)] <- seq(1:run[i])*slope+dat$Close[left] } return(dat) } Error in length(runsum) : object 'runsum' not found Error in length(runsum) : object 'runsum' not found [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [ESS] R window: Text is read only
If you activate the debugger (M-x toggle-debug-on-error), does it show anything meaningful? >> On Wed, Sep 27 2017 09:47, Christian wrote: > Hi, after start with M-x R the R window behaves OK, allowing writing to. BUT, > after navigating in the help system, R refuses writing with "Text is read > only". This makes R unusable. I found no help on the net, so, please, give me > some hints. > This behaviour started bothering me out of the blue. Some mis-manipulation? > TIA C. __ ESS-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/ess-help
[R] MSBVAR Package
dear sirs or madam, As I'm interested to search about the monetary transmission channel in our country by MSVAR model, I would be grateful if you help me and tell me how can I run MSVAR in R or send me the related code to run this model . Actually, I'm new user of R and I don't know how to run Markov Switching Var Model in R. Thank you very much in advance for your help. Best Regards, Ahmad Arabian [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[ESS] R window: Text is read only
Hi, after start with M-x R the R window behaves OK, allowing writing to. BUT, after navigating in the help system, R refuses writing with "Text is read only". This makes R unusable. I found no help on the net, so, please, give me some hints. This behaviour started bothering me out of the blue. Some mis-manipulation? TIA C. -- Christian Hoffmann Rigiblickstrasse 15b CH-8915 Hausen am Albis Switzerland Telefon +41-(0)44-7640853 __ ESS-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/ess-help
Re: [R] disturbed legend in ggplot2
Hi Ulrik - thanks a lot for pointing out this blunder - now it is fine! Best wishes Troels Den 27-09-2017 kl. 09:31 skrev Ulrik Stervbo: Hi Troels, Try to move the size argument out of the aesthetic. Best wishes, Ulrik On Mi., 27. Sep. 2017, 08:51 Troels Ringwrote: Dear friends - below is a subset of a much larger material showing two ways of generating two "lines". The intention is to have the colour reflect a variable, pH, but the legend is disturbed. The little part marked "3" above the colour scale is unwelcome. Why did it appear? How could I avoid it? I'm on Windows 7, R version 3.4.1 (2017-06-30) -- "Single Candle" All best wishes Troels Ring Aalborg, Denmark library(ggplot2) DF1 <- structure(list(P = c(0, 0.00222, 0.00444, 0.00667, 0.00889, 0.0111, 0.0133, 0.0156, 0.0178, 0.02, 0, 0.00222, 0.00444, 0.00667, 0.00889, 0.0111, 0.0133, 0.0156, 0.0178, 0.02), pH = c(12.3979595548431, 12.3129161148613, 12.2070984076445, 12.0669463736967, 11.8586790792785, 11.4437319273717, 7.64497330556925, 6.98905682614207, 6.63520883742788, 6.3229313658492, 12.176061323132, 12.0234712172719, 11.7861230637902, 11.2219147985144, 7.14240749824074, 6.53119941380901, 5.95522932117427, 3.25184520894594, 2.55614400932465, 2.30097494287507), BC = c(0.0576574111386315, 0.047331331067055, 0.037206026657832, 0.0268607893098731, 0.0166183791472022, 0.00639593998967551, 0.0033597279094, 0.00854377959176608, 0.00987464693654883, 0.00863636089604445, 0.0343718830720469, 0.0242985554593397, 0.0140710602077036, 0.00383913993097999, 0.00439784065436743, 0.00582135949288444, 0.00336240952299985, 0.00129948001017736, 0.00640073762860721, 0.0115158433720248), SID = c(25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15)), .Names = c("P", "pH", "BC", "SID"), row.names = c(NA, -20L), class = "data.frame") df1 <- subset(DF1,SID==25) df2 <- subset(DF1,SID==15) v <- ggplot() v <- v + geom_line(data=df1, aes(x=P, y=BC,col=pH,size=3)) v1 <- v + geom_line(data=df2, aes(x=P, y=BC,col=pH,size=3)) v <- ggplot() v <- v + geom_line(data=DF1, aes(x=P, y=BC,group=SID,col=pH,size=3)) __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] disturbed legend in ggplot2
Hi Troels, Try to move the size argument out of the aesthetic. Best wishes, Ulrik On Mi., 27. Sep. 2017, 08:51 Troels Ringwrote: > Dear friends - below is a subset of a much larger material showing two > ways of generating two "lines". The intention is to have the colour > reflect a variable, pH, but the legend is disturbed. The little part > marked "3" above the colour scale is unwelcome. Why did it appear? How > could I avoid it? > > I'm on Windows 7, R version 3.4.1 (2017-06-30) -- "Single Candle" > > All best wishes > > Troels Ring > Aalborg, Denmark > > library(ggplot2) > DF1 <- > structure(list(P = c(0, 0.00222, 0.00444, > 0.00667, 0.00889, 0.0111, > 0.0133, 0.0156, 0.0178, 0.02, > 0, 0.00222, 0.00444, 0.00667, > 0.00889, 0.0111, 0.0133, > 0.0156, 0.0178, 0.02), pH = c(12.3979595548431, > 12.3129161148613, 12.2070984076445, 12.0669463736967, 11.8586790792785, > 11.4437319273717, 7.64497330556925, 6.98905682614207, 6.63520883742788, > 6.3229313658492, 12.176061323132, 12.0234712172719, 11.7861230637902, > 11.2219147985144, 7.14240749824074, 6.53119941380901, 5.95522932117427, > 3.25184520894594, 2.55614400932465, 2.30097494287507), BC = > c(0.0576574111386315, > 0.047331331067055, 0.037206026657832, 0.0268607893098731, > 0.0166183791472022, > 0.00639593998967551, 0.0033597279094, 0.00854377959176608, > 0.00987464693654883, 0.00863636089604445, 0.0343718830720469, > 0.0242985554593397, 0.0140710602077036, 0.00383913993097999, > 0.00439784065436743, 0.00582135949288444, 0.00336240952299985, > 0.00129948001017736, 0.00640073762860721, 0.0115158433720248), > SID = c(25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 15, 15, 15, > 15, 15, 15, 15, 15, 15, 15)), .Names = c("P", "pH", "BC", > "SID"), row.names = c(NA, -20L), class = "data.frame") > > df1 <- subset(DF1,SID==25) > df2 <- subset(DF1,SID==15) > v <- ggplot() > v <- v + geom_line(data=df1, aes(x=P, y=BC,col=pH,size=3)) > v1 <- v + geom_line(data=df2, aes(x=P, y=BC,col=pH,size=3)) > > v <- ggplot() > v <- v + geom_line(data=DF1, aes(x=P, y=BC,group=SID,col=pH,size=3)) > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] disturbed legend in ggplot2
Dear friends - below is a subset of a much larger material showing two ways of generating two "lines". The intention is to have the colour reflect a variable, pH, but the legend is disturbed. The little part marked "3" above the colour scale is unwelcome. Why did it appear? How could I avoid it? I'm on Windows 7, R version 3.4.1 (2017-06-30) -- "Single Candle" All best wishes Troels Ring Aalborg, Denmark library(ggplot2) DF1 <- structure(list(P = c(0, 0.00222, 0.00444, 0.00667, 0.00889, 0.0111, 0.0133, 0.0156, 0.0178, 0.02, 0, 0.00222, 0.00444, 0.00667, 0.00889, 0.0111, 0.0133, 0.0156, 0.0178, 0.02), pH = c(12.3979595548431, 12.3129161148613, 12.2070984076445, 12.0669463736967, 11.8586790792785, 11.4437319273717, 7.64497330556925, 6.98905682614207, 6.63520883742788, 6.3229313658492, 12.176061323132, 12.0234712172719, 11.7861230637902, 11.2219147985144, 7.14240749824074, 6.53119941380901, 5.95522932117427, 3.25184520894594, 2.55614400932465, 2.30097494287507), BC = c(0.0576574111386315, 0.047331331067055, 0.037206026657832, 0.0268607893098731, 0.0166183791472022, 0.00639593998967551, 0.0033597279094, 0.00854377959176608, 0.00987464693654883, 0.00863636089604445, 0.0343718830720469, 0.0242985554593397, 0.0140710602077036, 0.00383913993097999, 0.00439784065436743, 0.00582135949288444, 0.00336240952299985, 0.00129948001017736, 0.00640073762860721, 0.0115158433720248), SID = c(25, 25, 25, 25, 25, 25, 25, 25, 25, 25, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15)), .Names = c("P", "pH", "BC", "SID"), row.names = c(NA, -20L), class = "data.frame") df1 <- subset(DF1,SID==25) df2 <- subset(DF1,SID==15) v <- ggplot() v <- v + geom_line(data=df1, aes(x=P, y=BC,col=pH,size=3)) v1 <- v + geom_line(data=df2, aes(x=P, y=BC,col=pH,size=3)) v <- ggplot() v <- v + geom_line(data=DF1, aes(x=P, y=BC,group=SID,col=pH,size=3)) __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.