From 82fca9583b9d9a5e84c989df7b297454edf4fa70 Mon Sep 17 00:00:00 2001 From: "Zhang, Dabao" Date: Fri, 20 May 2022 14:58:03 -0400 Subject: [PATCH] Update grpsnps.R --- script/cis-eQTL/grpsnps.R | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/script/cis-eQTL/grpsnps.R b/script/cis-eQTL/grpsnps.R index fe179fb..d7ccd54 100644 --- a/script/cis-eQTL/grpsnps.R +++ b/script/cis-eQTL/grpsnps.R @@ -17,10 +17,10 @@ write.table(paste0("$", commonidx), "common.snps.idx", row.names=F, col.names=F, raremap=map[rareidx,] write.table(raremap,"rare.Geno.map",row.names=F,col.names=F,quote=F,sep=" ") -print(paste0("There are ", nrow(raremap), " rare SNPs")) +print(paste0("There are ", nrow(raremap), " ultra-low-frequency SNPs")) lowmap=map[lowidx,] write.table(lowmap,"low.Geno.map",row.names=F,col.names=F,quote=F,sep=" ") -print(paste0("There are ", nrow(lowmap), " low SNPs")) +print(paste0("There are ", nrow(lowmap), " low-frequency SNPs")) commonmap=map[commonidx,] write.table(commonmap,"common.Geno.map",row.names=F,col.names=F,quote=F,sep=" ") print(paste0("There are ", nrow(commonmap), " common SNPs"))