R语言 for循环问题，用于WGCNA,求助大神！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在做一个WGCNA分析，里面有一段代码：后续需要建立一个矩阵，取TCGA_AML_Members和consMembers做fisher检验：

# Numbers of femaleand consensus modules

nTCGA_AML_Mods =length(TCGA_AML_Modules)

nConsMods =length(consModules)

# Initialize tablesof p-values and of the corresponding counts

pTable = matrix(0,nrow = nTCGA_AML_Mods, ncol = nConsMods);

CountTbl =matrix(0, nrow = nTCGA_AML_Mods, ncol = nConsMods);

# Execute all pairwaisecomparisons

for (Tmod in 1: nTCGA_AML_Mods)

for (Cmod in 1: nConsMods)

{

TCGA_AML_Members =(TCGA_AML_Colors == TCGA_AML_Modules[Tmod]);

consMembers =(moduleColors == consModules[Cmod]);

pTable[Tmod, Cmod]= -log10(fisher.test(TCGA_AML_Members, consMembers, alternative = "greater")$p.value);

CountTbl[Tmod, Cmod]= sum(TCGA_AML_Colors == TCGA_AML_Modules[Tmod] & moduleColors ==consModules[Cmod])

}

运行后显示Error in fisher.test(TCGA_AML_Members, consMembers, alternative = "greater") :

  'x'和'y'至少有两个层次

环境和变量如图所示，请问这个问题应该怎么解决呀？拜谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：for循环 求助大神 For CNA R语言 R语言 for循环