全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
557 0
2022-10-12
我在做一个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'至少有两个层次


环境和变量如图所示,请问这个问题应该怎么解决呀?拜谢!


附件列表
1.png

原图尺寸 176.19 KB

1.png

二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群