我目前有一套每月时间序列数据，可以追溯到1993年，涉及16个变量。我编写了在整个时间段内运行PCA的脚本，但是事后看来，我需要在连续12个月的时间段内运行计算（因此，第一次迭代将针对t = 1到t = 12）。我目前有：
path<- "C:/Research/Datasets"
setwd(path)
variables<- read.csv("Index (1993).csv")
prin_comp<- prcomp(variables)
prin_comp$rotation[1:16, 1]
std_dev<-prin_comp$sdev
pr_var<- stddev^2
prop_varex<- pr_var/sum(pr_var)
我还想将所有组件的所有迭代和变体上的负载保存在第一个组件上，但是我对如何做这些事情感到困惑。
任何帮助将不胜感激。预先感谢！