我在R中写了一个chi-sq stat test,我的标量chisq正在转换为矩阵,我无法弄清楚原因。
此处使用的变量包括
k是标量
f是kxk矩阵
n是要评估的数据点数
chisq是一个标量,但每次运行此代码时都会成为一个矩阵。
任何帮助将不胜感激
for(i in 1:k){
for(j in 1:k){
chisq<-chisq+(f[i,j]-(n/(k^2)))^2
}
}
答案 0 :(得分:2)
这是一个可重现的代码块:
chisq=0
k=3
f=matrix(runif(k*k),k,k)
n=3
for(i in 1:k){
for(j in 1:k){
chisq<-chisq+(f[i,j]-(n/(k^2)))^2
}
}
并且chisq
作为标量出现。如果我以chisq
作为矩阵开始,我会得到一个矩阵:
chisq = matrix(0,1,1)
所以,我怀疑你正在做后者,不管你说什么。我们和您一样都有可重复的例子。