有没有办法找到'R'中的卡方p值矩阵(属性之间有p值的矩阵)?
例如,考虑iris
数据集。我正在寻找一个矩阵如下:
| | Sepal length | Sepal width | Petal length | Petal width | Species |
|----------------|--------------|-------------|--------------|-------------|---------|
| Sepal length | | | | | |
| Sepal width | | | | | |
| Petal length | | | | | |
| Petal width | | | | | |
| Species | | | | | |
矩阵的元素将是iris
数据集的(i,j)变量的卡方值。
答案 0 :(得分:4)
如果这是你想要的,只考虑其中一列是分类变量,试试这个:
chisqmatrix <- function(x) {
names = colnames(x); num = length(names)
m = matrix(nrow=num,ncol=num,dimnames=list(names,names))
for (i in 1:(num-1)) {
for (j in (i+1):num) {
m[i,j] = chisq.test(x[,i],x[,j],)$p.value
}
}
return (m)
}
mat = chisqmatrix(iris)