[R] Comparing density plots using same axes or same axes scale

[Rui Barradas]

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
>