LDhat icon indicating copy to clipboard operation
LDhat copied to clipboard

Composite-likelihood and rho plot

Open Homap opened this issue 4 years ago • 3 comments

Hello,

Thank you for your software. I was wondering what does it mean when the composite likelihood plateaus for different values of rho as shown in the figure:

Figure_1

Thank you!

Homap avatar Nov 03 '20 21:11 Homap

You get such a profile when there is essentially no LD in the region - i.e. the data are consistent with a model in which there is 'infinite' recombination.

mcvean avatar Nov 04 '20 08:11 mcvean

Thank you very much for your reply. I worked a bit more with the data and I still get the same result. The data are of 2000 SNPs (genotypes) in a region of about 1 mega base. Does this mean there is an error with the data or could this have some biological meaning? Wouldn't we expect at least some level of LD in a region of about 1 Mb?

Thanks a lot again!

Homap avatar Nov 05 '20 21:11 Homap

Hi again,

It was indeed a too much large region. I took a smaller region of about 50 SNPs (about 35 Kb) and now I see a peak in my composite likelihood plot.

Thanks again, Homa

LD_black_50SNPs

Homap avatar Nov 10 '20 16:11 Homap