[R] curve

[Sarah Goslee]

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