数据库代码公开：GEO数据集差异分析
作者：半步博导 时间：2024年6月30日 05:09
阅读量：164
#GEO差异分析同样采用wilcox.test函数，采用scale函数进行zscore标准化
geneName = "MKI67"
tmp = read.csv("GEO_Tumor_Normal.csv",
row.names = NULL,
check.names = F,
header = T,
stringsAsFactors = F)
tmp = tmp[,c(geneName,"Type")]
tmp[,2] = scale(as.numeric(unlist(tmp[,2])),
center = TRUE,
scale = TRUE)
tmp = na.omit(tmp)
p = wilcox.test(tmp[which(tmp$Type == "Normal"),geneName],
tmp[which(tmp$Type == "Tumor"),geneName])$p.value