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]
批量读取文件
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 参考：https://blog.csdn.net/leo_888/article/details/80681184...
