如何解决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 举报,一经查实,本站将立刻删除。