library(funcharts)
data <- simulate_data_FRTM(n_obs = 100)
X <- sapply(1:100, function(ii)
data$x_true[[ii]])
x_fd <-
fda::smooth.basis(y = X,
argvals = data$grid,
fda::create.bspline.basis(c(0, 1), 30))$fd
H <- sapply(1:100, function(ii)
data$h[[ii]])
h_fd <-
fda::smooth.basis(y = H,
argvals = data$grid,
fda::create.bspline.basis(c(0, 1), 30))$fd
mod_mFPCA <- mFPCA(x_fd, h_fd, ncom = "ptv", par_ncom = 0.95)
plot(mod_mFPCA)
Run the code above in your browser using DataLab