先说一下你构想的存矩阵再循环方法,读入整个csv,形成矩阵A,之后a1<-A[1:53,],样本量不多的话,人工多做几次就好,样本量大的话,便一个循环,for (i in 1:((dim(A)[1])%%53+1)), ai<-A[((i*53)-52):(i*53)]。
不用你的构想,直接编个函数 function (x) B[x,]<-A[x+53,]-A[x,]
然后,用apply 族函数对A矩阵执行所编的函数,应该可以得到新的矩阵B。也可能在执行前需要设个空的矩阵B。
只是想法阿,细节可能会报错,我在外面没法验证,你自己调整代码吧