kernel density plot
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