Hello, I am using the following codes to plot a histogram and density line for x. For the density line, I just want it to show the two tails, eg, for x larger than 0.05 ans smaller than -0.05 hist (x, seq(-0.1,0.1,0.01),freq = FALSE) lines (density(x,bw="SJ"), x> 0.05 & x< (-0.05), col = "red") But is does not work, can anyone give me some advice? -- View this message in context: http://www.nabble.com/Histogram-tf3947281.html#a11197644 Sent from the R help mailing list archive at Nabble.com.
I expect there's a more elegant way of doing this, but this should work: set.seed(101) x <- rnorm(500,sd=.03) hist (x, seq(-0.1,0.1,0.01),freq = FALSE) d <- density(x,bw="SJ") lowt <- d$x < -.05 upt <- d$x > .05 lines (d$x[lowt],d$y[lowt], col = "red") lines(d$x[upt],d$y[upt], col = "red") On 19/06/07, livia <yn19832 at msn.com> wrote:> > Hello, I am using the following codes to plot a histogram and density line > for x. For the density line, I just want it to show the two tails, eg, for x > larger than 0.05 ans smaller than -0.05 > > hist (x, seq(-0.1,0.1,0.01),freq = FALSE) > lines (density(x,bw="SJ"), x> 0.05 & x< (-0.05), col = "red") > > But is does not work, can anyone give me some advice? > > > > -- > View this message in context: http://www.nabble.com/Histogram-tf3947281.html#a11197644 > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > R-help at stat.math.ethz.ch 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. >-- ================================David Barron Said Business School University of Oxford Park End Street Oxford OX1 1HP
Your subsetting expression in lines does not make any sense at all. Not tested but maybe something like: lines (density(subset(x, x> 0.05 & x< -0.05)bw=SJ), col='red") --- livia <yn19832 at msn.com> wrote:> Hello, I am using the following codes to plot a > histogram and density line > for x. For the density line, I just want it to show > the two tails, eg, for x > larger than 0.05 ans smaller than -0.05 > > hist (x, seq(-0.1,0.1,0.01),freq = FALSE) > lines (density(x,bw="SJ"), x> 0.05 & x< (-0.05), col > = "red") > > But is does not work, can anyone give me some > advice? > > > > -- > View this message in context: >http://www.nabble.com/Histogram-tf3947281.html#a11197644> Sent from the R help mailing list archive at > Nabble.com. > > ______________________________________________ > R-help at stat.math.ethz.ch 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. >