在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