我得到了一个4年时间序列的资产回报,我试图执行滚动窗口,以估计方差协方差矩阵的6个月的校准期。总的来说,我应该得到40个协方差矩阵。我试着运行下面写的代码,但是它是错误的。我如何修改这个R码?data <- matrix(rnorm(100), 20, 5) #data represents the returns of 5 assets over 20 days
我想要校准
从基因表达数据(40000个基因(变量)x 30个观察值)中,我想创建一个40000 x 40000协方差矩阵。这绝对比我的内存大。使用包'ff‘,我设法为相关性预先分配了一个40000x40000的空矩阵。然而,'cov‘或'cor’函数在我的系统上只能管理一个5000x5000的协方差矩阵,所以我必须按块进行1:5000,5001:10000等协方差计算</
Im有N对投资组合权重存储在一个numpy数组中,并且希望计算投资组合风险,这是w * E * w_T,其中w_T是权转置。我想出的方法是循环每个权重对并应用矩阵乘法。是否有一种向量化的方法,使得给定一个权重对(如果可能的话,N个权值,所有和为1),我对每一行都应用一个协方差矩阵来获得风险(即没有循环)?