我想对长格式的数据帧进行PCA:
时间1 id1数据11
时间1 id2数据12
时间2 id1数据21
等等。
有没有一种简单的方法可以做到这一点,或者reshape
和princomp
的标准方法。我的数据集相当大,大约有40000次和4000个ID。
对于这样一个简单的重塑,我想你只需要
m <- matrix(mydata[,3],nrow=ntimes,byrow=TRUE)
princomp(m)
这将为您提供一个ntimes
乘nIDs
矩阵。它将比reshape
快(可能快很多(。