This is a test version of Biostars. For the public version, visit https://www.biostars.org.
runtime error in LDhat

Hi,

I got a runtime error while running LDhat 'pairwise' program. I need some help to understand the issue so that I can take care of it. The command was:

pairwise -seq chr1.ldhat.sites -loc chr1.ldhat.locs -lk ../new_lk.txt

The error:

Run-time error

allocation failure 2 in imatrix()

On checking the error on github, I found it in this piece of code, but am unable to understand it: https://github.com/auton1/LDhat/blob/master/tools.c

To give a little background, I am working with water buffalo and have 79 samples. In order to generate a look-up table, I downloaded human-centric lookup table 'lk_n192_t0.001' from LDhat github website and used to lkgen program for generating look-up table for my species using the following command:

lkgen lk_n192_t0.001 -nseq 158

I used the generated look-up table on the 'pairwise' program. My final aim is to generate recombination maps for my species that can be used in iHS analysis downstream.

genome snp ldhat population genetics

0 answers

No answers yet.

Log in to answer this question.