PCAngsd: Difference between revisions
No edit summary |
No edit summary |
||
Line 24: | Line 24: | ||
python pcangsd.py -h | python pcangsd.py -h | ||
# | # Only estimate covariance matrix | ||
python pcangsd.py -beagle test.beagle.gz -o test | python pcangsd.py -beagle test.beagle.gz -o test | ||
# Estimate inbreeding coefficients | # Estimate inbreeding coefficients (and covariance matrix) | ||
python pcangsd.py -beagle test.beagle.gz -inbreed 1 -o test | python pcangsd.py -beagle test.beagle.gz -inbreed 1 -o test | ||
# Perform selection scan | # Perform selection scan (and covariance matrix) | ||
python pcangsd.py -beagle test.beagle.gz -selection 1 -o test | python pcangsd.py -beagle test.beagle.gz -selection 1 -o test | ||
</pre> | </pre> | ||
=Input= | =Input= | ||
The only input PCAngsd needs and accepts are genotype likelihoods in [http://faculty.washington.edu/browning/beagle/beagle.html | The only input PCAngsd needs and accepts are genotype likelihoods in [http://faculty.washington.edu/browning/beagle/beagle.html Beagle] format. [http://popgen.dk/angsd ANGSD] can be easily be used to compute genotype likelihoods and output them in the required Beagle format. | ||
<pre> | <pre> | ||
Line 45: | Line 45: | ||
=Using PCAngsd= | =Using PCAngsd= | ||
All the different options in PCAngsd are listed here. | All the different options in PCAngsd are listed here. **All required analyses must be run in the same command!** (*For now...*) | ||
PCAngsd will always compute the covariance matrix. It uses the computed principal components to estimate the individual allele frequencies in an iterative procedure. This procedure is performed until the individual allele frequencies have converged. | PCAngsd will always compute the covariance matrix. It uses the computed principal components to estimate the individual allele frequencies in an iterative procedure. This procedure is performed until the individual allele frequencies have converged. | ||
; -beagle [Beagle filename] | ; -beagle [Beagle filename] | ||
Path to file of the genotype likelihoods in | Path to file of the genotype likelihoods in Beagle format. | ||
; -beaglelist [filelist] | ; -beaglelist [filelist] | ||
Parse a file with a list of multiple | Parse a file with a list of multiple Beagle files, e.g. if the genotype likelihoods have been computed separately for each chromosome. | ||
; -M [int] | ; -M [int] | ||
Maximum number of iterations for covariance estimation. Only needed in rare cases. (Default: 100) | Maximum number of iterations for covariance estimation. Only needed in rare cases. (Default: 100) | ||
Line 118: | Line 118: | ||
==Example== | ==Example== | ||
<pre> | <pre> | ||
# Estimate covariance matrix, inbreeding coefficients and perform genome selection scan. | # Estimate covariance matrix, inbreeding coefficients, kinship matrix and perform genome selection scan including various filters. | ||
python pcangsd.py -beagle test.beagle.gz -inbreed 2 - | python pcangsd.py -beagle test.beagle.gz -inbreed 2 -kinship -selection 1 -o test -M_tole 1e-3 -inbreed_tole 1e-3 | ||
</pre> | </pre> | ||
==Citation== | ==Citation== |
Revision as of 12:02, 10 August 2017
This page contains information about the program PCAngsd, which estimates the covariance matrix for low depth NGS data in an iterative procedure based on genotype likelihoods. Based on the population structure inference PCAngsd is able to estimate individual allele frequencies. By incorporating these allele frequencies in Empirical Bayes approaches, PCAngsd can perform PCA (estimate covariance matrix), call genotypes, estimate inbreeding coefficients (per-individual and per-site) and perform a genome selection scan using principal components in structured populations. The entire program is written in Python 2.7.
Download
The program can be downloaded from Github: https://github.com/Rosemeis/pcangsd
git clone https://github.com/Rosemeis/pcangsd.git; cd pcangsd/
The following Python packages are needed to run PCAngsd (found in all popular distributions): numpy and pandas.
PCAngsd should work on all platforms meeting the requirements but server-side usage is recommended.
Quick start
# See all options in PCAngsd python pcangsd.py -h # Only estimate covariance matrix python pcangsd.py -beagle test.beagle.gz -o test # Estimate inbreeding coefficients (and covariance matrix) python pcangsd.py -beagle test.beagle.gz -inbreed 1 -o test # Perform selection scan (and covariance matrix) python pcangsd.py -beagle test.beagle.gz -selection 1 -o test
Input
The only input PCAngsd needs and accepts are genotype likelihoods in Beagle format. ANGSD can be easily be used to compute genotype likelihoods and output them in the required Beagle format.
./angsd -GL 1 -out genoLikes -nThreads 10 -doGlf 2 -doMajorMinor 1 -doMaf 2 -SNP_pval 1e-6 -bam bam.filelist
See ANGSD for more info on how to compute the genotype likelihoods and call SNPs.
Using PCAngsd
All the different options in PCAngsd are listed here. **All required analyses must be run in the same command!** (*For now...*)
PCAngsd will always compute the covariance matrix. It uses the computed principal components to estimate the individual allele frequencies in an iterative procedure. This procedure is performed until the individual allele frequencies have converged.
- -beagle [Beagle filename]
Path to file of the genotype likelihoods in Beagle format.
- -beaglelist [filelist]
Parse a file with a list of multiple Beagle files, e.g. if the genotype likelihoods have been computed separately for each chromosome.
- -M [int]
Maximum number of iterations for covariance estimation. Only needed in rare cases. (Default: 100)
- -M_tole [float]
Tolerance value for the iterative covariance matrix estimation. (Default: 1e-4)
- -EM [int]
Maximum number of EM iterations for computing the population allele frequencies. (Default: 200)
- -EM_tole [float]
Tolerance value in EM algorithm for population allele frequencies estimation. (Default: 1e-4)
- -e [int]
Manually select the number of eigenvalues to use in modelling of individual allele frequencies. (Default: Automatically tested)
- -reg
(Not fully tested!) Toogle to use Tikhonov regularization in modelling of individual allele frequencies to penalize lesser important PCs. May also help on convergence.
- -o [prefix]
Set the prefix for all output files created by PCAngsd (Default: "pcangsd").
Call genotypes
Genotypes can be called from posterior genotype probabilities incorporating the individual allele frequencies in prior.
- -geno [float]
Call genotypes with defined threshold.
- -genoInbreed [float]
Call genotypes with defined threshold also taken inbreeding into account. -inbreed [int] is required.
Inbreeding
Per-individual inbreeding coefficients incorporating population structure can be computed using three different methods:
- -inbreed 1
A maximum likelihood estimator computed by an EM algorithm. Only allows for F-values between 0 and 1. Based on [1].
- -inbreed 2
Simple estimator also computed by an EM algorithm. Based on [2].
- -inbreed 3
(Not recommended for low depth NGS data!) Estimator using the kinship matrix. Based on PC-Relate.
- -inbreed_iter [int]
Maximum number of iterations for the EM algorithm methods. (Default: 200)
- -inbreed_tole [float]
Tolerance value for the EM algorithms for inbreeding coefficients estimation. (Default: 1e-4)
Per-site inbreeding coefficients incorporating population structure alongside likehood ratio tests for HWE can be computed as follows:
- -inbreedSites
Selection
A genome selection scan can be computed using two different methods:
- -selection 1
Using an extended model of FastPCA. Produces a genome selection scan for all significant PCs.
- -selection 2
(Not fully tested!) Using an extended model of PCAdapt.
LD can also be taken into account when performing selection scans. LD regression has been implemented in PCAngsd but the functionality is not fully tested.
- -LD [int]
Select the window (in bases) of preceding sites to use in regression.
Relatedness
Work in progress.
Estimate kinship matrix:
- -kinship
Automatically estimated if -inbreed 3 has been selected.
Example
# Estimate covariance matrix, inbreeding coefficients, kinship matrix and perform genome selection scan including various filters. python pcangsd.py -beagle test.beagle.gz -inbreed 2 -kinship -selection 1 -o test -M_tole 1e-3 -inbreed_tole 1e-3