r中的卡方p值矩阵

时间:2015-09-23 06:34:12

标签: r chi-squared

有没有办法找到'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)变量的卡方值。

1 个答案:

答案 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)