一共七个细胞,CIBERSORT absolute score位于idx的七列中(设为i+1)
而IHC cell count就位于前一列(设为i列)
接下来就可以写个for循环做相关性分析了
idx...=1
tmpdata <- data[, c(1, idx[i], idx[i]+1)]
head(tmpdata)
colnames(tmpdata) <- c('Group', 'IHC', 'CIBERSORT...)+min(tmpdata$CIBERSORT))/2,
y = max(tmpdata$log10IHC)*0.95, label = text, size=2) + #...)+min(tmpdata$CIBERSORT))/2,
y = max(tmpdata$log10IHC)*0.95, label = text, size=2) +...xlim(0,max(tmpdata$CIBERSORT)) + ylim(0,max(tmpdata$log10IHC)); plt
outputPdf <- paste(c('..