setwd("E:/R/R-beginer-guide/data/RBook")Env <- read.table(file="RIKZENV.txt",header=TRUE)Env$MyTime <- Env$Year+Env$dDay3/365

library(lattice)

xyplot(SAL ~ MyTime | factor(Station),type="1",strip = function(bg,...)strip.default(bg='white',...),col.line=1,data=Env)xyplot(SAL ~ MyTime | factor(Station),data=Env)xyplot(SAL ~ MyTime | factor(Station),type="1",strip=TRUE,col.line = 1,data=Env)xyplot(SAL ~ MyTime | factor(Station),type="1",strip=FALSE,col.line = 1,data=Env)

bwplot(SAL ~ factor(Month)|Area ,strip=strip.custom(bg='white'),cex=0.5,layout=c(2,5),data=Env,xlab="Month",ylab="Salinity",par.settings=list(box.rectangle = list(col=1),box.umbrella = list(col = 1),box.symbol = list(cex = 0.5,col =1)))

dotplot(factor(Month)~SAL|Station,subset= Area=="OS",jitter.x=TRUE,col=1,data=Env,strip=strip.custom(bg='white'),cex=0.5,ylab="Month",xlab="Salinity")

histogram(~ SAL|Station,data=Env,subset=(Area == "OS"),layout=c(1,4),nint=30,xlab="Salinity",ylab="Frequencies",strip=FALSE,strip.left=TRUE)

xyplot(SAL~Month|Year,data=Env,type=c("p"),subset=(Station=="GROO"),xlim=c(0,12),ylim=c(0,30),pch=19,panel=function(...){panel.xyplot(...)panel.grid(...,h=-1,v=-1)panel.loess(...)})

dotplot(factor(Month)~SAL|Station,pch=16,subset=(Area=="OS"),data=Env,ylab="Month",xlab="Salinity",panel=function(x,y,...){Q<- quantile(x,c(0.25,0.5,0.75),na.rm=TRUE)R<- Q[3]-Q[1]L<- Q[2]-3*(Q[3]-Q[1])MyCex <- rep(0.4,length(y))MyCol<- rep(1,length(x))MyCex[x<L]<-1.5MyCol[x<L]<-2panel.dotplot(x,y,cex=MyCex,col=MyCol,...)})

setwd("E:/R/R-beginer-guide/data/RBook")Sparrows <- read.table(file="Sparrows.txt",header=TRUE)xyplot(Wingcrd~Tarsus|Species*Sex,xlab="Axis 1",ylab="Axis 2",data=Sparrows,xlim=c(-1.1,1.1),ylim=c(-1.1,1.1),panel=function(subscripts,...){zi <- Sparrows[subscripts,3:8]di <- princomp(zi,cor=TRUE)Load <- di$loadings[,1:2]Scor <- di$scores[,1:2]panel.abline(a=0,b=0,lty=0,col=1)panel.abline(h=0,v=0,lty=2,col=1)for(i in 1:6){llines(c(0,Load[i,1]),c(0,Load[i,2]),col=1,lwd=2)ltext(Load[i,1],Load[i,2],rownames(Load)[i],cex=0.7)}sc.max<-max(abs(Scor))Scor <- Scor/sc.maxpanel.points(Scor[,1],Scor[,2],pch=1,cex=0.5,col=1)})