This is a beta test.
Question: runtime error in LDhat
0
Entering edit mode

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.

Login before adding your answer.

Powered by the version 1.6