[R] trellis graphics in loops

Chuck Cleland ccleland at optonline.net
Sun Jan 30 17:21:04 CET 2005


Wrap the xyplot in print:

X <- data.frame(x=rnorm(10000), y=rnorm(10000), z=sample(c("foo1","foo2"),
10000,replace=T), year=sample(c(89:94), 10000,replace=T))

trellis.device("postscript", file="foo.eps")

for(i in unique(X$year)){
   temp <- X[X$year%in%i, ]
   print(xyplot(temp$x~temp$y|temp$z))
}

dev.off()

hope this helps,

Chuck Cleland

Jean Eid wrote:
> I have this awkward problem with trellis (lattice). I am trying to
> generate some plots through loops but the .eps file is empty. When I
> generate them in a list and print them outside the loop all is fine. this
> is an example below:( nothing shows up in foo.eps, but all show up in
> foo1.eps)
> R vesion 2.0.1, lattice version  0.10-16, on a debian 2.6.8-1 kernel.
> 
> 
> X <- data.frame(x=rnorm(10000), y=rnorm(10000), z=sample(c("foo1","foo2"),
> 10000,replace=T), year=sample(c(89:94), 10000,replace=T))
> trellis.device("postscript", file="foo.eps")
> for(i in unique(X$year)){
>   temp <- X[X$year%in%i, ]
>   xyplot(temp$x~temp$y|temp$z)
> }
> dev.off()
> 
> mylist <- list()
> for(i in unique(X$year)){
>   temp <- X[X$year%in%i, ]
>   mylist[[i]] <- xyplot(temp$x~temp$y|temp$z)
> }
> 
> trellis.device("postscript", file="foo1.eps")
> mylist[[89]]
> mylist[[90]]
> mylist[[91]]
> mylist[[92]]
> mylist[[93]]
> mylist[[94]]
> dev.off()
> 
> 
> Thank you,
> 
> Jean
> 
> ______________________________________________
> R-help at 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
> 

-- 
Chuck Cleland, Ph.D.
NDRI, Inc.
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212) 845-4495 (Tu, Th)
tel: (732) 452-1424 (M, W, F)
fax: (917) 438-0894




More information about the R-help mailing list