##-Finishedin6.300selapsed(6.030s cpu)# 计算TMB,就这1行代码,简单快速解决问题 coad.tmb<-tmb(maf,captureSize=38,logScale=T)##[1]-1.0-0.50.00.51.01.52.02.5 plot of chunk unnamed-chunk-2 上面的1行代码中的captureSize默认是50,有用30的,38的,35的,五花八门,我就用38了。logScale选择是否进行...
pValue=corT$p.value p1=ggplot(df1,aes(x,y))+ylab(j)+xlab(gene)+geom_point()+geom_smooth(method="lm",formula=y~x)+theme_bw()+stat_cor(method='spearman',aes(x=x,y=y))if(pValue<pFilter){pdf(file=paste0(j,".pdf"),width=5,height=5)print(p1)dev.off()outTab=rbind(outTab,...
# 计算TMB,就这1行代码,简单快速解决问题 coad.tmb <- tmb(maf, captureSize =38, logScale =T) ## [1] -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 plot of chunk unnamed-chunk-2 上面的1行代码中的captureSize默认是50,有用30的,38的,35的,五花八门,我就用38了。logScale选择是否进行log10转换。 结...
(method="lm",formula=y~x)+theme_bw()+stat_cor(method='spearman',aes(x=x,y=y))if(pValue<pFilter){pdf(file=paste0(j,".pdf"),width=5,height=5)print(p1)dev.off()outTab=rbind(outTab,cbind(Cell=j,pValue))}}}#相关性检验进行相关性检验,计算基因表达值(x)与各种免疫细胞类型比例(...