数据库代码公开：单肿瘤单基因表达差异

#示例基因为MKI67
geneName = "MKI67"
#示例肿瘤为肺腺癌
CancerName = "LUAD"
#读取表达数据
data = read.table(paste0("Expression.",geneName,".zscore.txt"),sep="\t",header=T,check.names=F)
#读取CancerType是LUAD的行
data = data[data$CancerType == CancerName,]
#保证基因表达量是数值型
data[,"geneName"] = as.numeric(unlist(data[,"geneName"]))
#去除空值
data = na.omit(data)
#wilcoxon计算正常组与肿瘤组是否存在显著性的差异
p = wilcox.test(data[which(data$Type == "Normal"),geneName],
                data[which(data$Type == "Tumor"),geneName])$p.value