Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
Next revision Both sides next revision
jlim [2016/11/18 16:16]
schun created
jlim [2016/11/18 16:24]
schun
Line 2: Line 2:
  
 JLIM is a cross-trait test of shared causal effect, which is described in Chun et al. ([[http://​biorxiv.org/​content/​early/​2016/​05/​12/​053165|BioRxiv]]). JLIM tests whether two traits – main and secondary – are driven by shared causal effect or not. Typically, the main trait is a large GWAS study, and secondary trait can be an expression Quantitative Trait Loci (eQTL) association study. For main trait, JLIM takes only summary-level association statistics, but for secondary trait, it requires genotype-level data to generate permutation-based null distribution. JLIM is a cross-trait test of shared causal effect, which is described in Chun et al. ([[http://​biorxiv.org/​content/​early/​2016/​05/​12/​053165|BioRxiv]]). JLIM tests whether two traits – main and secondary – are driven by shared causal effect or not. Typically, the main trait is a large GWAS study, and secondary trait can be an expression Quantitative Trait Loci (eQTL) association study. For main trait, JLIM takes only summary-level association statistics, but for secondary trait, it requires genotype-level data to generate permutation-based null distribution.
 +
 +===== Download =====
 +{{ :​jlim-1.0.2.tar.gz | JLIM release v1.0.2 }}
  
 ===== How to install ===== ===== How to install =====
Line 117: Line 120:
 </​code>​ </​code>​
   * Generate secondary trait genotypes (.ped.gz files).   * Generate secondary trait genotypes (.ped.gz files).
-  * Generate secondary trait association and mperm files. For example, ​``` plink --bfile region1 --pheno gene1tx.pheno --linear --out LCL.gene1 --covar yourcov\\ +  * Generate secondary trait association and mperm files. For example, ​ 
---mperm 1000 --merpm-save-all --chr yourCHR --from-bp yourSTARTBP --end-bp yourENDBP+<​code>​ 
 +plink --bfile region1 --pheno gene1tx.pheno --linear --out LCL.gene1 --covar yourcov --mperm 1000 --merpm-save-all --chr yourCHR --from-bp yourSTARTBP --end-bp yourENDBP
  
 gzip LCL.gene1.assoc.linear gzip LCL.gene1.assoc.linear
  
-gzip LCL.gene1.mperm.dump.all ​``` - Run jlim_gencfg.sh and run_jlim.sh similarly as the above.+gzip LCL.gene1.mperm.dump.all 
 +</​code>​ 
 +  * Run jlim_gencfg.sh and run_jlim.sh similarly as the above.