name <- rownames(rt_expr) #基因名
immune <- function(genename){
y <- as.numeric(rt_expr) #单一基因的表达量表格
do.call(rbind,lapply(colnames(rt),function(x){ #rt列明与rt_expr的行名相同
test <- cor.test(rt[,x],y,type="spearman")
cor <- test$estimate
p.value <- test$p.value
data.frame(gene=genename,symbol=x,correlation=cor,pvalue=p.value)
}))
}
data <- immune(name)
函数的参数的传入也可以多个,以及固定参数如下:
compute <- function(x,y,d=2){
....
}
链接表
-
近期文章
近期评论
- daizao发表在《TCGA与dbGaP账号申请(dbGaP申请指南)》
- zzz发表在《TCGA与dbGaP账号申请(dbGaP申请指南)》
- daizao发表在《TCGA与dbGaP账号申请(dbGaP申请指南)》
- liyuqing发表在《TCGA与dbGaP账号申请(dbGaP申请指南)》
- daizao发表在《TCGA与dbGaP账号申请(dbGaP申请指南)》
分类目录
分类目录
感谢博主分享