# 稍稍修改一下上面的计算程序即可。
pv=numeric();lpv=numeric() # 初始化结果向量
for (i in 1:10) {
tbdata=matrix(as.integer(dt[i,]),nrow=2,byrow=T)
pv[i]=fisher.test(tbdata)$p.value;lpv[i]=-log(pv[i])
}
cbind(pv,lpv)
输出
pv lpv
[1,] 0.23079919 1.4662073
[2,] 0.69675397 0.3613229
[3,] 0.62425086 0.4712030
[4,] 0.02086171 3.8698401
[5,] 0.02187766 3.8222892
[6,] 1.00000000 0.0000000
[7,] 0.02838429 3.5619194
[8,] 0.03646188 3.3114879
[9,] 0.77493610 0.2549747
[10,] 1.00000000 0.0000000