数据如下
TCGA_id type subtyp LL0YNC03-29C1.1 PLCXD1
1 TCGA-OR-A5J1-01 ACC 01 -9.9658 -9.9658
2 TCGA-OR-A5J2-01 ACC 01 -9.9658 -9.9658
3 TCGA-OR-A5J3-01 ACC 01 -9.9658 -9.9658
4 TCGA-OR-A5J5-01 ACC 01 -9.9658 -9.9658
5 TCGA-OR-A5J6-01 ACC 01 -9.9658 -9.9658
创造了一个这个函数
pancor <- function(gene1,gene2,data){
data1 <- split(data,data$type)
do.call(rbind,lapply(data, function(x){
dd <- cor.test(as.numeric(x[,gene1]),as.numeric(x[,gene2]),type="pearson")
data.frame(type=x$type[1],cor=dd$estimate,p.value=dd$p.value )
}))
}
运行
plotdf <- pancor("LL0YNC03-29C1.1 ","PLCXD1",finally
后
Error in h(simpleError(msg, call)) : error in evaluating the argument 'args' in selecting a method for function 'do.call': incorrect number of dimensions
出现这个报错,我要怎么解决呀?谢谢了