R小tip（十一）批量读取count文件

这次我们来介绍下用R批量RNA-seq的count文件

directory <- "/home/data/"   #count文件路径
Files <- grep(".htseqcount$",list.files(directory),value=TRUE) #正则表达式匹配文件 
celltype <- c(rep("A",3),rep("B",3)) #data.frame的colnames
multitable <- data.frame(Name = unlist(lapply(strsplit(Files,"\\."),"[[",1)),
                         fileName = Files,
                         celltype = celltype)
#读入一个data.frame，此时若干htseqcount表就合并成大的data.frame