使用3维阵列矩阵计算伽玛

Question

我有一个如下矩阵，样本量为10。

vt=matrix(1:30, ncol=3) 
n=10 
gamma<- array(0, dim=c(3,3,(n-2)))
for (j in 1:(n-2)){
gamma[,,j]<-(1/n)*(t(vt[(j+1):(n-1),])%*%vt[j:(n-2),]
}ter code here

我正在尝试使用如上所述的数组函数计算gamma，但它不起作用。 我想要的是使用上述代码的伽玛1到（n-2）的3乘3矩阵。

例如， gamma[,,1]=(1/n)*(t(vt[2:(n-1),])%*%vt[1:(n-2),] which is a 3 by 3 matrix. gamma[,,2]=(1/n)*(t(vt[3:(n-1),])%*%vt[2:(n-2),] which is a 3 by 3 matrix. gamma[,,3]=(1/n)*(t(vt[4:(n-1),])%*%vt[3:(n-2),] which is a 3 by 3 matrix. 等等....

您是否知道如何修复上述代码？

Answer 1

vt=matrix(1:30, ncol=3) 
n=10 
gamma<- array(0, dim=c(3,3,(n-2)))
for (j in 1:(n-2)){
gamma[,,j]<-(1/n)*t(vt[(j+1):(n-1),])%*%vt[j:(n-2),]
}

我想我明白了。 现在这段代码正在运行。