Code: Utilizes PBS Mapping
plots<-function(){
library(PBSmapping)
par(mfrow=c(2,1))
par(omi=c(.5,.5,.5,.5))
par(mai=c(0,0.5,0.0,.7))
par(ask=T)
nm<-c("Albacore Tuna","Bigeye Tuna","Yellowfin Tuna","Swordfish","Blue Marlin","Striped Marlin","Black Marlin","White Marlin","Bluefin Tuna","Southern Bluefin Tuna")
nm2<-c("ALB","BET","YFT","SWO","BUM","STM","BLM","WHM","BFT","SBT")
data(indLLhigh)
data(recoutglo30yr)
dfile<-recoutglo30yr
for(i in 1:levels(nm2)){
fl<-dfile[dfile$Sp==nm2[i],]
species<-nm[i]
nrec=length(fl$X)
Z<-vector(length=length(fl$X))
for (j in 1:nrec){
xx<-fl$X[j]
if(xx<20){
Z[j]=360+fl$X[j]
}else{
Z[j]=fl$X[j]
}
}
fl$X<-Z
makeplots(indLLhigh,fl,species,c(20,380),1)
makeplots(indLLhigh,fl,species,c(20,380),2)
outfile<-paste(nm2[i],"RN.eps","")
dev.copy2eps(file=outfile)
}
}
makeplots<-function(mfile,dfile,ttl,xl,zvariable){
yl<-c(-70,80)
if (zvariable==1) {
tlabel<-paste(ttl,"Recruitment",sep=" ")
} else {
tlabel<-paste(ttl,"Numbers",sep=" ")
}
plotMap(mfile,col="tan",bg="transparent",xlim=xl,ylim=yl,
main=tlabel,plt=NULL,type="n",projection="LL")
mtext("(Thousands)",side=3,line=0.5,cex=0.7)
grid<-makeGrid(x=seq(xl[1],xl[2],5),y=seq(yl[1],yl[2],5),projection="LL")
events<-dfile
if (zvariable==1) {
events$Z<-events$Rec
} else {
events$Z<-events$Ninit
}
locData<-findPolys(events,grid)
pdata<-combineEvents(events,locData,FUN=mean)
brksmax<-max(events$Z)
brksmin<-min(events$Z)
nbrks<-20
brks<-seq(brksmin,brksmax,by=(brksmax-brksmin)/nbrks)
#cols<-rev(grey(seq(0,1,by=1/nbrks)))
cols=rev(rainbow(nbrks,start=0,end=0.6))
cols[1]<-"darkblue"
pdata<-makeProps(pdata,brks,"col",cols)
addPolys(grid,polyProps=pdata,border=F)
addPolys(mfile,col="tan")
abline(h=0,lty="dashed",col="red") abline(h=40,lty="dashed",col="red") abline(h=-40,lty="dashed",col="red")
ltext<-vector(length=length(brks))
ltext<-NA
ltext[1]<-as.integer(brks[1]/1000);ltext[nbrks]<-as.integer(brks[nbrks]/1000);ltext[as.integer(nbrks/2)]<-as.integer(brks[as.integer(nbrks/2)]/1000)
par(xpd=NA)
tmp<-legend(xl[2]+.1,yl[2]-.1,legend=rep("",nbrks),col=rev(cols),pch=15,pt.cex=1.2,bty="n",y.intersp=.45)
text(tmp$rect$left+10,tmp$text$y,pos=4,rev(ltext),cex=.7)
par(xpd=FALSE)
}
-- _____________________________ Robert Ahrens Department of Zoology University of British Columbia 6270 University Blvd. Vancouver, B.C. V6T 1Z4 604.822.0046 [EMAIL PROTECTED]
______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
