Numpy error during MavericK normalization
*Created by: cathynewman*
I am running MavericK via Structure_threader on a Linux HPC server (RHEL 6), with GCC 6.4.0 and python 3.5.2 (Anaconda), which includes Numpy v. 1.11.3. I did a test run on a data set with 11 individuals and 1,000 SNPs, running K = 1-3, with 3 reps each K. I think all of the runs finished successfully and the outputEvidence.csv file was generated. But I get the following error during what looks like the normalization step:
[EDITED to clarify: the error is during the normalization step of _merging_, after all of the runs have completed and the _merged_ outputEvidence.csv file has been created.]
```
INFO: [32mAll 3 jobs finished successfully.[0m
Traceback (most recent call last):
File "/home/cenewman/.local/bin/structure_threader", line 11, in <module>
sys.exit(main())
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/structure_threader.py", line 332, in main
full_run(arg)
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/structure_threader.py", line 286, in full_run
arg.notests)
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/wrappers/maverick_wrapper.py", line 312, in maverick_merger
bestk = _write_normalized_output(evidence, k_list, ti_in_use)
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/wrappers/maverick_wrapper.py", line 258, in _write_normalized_output
k_list))
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/wrappers/maverick_wrapper.py", line 335, in maverick_normalization
for _ in range(draws)])
File "/home/cenewman/.local/lib/python3.5/site-packages/structure_threader/wrappers/maverick_wrapper.py", line 335, in <listcomp>
for _ in range(draws)])
File "mtrand.pyx", line 1902, in mtrand.RandomState.normal (numpy/random/mtrand/mtrand.c:19771)
ValueError: scale <= 0
```
Here is how I call Structure_threader:
```
export HOME_DIR=/home/cenewman
export WORK_DIR=/work/cenewman/popgen/maverick
cd $WORK_DIR
module load gcc/6.4.0
module load python/3.5.2-anaconda-tensorflow
structure_threader run -K 3 -i data.txt -o output --params parameters.txt -t 20 -mv $HOME_DIR/MavericK-1.0.5/MavericK
```
issue