R適合しない引数

Question

何千もの質問がありましたが、私はまだ問題を解決できないようです。私は、この二重のループを実行しようとしています

：

rho=matrix(0,12,1) 

for (i in 1:k){ 
    for (j in 1:n-k){ 
    temp=((1/(sigma^2))*(1/(n-i)))*((as.matrix(vR)[j,]-mean)%*%(as.matrix(vR)[i+j,]-mean)) 
    rho[i,]=temp 
    } 
}  

を私はエラーを取得しておいてください。

Error in (as.matrix(vR)[j, ] - mean) %*% (as.matrix(vR)[i + j, ] - mean) : non-conformable arguments

また、私の式は次のとおりです。

(1/sigma^2) * (1/(t-k)) * sum[(vR(t)-mean)*(vR(t+k)-mean)] 

合計されます " t = 1からTkに行く大きなシグマ "

私はやっているように動かすことは可能ですか？

Answer 1

crossprod関数のベクトルの「行列乗算」を処理するためのより安全な方法があります。さらに、私はあなたのエラーは、実際には、配列関数の演算子の優先順位から、大部分の幹かもしれないと思う「：」オーバー内側のループインデックスの設定で「 - 」機能：

for (i in 1:k){ 
for (j in 1:(n-k)){ # note inner parentheses, safer to always use with "(" 
temp=((1/(sigma^2))*(1/(n-i)))*crossprod(vR[j,]-mean,vR[i+j,]-mean) 
rho[i,]=temp 
} 
}