Estimate recombination rates from population genetic data
C PostScript TeX C++ R Makefile
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Permalink
Failed to load latest commit information.
Example Initial Commit. Version 2.2a Jul 23, 2015
lk_files Initial Commit. Version 2.2a Jul 23, 2015
manual_source Update manual.tex Jul 24, 2015
MCMC.h Initial Commit. Version 2.2a Jul 23, 2015
README.md Update README.md Aug 11, 2015
block.h Initial Commit. Version 2.2a Jul 23, 2015
coalescent.r Added coalescent.r Jul 23, 2015
compLK.cpp Initial Commit. Version 2.2a Jul 23, 2015
compLK.h Initial Commit. Version 2.2a Jul 23, 2015
complete.c Initial Commit. Version 2.2a Jul 23, 2015
convert.c Initial Commit. Version 2.2a Jul 23, 2015
data.h Initial Commit. Version 2.2a Jul 23, 2015
fin.c Initial Commit. Version 2.2a Jul 23, 2015
fin.h Initial Commit. Version 2.2a Jul 23, 2015
gamma.c Initial Commit. Version 2.2a Jul 23, 2015
gamma.h Initial Commit. Version 2.2a Jul 23, 2015
hotspot.h Initial Commit. Version 2.2a Jul 23, 2015
ldhat.h Initial Commit. Version 2.2a Jul 23, 2015
likelihood.h Initial Commit. Version 2.2a Jul 23, 2015
lkgen.c Initial Commit. Version 2.2a Jul 23, 2015
makefile Initial Commit. Version 2.2a Jul 23, 2015
manual.pdf Updated manual. Jul 24, 2015
pair_int.c Initial Commit. Version 2.2a Jul 23, 2015
pairdip.c Initial Commit. Version 2.2a Jul 23, 2015
pairlk.c Initial Commit. Version 2.2a Jul 23, 2015
params.h Initial Commit. Version 2.2a Jul 23, 2015
rhomap.cpp Initial Commit. Version 2.2a Jul 23, 2015
rhomap.h Initial Commit. Version 2.2a Jul 23, 2015
rhomap_tools.cpp Initial Commit. Version 2.2a Jul 23, 2015
rhomap_tools.h Initial Commit. Version 2.2a Jul 23, 2015
seqtools.c Initial Commit. Version 2.2a Jul 23, 2015
seqtools.h Initial Commit. Version 2.2a Jul 23, 2015
snp_sim.c Initial Commit. Version 2.2a Jul 23, 2015
snp_sim.h Initial Commit. Version 2.2a Jul 23, 2015
stat.c Initial Commit. Version 2.2a Jul 23, 2015
tools.c Initial Commit. Version 2.2a Jul 23, 2015
tools.h Initial Commit. Version 2.2a Jul 23, 2015

README.md

LDhat

LDhat is a package written in the C and C++ languages for the analysis of recombination rates from population genetic data. The package is available either as C/C++ source code.

The following programs are included:

  • convert: Simple manipulation of DNA sequence data
  • pairwise: Parametric and nonparametric analyses of recombination
  • interval: Estimation of variable recombination rates
  • rhomap: Estimation of variable recombination rates in the presence of hotspots
  • fin: Simulation of genetic variation data in LDhat format

Other programs in the package can be used to generate lookup tables and summarise the output of analyses. Programs are run from the command line. A PDF manual containing information about installing, running and interpreting results from the package should be read before proceeding. An example data set - LPL sequences from a Finnish population (Nickerson et al., 1998) - is included in the package.

Installation

A makefile is included, so it should be possible to compile LDhat by just typing make. The resulting executables are run from the command line.

Likelihood Lookup Tables

A number of precomputed likelihood lookup tables are available for download. Using these tables in conjunction with lkgen will speed up analyses considerably, as the calculation of 2-locus coalescent likelihoods is the most computationally intesive aspect of the core LDhat algorithms. These tables all use the same grid for 4Ner (rho), 101 points evenly spaced between 0 and 100. They differ only in the number of sequences (n) and the value of population mutation rate parameter (theta). Note that lkgen can be used to generated lookup tables from these for smaller numbers of chromosomes, and that minor differences in theta do not appear to strongly influence the results.

All of the likelihood lookup files are gzip compressed, and should be decompressed before use.

Getting help

Send questions, comments, or suggestions to ldhat-help@lists.sourceforge.net