[R] kernel density plot

[David Winsemius]

On Apr 3, 2011, at 6:56 AM, Muzna Alvi wrote:

> 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?

Why not just set ylab="" in subsequent calls to plot?

--

David Winsemius, MD
West Hartford, CT