mfpcaPlot: Plot multivariate functional pca

In Funclustering: Functional Data Clustering

Description

This function plots the functional pca.

Usage

mfpcaPlot(pca, grid = c())

Arguments

pca	is the result of mfpca. In the univariate case mfpcaPlot use the package fda and will be similar to it's function "plot.pca.fd". In multivariate functional pca, we will make a graphic window for each dimension.
grid	specify how to divide the graphics window. grid=c(n,m) divided the widow in to n lines and m columns. If user don't specify grid then he must enter <Enter> to pass to the next graphic.

Examples

# Multivariate
# ---------  CanadianWeather (data from the package fda) --------
CWtime<- 1:365
CWrange<-c(1,365)
CWbasis <- create.fourier.basis(CWrange, nbasis=65)
harmaccelLfd <- vec2Lfd(c(0,(2*pi/365)^2,0), rangeval=CWrange)

# -- Build the curves ---
temp=CanadianWeather$dailyAv[,,"Temperature.C"]
CWfd1 <- smooth.basisPar(
CWtime, CanadianWeather$dailyAv[,,"Temperature.C"],CWbasis, 
Lfdobj=harmaccelLfd, lambda=1e-2)$fd
precip=CanadianWeather$dailyAv[,,"Precipitation.mm"]
CWfd2 <- smooth.basisPar(
CWtime, CanadianWeather$dailyAv[,,"Precipitation.mm"],CWbasis, 
Lfdobj=harmaccelLfd, lambda=1e-2)$fd

CWfd=list(CWfd1,CWfd2) 

pca=mfpca(CWfd,nharm=4)
mfpcaPlot(pca,grid=c(2,2))

Funclustering documentation built on May 2, 2019, 5:05 p.m.