Skip to content

Commit

Permalink
Update grpsnps.R
Browse files Browse the repository at this point in the history
  • Loading branch information
zhangdb authored May 20, 2022
1 parent 3387272 commit 82fca95
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions script/cis-eQTL/grpsnps.R
Original file line number Diff line number Diff line change
Expand Up @@ -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"))

0 comments on commit 82fca95

Please sign in to comment.