功能主成分分数的双标图

Question

我正在尝试获取两个 F. 主成分（或谐波）之间的双标图。我提供了 fda 包文档中的示例。解谜：

library(fda)
#BASIS FUNCTIONS
daybasis65 <- create.fourier.basis(c(0, 365), nbasis=65, period=365)
#DATA
harmaccelLfd <- vec2Lfd(c(0,(2*pi/365)^2,0), c(0, 365))
#SMOOTHING BASIS TO THE DATA
harmfdPar <- fdPar(daybasis65, harmaccelLfd, lambda=1e5)
daytempfd <- smooth.basis(day.5, CanadianWeather$dailyAv[,,"Temperature.C"],
daybasis65, fdnames=list("Day", "Station", "Deg C"))$fd
#FPCA
daytemppcaobj <- pca.fd(daytempfd, nharm=4, harmfdPar)
#RETRIEVING COMPONENTS
f_loadings <- daytemppcaobj $harmonics
#BIPLOT BETWEEN 1ST AND 2ND HARMONIC
biplot(f_loadings[1],f_loadings[2])
+Error in 1L:n : argument of length 0

请注意 f_loadings 和 daytempfd 是 fd 对象，而 daytemppcaobj 是 pca.fd 对象。

提前谢谢你。欢迎并奖励所有类型的帮助！

如果您有任何关于使用 ggplot2 实现此功能的想法，请不要犹豫，写下答案。

Answer 1

对于这种情况，biplot 不是正确的函数，因为它不接受 fd 和 pcafd objects。 FDA 包有 plotscores 允许 BIplot Principal Component Scores 或者换句话说，在成对的 principal components

上绘制曲线分数

plotscores(daytemppcaobj, scores=c(1, 2))

功能主成分分数的双标图

Biplots for Functional Principal Component Scores

r

pca

biplot