plotANM {anm} | R Documentation |
Three plots are provided: a plot of the minimum distances versus time, a plot comparing the analogs from EOF and observations and a plot of errors versus time.
plotANM(x,tmp,station,eof.file,leps)
x |
the anm object inheriting from anm routine and for which prediction is desired. |
tmp |
True if the analysis is on temperature, False if on precipitation. |
station |
the name of the station. |
eof.file |
string giving the name of the eof file used for the study. |
leps |
if true, postscripts are created for the plots. |
Alexandra Imbert
anm, stepANM, predictAnm, printAnm
library(survival) library(clim.pact) data(temp.era) data(susendal) y<-susendal$V6 # temperatures X<- eof$PC[,c(1,2)] calibration <- c(susendal$V4>1979 & susendal$V4<1990 & (susendal$V3==1 | susendal$V3==2 | susendal$V3==12)) evaluation <- c((susendal$V4>1990 & susendal$V4<1993 | susendal$V4==1990) & (susendal$V3==1 | susendal$V3==2 | susendal$V3==12)) y.calib <- y[calibration] y.eval <- y[evaluation] eof.calib <- c(eof$yy>1979 & eof$yy<1990) eof.eval <- c(eof$yy> 1990 & eof$yy<1993| eof$yy==1990) period <- c(calibration, evaluation) y.period <- y[(susendal$V4>1979 & susendal$V4<1993) & (susendal$V3==1 | susendal$V3==2 | susendal$V3==12)] test.data <-data.frame(y=y.period, X1=X[eof$yy< 1993 & eof$yy> 1979,1], X2=X[eof$yy< 1993 & eof$yy> 1979,2], yy=eof$yy[eof.calib | eof.eval], mm=eof$mm[eof.calib | eof.eval], dd=eof$dd[eof.calib | eof.eval]) test.anm<-anm(y ~ X1 + X2,data=test.data) plotANM(test.anm,TRUE,"Susendal","eof_ERA-15_TEM_16E31E-64N73N_DJF",FALSE)