下载芯片注释文件
图片1.png
图片2.png
library(data.table)
ids=fread('GPL570.annot.gz',data.table = F)[,c(1:4)]
ids$gene=str_split(ids$"Gene symbol",'///',simplify = T)[,1]
row.names(ids) <- ids$ID
ids <- ids[,-1]
合并表达矩阵
将探针注释矩阵ids按差异基因data1行名排序，然后合并
ids<-ids[match(rownames(data1),rownames(ids)),]
data2<-cbind(data1,ids)
筛选差异基因
data3 <- subset(data2,P.Value<0.05)
data4 <- subset(data3,abs(logFC) > 1)
保存结果
write.csv(data2,"GSE23610_10μM_all.csv")
write.csv(data3,"GSE23610_10μM_padj0.05.csv")
write.csv(data4,"GSE23610_10μM_logFC1.csv")
网友评论