| mergeEOF {clim.pact} | R Documentation |
Merges two EOF objects by making the first one resembling
the second (a bit similar to mergeStation). Can be used
for reconstructing fields.
mergeEOF(eof1,eof2,plot=TRUE,silent=FALSE,method="lm",
match="time",cut.off=8,adjust=TRUE)
eof1 |
The EOF object (see EOF) for the earliest data (assuming this
data set contains most errors). |
eof2 |
The EOF object (see EOF) for the most recent data (assuming this
data set contains less errors) |
plot |
if 'TRUE' then plot diagnostics |
silent |
if 'TRUE' then do not print out diagnostics. |
method |
"lm" or "matrix projection" (G. Strang (1988), "Linear algrebra and its applications", Hartcourt Brace & Company, 3rd ed.(p.147)) |
match |
"time" or "space" |
cut.off |
How many EOFs to include |
adjust |
TRUE: apply adjust.eof to results. |
An 'eof' object.
R.E. Benestad
## Not run:
data(DNMI.slp)
NCEP.slp<-retrieve.nc("~/data/ncep/slp.mon.mean.nc",
x.rng=c(-90,50),y.rng=c(0,75))
#--------------------------------------------------
# Need to fix some details of the NCEP.slp object (only for
# proper 'housekeeping')
NCEP.slp$dd[] <- 15
attr(NCEP.slp$tim,unit') <- "month"
class(NCEP.slp) <- c("field","monthly.field.object")
#---------------------------------------------------
eof1 <- EOF(DNMI.slp,mon=1)
eof2 <- EOF(NCEP.slp,mon=1)
eof <- mergeEOF(eof1,eof2)
## End(Not run)