mergeEOF              package:clim.pact              R Documentation

_M_e_r_g_e _E_O_F_s

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

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

_A_r_g_u_m_e_n_t_s:

    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.

_V_a_l_u_e:

     An 'eof' object.

_A_u_t_h_o_r(_s):

     R.E. Benestad

_E_x_a_m_p_l_e_s:

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

