matrix ST=J(3,3,0)
forvalues i=1(10)40 {
matrix C=A[`i'..(`i'+10-1),1...]
svmat C
matrix accum cov = C1-C3, deviations noconstant
matrix cov=cov/(r(N)-1)
matrix ST=ST+cov
drop C1-C3
}
然后本人用excel验证时,把每个时间的协差阵分别算出来,然后再算这些协差阵的平均数,发现两种方法得出的结果不一致。
因为原文中的stata代码看不懂,也不清楚到底是哪里的问题,有没有朋友能给解释一下?谢谢
附件中是我模拟的数据和计算过程



雷达卡




京公网安备 11010802022788号







