mergeEOF {clim.pact}R Documentation

Merge EOFs

Description

Merges two EOF objects by making the first one resembling the second (a bit similar to mergeStation). Can be used for reconstructing fields.

Usage

mergeEOF(eof1,eof2,plot=TRUE,silent=FALSE,method="lm",
                match="time",cut.off=8,adjust=TRUE)

Arguments

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.

Value

An 'eof' object.

Author(s)

R.E. Benestad

Examples

## 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)

[Package clim.pact version 2.2-15 Index]