Skip to content

Commit

Permalink
user output refine
Browse files Browse the repository at this point in the history
  • Loading branch information
Zhongli Jiang committed May 2, 2022
1 parent 92fc3e2 commit 37d5dbe
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 3 deletions.
2 changes: 1 addition & 1 deletion script/cis-eQTL/common.ciseQTL.r
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
2 changes: 1 addition & 1 deletion script/cis-eQTL/low.ciseQTL.r
Original file line number Diff line number Diff line change
Expand Up @@ -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=" ")
Expand Down
2 changes: 1 addition & 1 deletion script/cis-eQTL/rare.ciseQTL.r
Original file line number Diff line number Diff line change
Expand Up @@ -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=" ")
Expand Down

0 comments on commit 37d5dbe

Please sign in to comment.