用于chisq.test的循环

Question

我有一个包含多个分类列的数据集。我想使用chisq.test为每对列计算p值。正在尝试这个，但它给出了一个错误。

e.g.
dataset
v1  v2  v3  v4
a   q   e   f
a   w   g   h
b   z   e   i
c   q   e   h

a = dataset
value = matrix(nrow=ncol(a), ncol=ncol(a))
for (i = 1:ncol(a)){
  for (j = 1:ncol(a)){
    tbl = table(a[,i],a[,j])
    tst = chisq.test(tbl)
    value[i,j] == tst$p.value
  }
}

任何潜在客户都将获得奖励

Answer 1

您需要对代码进行一些更改，以便将p值存储在value。

中

for (i in 1:ncol(a)) {
  for (j in 1:ncol(a)) {
    tbl <- table(a[, i], a[, j])
    tst <- chisq.test(tbl)
    value[i, j] <- tst$p.value
  }
}