Message-ID: <BANLkTinVWYgNYH964Ufr4ZaimHQ_NHx5gQ@mail.gmail.com>
Date: 2011-05-24T03:52:29Z
From: Steve Lianoglou
Subject: histogram with density
In-Reply-To: <BANLkTi=B2SqB36SqAavF60DSATK+T=Ychg@mail.gmail.com>
Hi,
On Mon, May 23, 2011 at 11:41 PM, Rekha <chithralekha86 at gmail.com> wrote:
> Hello All,*
>
> *I want to draw a histogram with density curve. *
>
> *For that simply i created a data called*"x" *and i used the function called
> * ? ?hist(x, col = "blue", freq = FALSE),** ? *from this function i got a
> histogram*.
>
> *After that , i tried this function* ?** lines(density(x), col = "red", lwd
> = 2 ?) ?*but i could not get the density curve.
>
> So, Again i used ?other function* ? ?curve(dnorm, col = 2, add = TRUE)**
> *instead of the previous function for density curve.
>
> but still am getting error in R ?like this:
>
> ?Error in plot.xy(xy.coords(x, y), type = type, ...)
> ?plot.new has not been called yet"
>
>
> * ?I have used below codes to have a data and histogram with density curve"
>
> x <- rnorm(200)
> hist(x, col = "blue", freq = FALSE)
> lines(density(x), col = "red", lwd = 2)
You should use the code above.
I don't get it though -- this is working fine for me. Are you getting
the "plot.new" error for both versions of your code?
Although I feel like in 99% of the time this won't be necessary, but
what happens if you call `plot.new()` before your call to hist?
eg:
R> x <- rnorm(200)
R> plot.new()
R> hist(x, col = "blue", freq = FALSE)
R> lines(density(x), col = "red", lwd = 2)
Does that work? If not -- what if you use `dev.new()` in place of `plot.new()`?
-steve
--
Steve Lianoglou
Graduate Student: Computational Systems Biology
?| Memorial Sloan-Kettering Cancer Center
?| Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact