Cochran检验

提示：文章写完后，目录可以自动生成，如何生成可参考右边的帮助文档

---

前言

非参数统计

例4.14

---

将A、B、C三种榨汁机分给10位家庭主妇使用，用于比较三种榨汁机的受喜爱程度是否相同。

candid1=c(0,0,0,1,0,0,0,0,0,1)
candid2=c(1,1,0,1,0,1,0,0,1,1)
candid3=c(1,1,1,1,1,1,1,1,1,0)
candid=matrix(c(candid1,candid2,candid3),10,3)
ni.=apply(candid,1,sum)
n.j=apply(candid,2,sum)
k=ncol(candid)
Q=(k-1)*((k*sum(n.j^2)-sum(n.j)^2))/(k*sum(ni.)-sum(ni.^2))
pvalue=pchisq(Q,k-1,lower.tail=F)
pvalue
############自编函数######################
candid1=c(0,0,0,1,0,0,0,0,0,1)
candid2=c(1,1,0,1,0,1,0,0,1,1)
candid3=c(1,1,1,1,1,1,1,1,1,0)
candid=matrix(c(candid1,candid2,candid3),10,3)
calculate_pvalue <- function(candid) {
  ni. <- apply(candid, 1, sum)
  n.j <- apply(candid, 2, sum)
  k <- ncol(candid)
  Q <- (k-1)*((k*sum(n.j^2)-sum(n.j)^2))/(k*sum(ni.)-sum(ni.^2))
  pvalue <- pchisq(Q, k-1, lower.tail = FALSE)
  return(pvalue)
}

# 使用自定义函数计算p值
pvalue <- calculate_pvalue(candid)
pvalue