Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
NLykoskoufis committed Jan 26, 2021
1 parent cb7b486 commit cf500a3
Showing 1 changed file with 2 additions and 3 deletions.
5 changes: 2 additions & 3 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -21,10 +21,9 @@ sort -V -k1 -k2,2n -k3,3n H3K4me3_data.bed > H3K4me3_data_sorted.bed

### Step2: Create null file

Before creating the null file you also need to provide a list of variants that are nominally significant eQTLs. These can be provide by using QTLtoos cis --nominal 1. Then you need to only keep the significant snps, like so

Before creating the null file you also need to provide a list of variants that are nominally significant eQTLs. These can be provided by using QTLtools cis --nominal 0.05. Then you need to only keep the var_id. It should be the 8th column
```bash
zcat nominal1_chrALL.txt.gz | awk '{if($12 <0.05) {print $8}}' > nominal_only_significant_snps.txt
zcat nominal005_chrALL.txt.gz | cut -d" " -f8 | sort | uniq > nominal_only_significant_snps.txt
```

```bash
Expand Down

0 comments on commit cf500a3

Please sign in to comment.