采用的数据是GEO上的200079973
RMA法预处理normal样本
setwd(“”)
library(affyPLM)
library(affy)
Data<-ReadAffy()
sampleNames(Data)
N=length(Data)
#用RMA预处理数据
<-rma(Data)
#获取表达数据并输出到表格
normal_exprs<-exprs()
probeid<-rownames(normal_exprs)
normal_exprs<-cbind(probeid,normal_exprs)
(normal_exprs,file=“”,sep=‘\t’,quote=F,=F)
RMA法预处理tumor样本
setwd(“”)
library(affyPLM)
library(affy)
Data<-ReadAffy()
sampleNames(Data)
N=length(Data)
#用RMA预处理数据
<-rma(Data)
#获取表达数据并输出到表格
normal_exprs<-exprs()
probeid<-rownames(normal_exprs)
normal_exprs<-cbind(probeid,normal_exprs)
(normal_exprs,file=“”,sep=‘\t’,quote=F,=F)
合并N和T的数据
#setwd(" “)
normal_exprs<-(“”,header=T,sep=”\t")
tumor_exprs<-(“”,header=T,sep=“\t”)
#讲T和N合并
probe_exprs<-merge(normal_exprs,tumor_exprs,by=“probeid”)
(probe_exprs,file=“”,sep=‘\t’,quote=F,=F)