[R] curve
Sarah Goslee
sarah.goslee at gmail.com
Mon Dec 13 20:01:28 CET 2010
Here's one way to do what I think you want:
test<- rnorm(5000,1000,100)
test1 <- subset(test, subset=(test > 1100))
d <- density(test)
plot(d, main="Density of production", xlab="")
lines(d$x[d$x > 1100], d$y[d$x > 1100], col="blue", lwd=2)
curveheight <- d$y[abs((d$x - mean(test1))) == min(abs((d$x - mean(test1))))]
segments(x0=mean(test1), y0=0, y1=curveheight)
Sarah
On Mon, Dec 13, 2010 at 1:44 PM, Val <valkremk at gmail.com> wrote:
> Hi All,
>
> I generated 5000 samples using the following script
>
> test<- rnorm(5000,1000,100)
> test1 <- subset(test, subset=(test > 1100))
> d <- density(test)
> plot(d, main="Density of production")
> abline(v=mean(test1)
>
> I wanted to do the following but faced difficulties
> 1. to shade or color (blue) the curve using the criterion that any values
> greater than 11,000
> 2. I drew a vertical line but I wanted the v-line within the curve not to
> stick outside the curve
> 3. to suppress the output produced at the bottom of the curve( N=5000 and
> bandwidth =16.22)
>
> Thanks in advance
> Val
>
--
Sarah Goslee
http://www.functionaldiversity.org
More information about the R-help
mailing list