McClatchie, Sam (PIRSA-SARDI)
2005-Mar-16 05:31 UTC
[R] Fixed/ trellis/ panel.superpose/ passing in superscripts/ groups
Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues I have solved the miscoding with subscript and groups in panel.superpose. I was neglecting to pass in the variables correctly. "sardine.egg.T.S.space.2001.and.2002.exp" <- function() { library(lattice) library(grid) ## trellis.device(postscript, ## file="../figures/sardine.egg.T.S.space.2001.and.2002.ps", ## horizontal=FALSE, color=TRUE) year <- as.factor(rep(2001, dim(mn.ts.e.2001)[1])) year.2001 <- cbind(mn.ts.e.2001,year) year <- as.factor(rep(2002, dim(mn.ts.e.2002)[1])) year.2002 <- cbind(mn.ts.e.2002,year) mn.ts.e.both <- rbind(year.2001, year.2002) #browser() ###trellis plot ## define conditioning interval int <- matrix(c(0,2,3,4,5,8,9,16,17,32,33,64), ncol=2, byrow=TRUE) egg.counts <- shingle(mn.ts.e.both$eggs2.Pilch.Eggs, intervals = int) ##larvae.counts <- shingle(mn.ts.e.both$eggs2.Pilch.Larv, intervals int) ## define plot structure out1 <- xyplot(mn.t ~ mn.s | egg.counts, data = mn.ts.e.both, groups = year, xlim = c(35,38), ylim = c(12,24), xlab = "mean salinity", ylab = "mean temperature (deg. C)", main = "2001 and 2002 egg densities in Temperature-Salinity space", aspect = "xy", jitter = T, layout = c(1,6), auto.key=TRUE, ##add to individual panels panel = function(x, y, subscripts, groups){ ## plot all the T-S data for both years ## using explicit variables, differing from x,y panel.xyplot(data.2001$Salinity, data.2001$Temperature.oC, pch=".", col="yellow"); panel.xyplot(data.2002$Salinity, data.2002$Temperature.oC, pch=".", col="orange"); ## add some labels for watermasses panel.text(37.8, 15, "Upwelling"); panel.text(37.8, 17.5, "Warm pool"); panel.text(37.8, 19.75, "Shelf break"); panel.text(37.8, 21.5, "Spencer Gulf"); panel.abline(h = c(16,19,20.5), v = 36.5, col="red", lty=2); ## overlay the conditioned egg data on the TS plots panel.superpose(x, y, subscripts, groups); ## add a legend auto.key = T } ) print(out1) ## graphics.off() } ---- Sam McClatchie, Biological oceanography South Australian Aquatic Sciences Centre PO Box 120, Henley Beach 5022 Adelaide, South Australia email <mcclatchie.sam at saugov.sa.gov.au> Telephone: (61-8) 8207 5448 FAX: (61-8) 8200 2481 Research home page <http://www.members.iinet.net.au/~s.mcclatchie/> /\ ...>><xX(?> //// \\\\ <?)Xx><< ///// \\\\\\ ><(((?> >><(((?> ...>><xX(?>O<?)Xx><<