Hello,
Something like this?
sag <- matrix(rnorm(1e3 * 9), ncol = 9)
d.list <- apply(sag, 2, density)
xrange <- range(sapply(d.list, function(d) range(d$x)))
ymax <- max(sapply(d.list, function(d) max(d$y)))
op <- par(mfrow=c(3,3))
for(i in 1:9)
plot(d.list[[i]], xlim = xrange, ylim = c(0, ymax))
par(op)
Hope this helps,
Rui Barradas
Em 27-09-2012 01:53, Meredith Ballard LaBeau escreveu:
Good Evening-
I have a set of nine scenarios I want to plot to see how the distribution
is changing, if one tail is getting larger in certain scenario, currently I
am using this code
colnames<-dimnames(sag_pdfs)[[2]]
par(mfrow=c(3,3))
for(i in 1:9) {
d<-density(sag[,i])
plot(d,type="n", main=colnames[i])
polygon(d,col="red",border="grey")}
where sag is a 7305x9 double matrix and 9 different scenarios. I want to be
able to compare the distribution using the same axes scale.
Can anyone help?
Thanks
Meredith LaBeau
______________________________________________
R-help@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.
