在GEO中下載MINiML formatted family file(s)
setwd("C:/AW")
PROBE=c("26540_at","76571_s_at"）#目標(biāo)基因的兩個(gè)探針
PROBE=as.data.frame(PROBE)
colnames(PROBE)="probe"
a=list.files("GSE39582_family") #提取GSE39582_family文件夾下所有文件名
n=length(a) #GSE39582_family文件夾文件數(shù)量
for(i in 1:n){
? b[i]=substring(a[i],1,10) #在文件名中取前10即樣本名
}
dir=paste("./GSE39582_family",a,sep="/")?
merge.data=read.delim(file=dir[1],header=F) #提取第一個(gè)文件
colnames(merge.data)=c("probe",b[1]) #列名
for(i in 2:n){
? new.data=read.delim(file=dir[i],header=F)
? colnames(new.data)=c("probe",b[i])
? merge.data=merge(merge.data,new.data,by="probe") #根據(jù)“probe”列合并文件
}
write.csv(merge.data,file = "./GSE39582_family/merge.csv",row.names=F) #輸出合并文件
expression=merge(merge.data,PROBE,by="probe",all=F) #提取合并文件中目標(biāo)基因列
write.csv(expression,file = "./GSE39582_family/GSE39582C1.csv",row.names=F) #輸出目標(biāo)基因文件
參考：https://blog.csdn.net/weixin_40645816/article/details/82110057
