eofRecon | R Documentation |
This function reconstructs the original field from an EOF object of the function eof.
eofRecon(EOF, pcs = NULL, newpcs = NULL, uncenter = TRUE, unscale = TRUE)
EOF |
An object resulting from the function eof. |
pcs |
The principal components (PCs) to use in the reconstruction
(defaults to the full set of PCs: |
newpcs |
An (optional) matrix of new principal coordinates to use in the reconstruction. This would typically come from a gappy dataset whose missing values are to be predicted based on the EOF loadings of a the EOF object (see eofPred). |
uncenter |
Logical. Should reconstructed matrix be un-centered (e.g. if
|
unscale |
Logical. Should reconstructed matrix be un-centered (e.g. if
|
set.seed(1) iris.gappy <- as.matrix(iris[,1:4]) iris.gappy[sample(length(iris.gappy), 0.25*length(iris.gappy))] <- NaN Er <- eof(iris.gappy, method="svd", recursive=TRUE) # recursive (RSEOF) Enr <- eof(iris.gappy, method="svd", recursive=FALSE) # non-recursive (LSEOF) iris.gappy.recon.r <- eofRecon(Er) iris.gappy.recon.nr <- eofRecon(Enr) # Reconstructed values vs. observed values op <- par(mfrow=c(1,2)) lim <- range(iris.gappy, na.rm=TRUE) plot(iris.gappy, iris.gappy.recon.r, col=c(2:4)[iris$Species], main="recursive=TRUE", xlim=lim, ylim=lim) abline(0, 1, col=1, lwd=2) plot(iris.gappy, iris.gappy.recon.nr, col=c(2:4)[iris$Species], main="recursive=FALSE", xlim=lim, ylim=lim) abline(0, 1, col=1, lwd=2) par(op) # Reconstructed values from gappy data vs. all original values op <- par(mfrow=c(1,2)) plot(as.matrix(iris[,1:4]), iris.gappy.recon.r, col=c(2:4)[iris$Species], main="recursive=TRUE") abline(0, 1, col=1, lwd=2) plot(as.matrix(iris[,1:4]), iris.gappy.recon.nr, col=c(2:4)[iris$Species], main="recursive=FALSE") abline(0, 1, col=1, lwd=2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.