Re: [R] Plotting rows (or columns) from a matrix in different graphs, not using "par"

From: Gabor Grothendieck <ggrothendieck_at_gmail.com>
Date: Wed 15 Jun 2005 - 14:02:22 EST

On 6/14/05, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> On 6/14/05, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> > On 6/14/05, Camarda, Carlo Giovanni <Camarda@demogr.mpg.de> wrote:
> > > Dear R-users,
> > > I would like to ask whether it's possible (for sure it would be), to
> > > plot each rows (or columns) in different graphs and in the same figure
> > > region without using the function "par" and then struggling around with
> > > "axes" and labels etc.
> > > Luckily, I would always have "rows + columns = even number" and the same
> > > "ylim".
> > >
> > > The next one could be a sort of example on what I would like to avoid
> > > plotting the rows of the matrix "mat":
> > >
> > > ########### EXAMPLE ######################
> > > dat <- sort(runif(16, 1, 1000))
> > > mat <- matrix(dat, ncol=4, byrow = T)
> > > y <- seq(1950, 1953)
> > > par(mfrow=c(2,2))
> > > plot(y, mat[1,], ylim=c(1,1000), axes=F, xlab="", ylab="simul.")
> > > box()
> > > axis(side=2, tick=T, labels=T)
> > > axis(side=3, tick=T, labels=T)
> > > plot(y, mat[2,], ylim=c(1,1000), axes=F, xlab="", ylab="")
> > > box()
> > > axis(side=3, tick=T, labels=T)
> > > plot(y, mat[3,], ylim=c(1,1000), axes=F, xlab="years", ylab="simul.")
> > > box()
> > > axis(side=1, tick=T, labels=T)
> > > axis(side=2, tick=T, labels=T)
> > > plot(y, mat[4,], ylim=c(1,1000), axes=F, xlab="years", ylab="")
> > > box()
> > > axis(side=1, tick=T, labels=T)
> > > ########### END EXAMPLE ########################
> > >
> > > Naturally something more compact would be even nicer.
> > >
> >
> > plot(ts(mat, start = 1950), nc = 2)
> >
>
> I just looked at this again and noticed that it did not have
> identical y axes. I tried with ylim= too but plot.ts seemed
> to ignore it.
>
> Here is another solution that uses the zoo library.
> Be sure to use zoo 1.0-1. It uses the screens= argument of
> plot.zoo which was not available in earlier versions of zoo.
>
> Instead of creating a ts series with 4 columns create a
> similar zoo series. Append to that 4 dummy series all
> with the same range. Use the screens= argument
> on plot.zoo to plot one real series and one dummy series
> in each graph. Use type = "l" to plot the real series
> and type = "n" for the dummy series. This causes the
> real series to be plotted with lines and the dummy
> series to be plotted invisibly yet still affect the
> y axis range. Note that we have used the fact that
> mat is square so if the real mat is not square be sure
> to modify this accordingly.
>
> library(zoo) #
>
> n <- nrow(mat)
> z <- zooreg(cbind(mat, max(mat) * diag(n)), start = 2000)
> plot(z, screens = c(1:4, 1:4), type = rep(c("l","n"), each = 4), nc = 2)
>

The problem with ylim= that seems to affect both plot.ts and plot.zoo has been fixed in zoo so in the next version of zoo the above will be reducable to just this:

   plot(zooreg(mat, start = 2000), ylim = range(mat), nc = 2)



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jun 15 14:05:09 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:39 EST