Болты для оценки функциональных основных компонентов

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 или, другими словами, составлять оценки кривых по парам основных компонентов

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