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.
RealSFS: Difference between revisions
Jump to navigation
Jump to search
No edit summary |
|||
| Line 1: | Line 1: | ||
This program will estimate the SFS based on a . | This program will estimate the SFS based on a .saf file generated from the '''./angsd [options] -doSaf '''. | ||
=Brief overview= | =Brief overview= | ||
<pre> | <pre> | ||
./emOptim2 afile. | ./emOptim2 afile.saf nChr [-start FNAME -P nThreads -tole tole -maxIter -nSites ] | ||
nChr is the number of chromosomes. (twice the number of diploid invididuals) | nChr is the number of chromosomes. (twice the number of diploid invididuals) | ||
</pre> | </pre> | ||
Revision as of 12:54, 17 December 2013
This program will estimate the SFS based on a .saf file generated from the ./angsd [options] -doSaf .
Brief overview
./emOptim2 afile.saf nChr [-start FNAME -P nThreads -tole tole -maxIter -nSites ] nChr is the number of chromosomes. (twice the number of diploid invididuals)
Program defaults to use the EM algorithm for the optimisation, the BFGS method can used by supplying BFGS as the first argument after ./emOptim2
Output
Main results are printed to the stdout. Values are in logspace.
NB
Use as many sites as possible, for more reliable estimates.