批量读取文件

setwd("G:/wz/hg38/ovary_aging/SRP075354_counts_GC/dailt")
rm(list=ls())
gc()
samples=list.files("G:/wz/hg38/ovary_aging/SRP075354_counts_GC/dailt")
samples
files <- file.path(samples)
names(files) <- samples
rm(list="samples")

library(limma)
library(sva)

#获取交集基因
fs=list.files('G:/wz/hg38/ovary_aging/SRP075354_counts_GC/dailt')
a=do.call(cbind,lapply(list.files('G:/wz/hg38/ovary_aging/SRP075354_counts_GC/dailt'), function(x){
  read.table(file.path('G:/wz/hg38/ovary_aging/SRP075354_counts_GC/dailt',x),
             header = T,sep = ',',row.names = 1)
}))
a[1:4,1:4]


# 下面是合并后的表达矩阵添加行名和列名
rownames(a)=rownames(read.table('GSM2157535_Younger_1_CPM.csv.gz',
                                header = T,sep = ',',row.names = 1))
d1=seq(2,40,2)
d1=c(1,d1)
a=a[,d1]