Мой вопрос о функциональном анализе главных компонент в R.
Я работаю с многомерным временным рядом, который выглядит примерно так:
Моя цель состоит в том, чтобы уменьшить размеры, применяя функциональный PCA, а затем построить первый основной компонент следующим образом:
Я уже использовал функцию FPCA пакета fdapace в наборе данных. К сожалению, я не понимаю, как интерпретировать полученную матрицу оценок FPCA (xiEst).
В моем понимании значения основных компонентов хранятся в столбцах матрицы.
К сожалению, количество столбцов не соответствует количеству временных интервалов моего многомерного временного ряда.
Я не знаю, как значения в матрице соответствуют значениям исходных данных и как построить первый основной компонент как уменьшение размеров исходных данных.
Если вам нужен какой-то код для воспроизведения ситуации, вы можете использовать набор данных medfly пакета:
library(fdapace)
data(medfly25)
Flies <- MakeFPCAInputs(medfly25$ID, medfly25$Days, medfly25$nEggs)
pfcaObjFlies <- FPCA(Flies$Ly, Flies$Lt)
когда я строю первый главный компонент через
plot(fpcaObjFlies$xiEst[,1], type = "o")
график не совсем соответствует моим ожиданиям:
Я бы ожидал график с 25 наблюдениями, похожий на графики набора данных medfly.