| printAnm {anm} | R Documentation |
Prints the coefficients of an anm object.
printAnm(x, digits = max(3, getOption("digits") - 3), ...)
x |
the anm object. |
digits |
the vector defining the format of printing. |
... |
currently disregarded. |
Alexandra Imbert
library(survival)
library(clim.pact)
data(susendal)
data(temp.era)
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)
printAnm(test.anm)