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=" ")