问题描述
我正在使用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 (将#修改为@)