Re: [R] kernel density plot

From: Greg Snow <Greg.Snow_at_imail.org>
Date: Sun, 03 Apr 2011 13:49:40 -0600

It would be easier to test and example if we knew what s22, s33, and s44 were, but you can try:

density(s22$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->t
plot(t, xlim=c(-30000,40000), main="Net Income Distribution", axes=F, ylim=c(0,0.00035). xlab="Value in Rupees") density(s33$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->u
lines(u, col="red")
density(s44$Net_income_Total.1, bw="nrd0",adjust=1, kernel=c("gaussian"))->v
lines(v, col="blue")

--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.snow_at_imail.org
801.408.8111

From: Muzna Alvi [mailto:muzna.alvi_at_gmail.com]
Sent: Sunday, April 03, 2011 12:52 PM
To: Greg Snow
Cc: r-help_at_r-project.org
Subject: Re: [R] kernel density plot

i am sorry greg, can you explain that with an example?
On Mon, Apr 4, 2011 at 12:11 AM, Greg Snow <Greg.Snow_at_imail.org<mailto:Greg.Snow_at_imail.org>> wrote:
It is better to replace your later calls to plot with calls to lines instead, then you don't need to use par(new=T) which as you see tends to cause problems.

--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.snow_at_imail.org<mailto:greg.snow_at_imail.org>
801.408.8111



> -----Original Message-----
> From: r-help-bounces_at_r-project.org<mailto:r-help-bounces_at_r-project.org> [mailto:r-help-bounces_at_r-<mailto:r-help-bounces_at_r->
> project.org<http://project.org>] On Behalf Of Muzna Alvi
> Sent: Sunday, April 03, 2011 4:56 AM
> To: r-help_at_r-project.org<mailto:r-help_at_r-project.org>
> Subject: [R] kernel density plot
>
> I am using the following commands for plotting kernel density for three
> kinds of crops
>
> density(s22$Net_income_Total.1, bw="nrd0",adjust=1,
> kernel=c("gaussian"))->t
> plot(t, xlim=c(-30000,40000), main="Net Income Distribution", axes=F,
> ylim=c(0,0.00035). xlab="Value in Rupees")
> par(new=T)
> density(s33$Net_income_Total.1, bw="nrd0",adjust=1,
> kernel=c("gaussian"))->u
> plot(u, xlim=c(-30000,40000), axes=F, main="", col="red",
> ylim=c(0,0.00035))
> par(new=T)
> density(s44$Net_income_Total.1, bw="nrd0",adjust=1,
> kernel=c("gaussian"))->v
> plot(v, xlim=c(-30000,40000), col="blue", axes=F, main="",
> ylim=c(0,0.00035))
>
> the problem is that in the graph that is drawn
>
> 1. the xlab gets hidden with the [N= and the bandwidth=] values
> 2. when i do par(new=T) this N and bandwidth value appears multiple
> times..overlapping each time and making the graph look untidy..
>
>
> Is there any way of making these N and Bandwidth values not appear in
> the
> graph?
>
> Thanks
>
> --
> --
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help_at_r-project.org<mailto:R-help_at_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.
-- -- [[alternative HTML version deleted]] ______________________________________________ R-help_at_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.
Received on Sun 03 Apr 2011 - 19:52:55 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 03 Apr 2011 - 20:30:26 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive