2022年2月18日

按分组,每个组画图 

摘要: data %>%group_by(VARIANT_TYPE) %>%do({ p<-ggplot(.,aes(x=TRUTH_SAMPLE,y=NON_REF_GENOTYPE_CONCORDANCE))+geom_point() ggsave(p,filename = paste0(unique( 阅读全文

posted @ 2022-02-18 16:05 BioinformaticsMaster 阅读(32) 评论(0) 推荐(0) 编辑

R 遍历目录下所有子目录中的某类文件

摘要: testvcfdir="/jdfssz1/ST_HEALTH/P18Z10200N0124/qdNB/" testvcfList <- list.files(testvcfdir, pattern = "*.genotype.vcf.gz$", include.dirs = FALSE, full. 阅读全文

posted @ 2022-02-18 14:00 BioinformaticsMaster 阅读(113) 评论(0) 推荐(0) 编辑

导航