predict.MFPCAfit: Function prediction based on MFPCA results

In ClaraHapp/MFPCA: Multivariate Functional Principal Component Analysis for Data Observed on Different Dimensional Domains

Description

Predict functions based on a truncated multivariate Karhunen-Loeve representation:

\hat x = \hat mu + ∑_{m = 1}^M ρ_m \hat ψ_m

with estimated mean function \hat μ and principal components ψ_m. The scores ρ_m can be either estimated (reconstruction of observed functions) or user-defined (construction of new functions).

Usage

## S3 method for class 'MFPCAfit'
predict(object, scores = object$scores, ...)

Arguments

object	An object of class MFPCAfit, typically resulting from a MFPCA function call.
scores	A matrix containing the score values. The number of columns in scores must equal the number of principal components in object. Each row represents one curve. Defaults to the estimated scores in object, which yields reconstructions of the original data used for the MFPCA calculation.
...	Arguments passed to or from other methods.

Value

A multiFunData object containing the predicted functions.

See Also

MFPCA

Examples

#' # Simulate multivariate functional data on one-dimensonal domains
# and calculate MFPCA (cf. MFPCA help)
set.seed(1)
# simulate data (one-dimensional domains)
sim <-  simMultiFunData(type = "split", argvals = list(seq(0,1,0.01), seq(-0.5,0.5,0.02)),
                       M = 5, eFunType = "Poly", eValType = "linear", N = 100)
# MFPCA based on univariate FPCA
PCA <- MFPCA(sim$simData, M = 5, uniExpansions = list(list(type = "uFPCA"),
                                                     list(type = "uFPCA")))

# Reconstruct the original data
pred <- predict(PCA) # default reconstructs data used for the MFPCA fit

# plot the results: 1st element
plot(sim$simData[[1]]) # original data
plot(pred[[1]], add = TRUE, lty = 2) # reconstruction

# plot the results: 2nd element
plot(sim$simData[[2]]) # original data
plot(pred[[2]], add = TRUE, lty = 2) # reconstruction

ClaraHapp/MFPCA documentation built on Nov. 9, 2021, 11:35 p.m.