※ 引述《qwer098 (^^)》之銘言:
: 若Y_t=(y_t,....,y_(t-p+1))'
: 現在想算Y_t*Y_(t-1)'的和(t=p+1~n)
a=(y_n,y_(n-1),...,y_1) 這是你的實際data
n= n跟p是你要給定的一個值
p=
下面是不變的code
b=a%*%t(a)
c=0
for( i in 1:(n-p))
{
c=c+sum(b[i:(i+p-1),(i+1):(i+p)])
}
然後c就是你要的東西了,如果我沒理解錯誤的話
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.109.40.59