You are most likely to get reasonable help if you provide at a minimum your OS, the version of R you're using, and the error message you get.
And also if you send your replies to the R-help list, not just me. Sarah On Tue, Dec 14, 2010 at 10:58 AM, Val <valkr...@gmail.com> wrote: > Thank you Sarah, > > > It worked in my Linux machine as well but not in Windows. > > > > > > On Tue, Dec 14, 2010 at 9:37 AM, Sarah Goslee <sarah.gos...@gmail.com> > wrote: >> >> Val, >> >> Here's the complete console output. The graph produced is at: >> http://www.functionaldiversity.org/temp/curve.png >> >> R version 2.12.0 (2010-10-15) >> Copyright (C) 2010 The R Foundation for Statistical Computing >> ISBN 3-900051-07-0 >> Platform: x86_64-redhat-linux-gnu (64-bit) >> >> R is free software and comes with ABSOLUTELY NO WARRANTY. >> You are welcome to redistribute it under certain conditions. >> Type 'license()' or 'licence()' for distribution details. >> >> Natural language support but running in an English locale >> >> R is a collaborative project with many contributors. >> Type 'contributors()' for more information and >> 'citation()' on how to cite R or R packages in publications. >> >> Type 'demo()' for some demos, 'help()' for on-line help, or >> 'help.start()' for an HTML browser interface to help. >> Type 'q()' to quit R. >> >> > ls() >> character(0) >> > >> > test<- rnorm(5000,1000,100) >> > test1 <- subset(test, subset=(test > 1100)) >> > d <- density(test) >> > >> > png("curve.png") >> > plot(d, main="Density of production", xlab="") >> > >> > xarea <- c(1100, d$x[d$x > 1100], c(max(test))) >> > yarea <- c(0, d$y[d$x > 1100], 0) >> > >> > >> > polygon(xarea, yarea, col="blue") >> > >> > curveheight <- d$y[abs((d$x - mean(test1))) == min(abs((d$x - >> > mean(test1))))] >> > segments(x0=mean(test1), y0=0, x1=mean(test1), y1=curveheight, >> > col="red", lwd=2) >> > dev.off() >> null device >> 1 >> > >> >> Sarah >> >> >> On Tue, Dec 14, 2010 at 8:59 AM, Val <valkr...@gmail.com> wrote: >> > Hi Sara, >> > >> > Could you please send me your output as an attached file if it is >> > possible? >> > >> > Thanks >> > >> > >> >> >> 1. to shade or color (blue) the curve using the criterion that any >> >> >> values >> >> >> greater than 11,000 >> >> > >> >> > I think I was not clear in the above point. I want shade not the line >> >> > but >> >> > the area under the curve, >> >> >> >> Here's an example of how to do that using polygon: >> >> http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=7 >> >> >> >> > and >> >> > Your last line of code, >> >> > segments(x0=mean(test1), y0=0, y1=curveheight) >> >> > >> >> > gave me the following error message >> >> > >> >> > Error in segments(x0 = mean(test1), y0 = 0, y1 = curveheight) : >> >> > element 3 is empty; >> >> > the part of the args list of '.Internal' being evaluated was: >> >> > (x0, y0, x1, y1, col = col, lty = lty, lwd = lwd, ...) >> >> > >> >> > could you check it please >> >> >> >> I checked it before I sent it to you. The code I provided works >> >> correctly >> >> on my computer. (R 2.12.0, Linux). >> >> >> >> You could try this statement instead: >> >> segments(x0 = mean(test1), y0 = 0, x1=mean(test1), y1 = curveheight) >> >> >> >> Sarah >> >> >> >> -- >> >> >> >> -- >> Sarah Goslee >> http://www.functionaldiversity.org > > -- Sarah Goslee http://www.stringpage.com http://www.astronomicum.com http://www.functionaldiversity.org ______________________________________________ R-help@r-project.org mailing list 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.
