cellcall分析报错

n_celltype<-n[,  57:64] #1:9,  10:18,  19:27,  28:36,  37:45,  46:54,  55:63,  64:72,  73:81    #1:8  9:16  17:24  25:32  33:40  41:48  49:56  57:64

pathway.hyper.list <- lapply(colnames(n_celltype), function(i){ 
  print(i)
  tmp <- getHyperPathway(data = n_celltype, object = mt, cella_cellb = i, Org="Homo sapiens")   #data:A dataframe of communication score where row name is ligand-receptor and column names is cellA-cellB, stored in the data$expr_l_r_log2_scale slot of S4 object;object:A Cellcall S4 object, the result of function CreateNichConObject() and TransCommuProfile().;cella_cellb:If explore the pathway enriched by paired ligand-receptor dataset between sender cellA and receiver cellB, user can set cella_cellb="A-B".
  return(tmp)
})   #A list of enrichment result
myPub.df <- getForBubble(pathway.hyper.list, cella_cellb=colnames(n_celltype))  #If explore the pathway enriched by paired ligand-receptor dataset between sender cellA and receiver cellB, user can set cella_cellb="A-B".
png("1.png")
plotBubble(myPub.df)
dev.off()
[1] "Erythrocytes-CD14+ Monocytes"
Error in apply(path.list.tmp.df, 1, function(x) { : 
  dim(X) must have a positive length

原因:因为n_celltype的值全部为0?

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值