R:从for循环中的卡方检验中提取p值

问题描述

我正在使用for循环在R中进行大约10,000个卡方检验。测试本身似乎运行良好,但是我在提取p值时遇到问题。理想情况下,我想将p值放入预先存在的数据框中。 (我也尝试过使用Magic For包,但是(我想)是因为卡方测试结果的输出方式,所以整个输出都在一列中,但这仍然不是我真正想要的。)

这是我到目前为止所拥有的。

Pat1results<-data.frame(Markers=Pat1_markers,p=NA) #dataframe I want the p-values to end up in

test<-for (i in 1:nrow(Pat1_obs)){
  print(i)
  ChiSq<-chisq.test(c(Pat1_obs[i,1],Pat1_obs[i,2]),correct=FALSE,p=AllHetHomoexpected)$p.value
  Pat1results[i,"p"]<-ChiSq
}

但是,应该包含p值的列是所有NA。如何在其中获得p值?

(很抱歉,我很抱歉,我仍然很熟悉R)

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)