graphics.off() trellis.device(color=F) xyplot(power~delta|which, data=power.df, strip= function(...) strip.default(..., style = 1), panel=function(x,y,subscripts=T){ ind<-max(subscripts)/401 indModel<-ceiling(ind/4) indPar<-ind%%4 indPar[indPar==0]<-4 cat(c(indModel,indPar), fill=T) panel.xyplot(x,y,type="l") segments(ds,tbeta[indModel,,indPar]-1.96*tsd[indModel,,indPar],ds,tbeta[indModel,,indPar]+1.96*tsd[indModel,,indPar],lwd=3) }, layout = c(3,4), xlab="delta", ylab="power")