功能主成分分数的双图答案

【问题标题】：Biplots for Functional Principal Component Scores功能主成分分数的双图
【发布时间】：2019-09-28 19:39:05
【问题描述】：

我正在尝试获取两个 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 实现此功能，请不要犹豫，写一个答案。

【问题讨论】：

标签： r pca biplot

【解决方案1】：

对于这种情况，biplot 不是正确的函数，因为它不承认 fd 和 pcafd objects。 FDA 软件包有plotscores，它允许对主成分分数进行 BIplot，或者换句话说，在主成分对上绘制曲线的分数

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

【讨论】：