Skip to content

Commit

Permalink
Bump release to version 1.0.4
Browse files Browse the repository at this point in the history
  • Loading branch information
mkerin committed Sep 23, 2022
1 parent 44120cb commit 5acc9f8
Show file tree
Hide file tree
Showing 4 changed files with 14 additions and 14 deletions.
2 changes: 1 addition & 1 deletion CMakeLists.txt
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
cmake_minimum_required(VERSION 3.16.9) # "CMake should be at least as new as your compiler"
project(LEMMA VERSION 1.0.3)
project(LEMMA VERSION 1.0.4)
include(CTest)

cmake_policy(SET CMP0110 OLD)
Expand Down
6 changes: 3 additions & 3 deletions docs/gplemma.md
Original file line number Diff line number Diff line change
Expand Up @@ -18,7 +18,7 @@ The phenotype has been simulated to have:
## Getting started
To run the GPLEMMA method on the dataset given above, run the following commands
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--gplemma --random-seed 1 \
--pheno example/pheno.txt.gz \
--streamBgen example/n5k_p20k_example.bgen \
Expand All @@ -29,7 +29,7 @@ This should return heritability estimates of h2-G = 0.229 (0.028) and h2-GxE = 0

To run the MEMMA method on the same dataset, run the following commands
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--RHEreg --random-seed 1 \
--pheno example/pheno.txt.gz \
--streamBgen example/n5k_p20k_example.bgen \
Expand All @@ -54,6 +54,6 @@ rather than being restricted to a single node.

To set the number of cores on the commandline explicitly, use
```
mpirun -n <cores> build/lemma_1_0_3
mpirun -n <cores> build/lemma_1_0_4
```

2 changes: 1 addition & 1 deletion docs/install.md
Original file line number Diff line number Diff line change
Expand Up @@ -34,7 +34,7 @@ ctest --test-dir build

View the commandline options supported by LEMMA
```
./build/lemma_1_0_3 -h
./build/lemma_1_0_4 -h
```


Expand Down
18 changes: 9 additions & 9 deletions docs/lemma.md
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@ for cc in `seq 1 22`; do
echo "example/n5k_p20k_example_chr${cc}.bgen" >> example/bgen_filenames.txt;
done
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--pheno example/pheno.txt.gz \
--environment example/env.txt.gz \
--VB \
Expand Down Expand Up @@ -63,7 +63,7 @@ The LEMMA algorithm is modular, and so each step can be performed separately as

### Running the LEMMA variational inference algorithm
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--VB \
--pheno example/pheno.txt.gz \
--environment example/env.txt.gz \
Expand All @@ -74,7 +74,7 @@ In this case the algorithm should converge in 59 iterations.

### Association testing with imputed SNPs
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--singleSnpStats --maf 0.01 \
--pheno example/pheno.txt.gz \
--resid-pheno example/inference_converged_yhat.out \
Expand All @@ -88,7 +88,7 @@ For analyses of large genomic datasets it may be useful to parallelize associati

### Heritability estimation
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--RHEreg --random-seed 1 \
--pheno example/pheno.txt.gz \
--mStreamBgen example/bgen_filenames.txt \
Expand All @@ -113,7 +113,7 @@ rather than being restricted to a single node.

To set the number of cores on the commandline explicitly, use
```
mpirun -n <cores> build/lemma_1_0_3
mpirun -n <cores> build/lemma_1_0_4
```

### Precomputing the dXtEEX array
Expand All @@ -139,7 +139,7 @@ zcat example/dxteex_chr*.out.gz > example/dxteex.out.gz
```
Then provide the file `example/dxteex.out.gz` to LEMMA with the commandline flag `--dxteex`.
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--VB \
--pheno example/pheno.txt.gz \
--environment example/env.txt.gz \
Expand All @@ -158,7 +158,7 @@ To convert into the file format expected by LEMMA we have provided a brief Rscri

Then run the heritability analysis as follows
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--RHEreg --n-RHEreg-samples 20 --n-RHEreg-jacknife 100 --random-seed 1 \
--pheno example/pheno.txt.gz \
--bgen example/n5k_p20k_example.bgen \
Expand All @@ -170,15 +170,15 @@ mpirun build/lemma_1_0_3 \
### Resuming from a previous parameter state
In case of runtime crashes, LEMMA can save the parameter state at periodic intervals by providing the commandline flag `--resume-from-state`. LEMMA can then subsequently resume inference from this saved state. For example
```
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--VB \
--pheno example/pheno.txt.gz \
--environment example/env.txt.gz \
--bgen example/n5k_p20k_example.bgen \
--state-dump-interval 10 \
--out example/inference.out
mpirun build/lemma_1_0_3 \
mpirun build/lemma_1_0_4 \
--VB \
--pheno example/pheno.txt.gz \
--environment example/env.txt.gz \
Expand Down

0 comments on commit 5acc9f8

Please sign in to comment.