From 37d5dbe08ecefb6cefbc7eb8b3e95f3be7473732 Mon Sep 17 00:00:00 2001 From: Zhongli Jiang Date: Sun, 1 May 2022 20:11:10 -0400 Subject: [PATCH] user output refine --- script/cis-eQTL/common.ciseQTL.r | 2 +- script/cis-eQTL/low.ciseQTL.r | 2 +- script/cis-eQTL/rare.ciseQTL.r | 2 +- 3 files changed, 3 insertions(+), 3 deletions(-) diff --git a/script/cis-eQTL/common.ciseQTL.r b/script/cis-eQTL/common.ciseQTL.r index 88e9ab5..3dd3b8d 100755 --- a/script/cis-eQTL/common.ciseQTL.r +++ b/script/cis-eQTL/common.ciseQTL.r @@ -19,7 +19,7 @@ for (i in 1:len){ fit=lm(y[,idx[i,1]]~x[,idx[i,2]]) p[i]=summary(fit)$coefficients[2,4] } - if( (i%%100)==0) print(i) + if( (i%%10000)==0) print(i) } p=cbind(idx,p) # Col 1 is index of gene, Col 2 is index of its cis-SNP, Col 3 is p-value diff --git a/script/cis-eQTL/low.ciseQTL.r b/script/cis-eQTL/low.ciseQTL.r index 5062644..2e9bab9 100755 --- a/script/cis-eQTL/low.ciseQTL.r +++ b/script/cis-eQTL/low.ciseQTL.r @@ -93,7 +93,7 @@ for (i in YYstartYY:YYendYY){ } } - if((i%%100)==0) print(i) + if((i%%10000)==0) print(i) } write.table(w,"low.ciseQTL.weightYYY",row.names=F,col.names=F,quote=F,sep=" ") diff --git a/script/cis-eQTL/rare.ciseQTL.r b/script/cis-eQTL/rare.ciseQTL.r index bf7b2bf..7e5997d 100644 --- a/script/cis-eQTL/rare.ciseQTL.r +++ b/script/cis-eQTL/rare.ciseQTL.r @@ -94,7 +94,7 @@ for (i in YYstartYY:YYendYY){ } } - if((i%%100)==0) print(i) + if((i%%10000)==0) print(i) } write.table(w,"rare.ciseQTL.weightYYY",row.names=F,col.names=F,quote=F,sep=" ")