ANGSD: Analysis of next generation Sequencing Data

Latest tar.gz version is (0.938/0.939 on github), see Change_log for changes, and download it here.

HWE and Inbreeding estimates

From angsd
Revision as of 13:32, 25 June 2013 by Thorfinn (talk | contribs)
Jump to navigation Jump to search

Filipe G. Vieira has been extending the ANGSD feature set by

  1. A method for estimating inbreeding (deviation from Hardy-Weinberg) based on output from ANGSD.
  2. Some parts of the ANGSD code will work for inbreeding estimates.

To compile angsd with support for deviations from HWE in the frequency estimates you need to

cd angsd0.539/
make inbreed

Verify that it worked by

./angsd_inbreed -doMaf
ommand:
./angsd_inbreed -doMaf 
	-> angsd version: 0.539	 build(Jun 25 2013 13:27:43)
	-> Analysis helpbox/synopsis information:
------------------------
analysisMaf.inbreed.cpp:
-doMaf	0
	1: BFGS frequency (known major minor)
	2: EM frequency (known major minor)
	4: BFGS frequency (unknown major minor)
	8: EM frequency (unknown major minor)
	16: frequency from genotype probabilities
	32: alleleCounts based method (known major minor)
	-doSNP	0
	-minMaf	0.010000 0
	-minLRT	24.000000 0
	-ref	(null)
	-anc	(null)
	-doZ	0
	-eps	0.001000 [Only used for -doMaf &32]
	-indFname	(null) (file containing individual inbreedcoeficients)
NB these frequency estimators requires major/minor -doMajorMinor

HWE and Inbreeding estimates

For method details see here