差异基因分析筛选出差异基因表达矩阵

setwd("C:/Users/linqingquan/Desktop/GSE_121787");
Sys.setlocale('LC_ALL','Chinese');

GPL_table = read.table('GPL21185-21174.txt',sep="\t",comment.char="#",
                       stringsAsFactors=F,header=T,fill=TRUE,quote="");

GSE_matrix <- read.table('GSE121787_series_matrix.txt',sep="\t",comment.char="!",
                         stringsAsFactors=F,header=T,fill=TRUE);
ID_Sybmol = GPL_table[,c(1,6)];
colnames(ID_Sybmol)[2]="Symbol";

Exp=merge(ID_Sybmol,GSE_matrix,by.x="ID",by.y="ID_REF",all=T);
Exp=Exp[,-1];
View(Exp)

Exp=Exp[Exp$Symbol != "",];
Exp=na.omit(Exp);

Exp1=data.frame(Exp[-grep("/",Exp$"Symbol"),]);
meanfun <- function(x) {
  x1 <- data.frame(unique(x[,1]));
  colnames(x1) <- c("Symbol");
  for (i in 2:ncol(x)) {
    x2