学习果子老师的操作,目的是掌握循环和dplyr的应用
(ftp://ftp.ncbi.nlm.nih.gov/geo/series/GSE48nnn/GSE48213/suppl/GSE48213_RAW.tar))首先下载源文件,共56个文件,下载完成后,解压到GSE48213_RAW。
nameList <- list.files("GSE48213_RAW/")
此代码是看在GSE48213_RAW/下有多少个文件。
新创建一个矩阵,将文件逐个放到矩阵里
matrix <- read.table(paste0("GSE48213_RAW/",nameList[1]),header = T)
先试第一个文件。
for (i in 2:length(nameList)){
matrix <- inner_join(matrix,
read.table(paste0("GSE48213_RAW/",nameList[i]),header = T),
by="EnsEMBL_Gene_ID")
}
操作就完成了
image.png
网友评论