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

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))