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.

SFS Estimation: Difference between revisions

From angsd
Jump to navigation Jump to search
No edit summary
Line 4: Line 4:
This is a 2 step procedure first generate a ".sfs" file using -realSFS 1. And then do an optimization of the .sfs file which will estimate the Site frequency spectrum.
This is a 2 step procedure first generate a ".sfs" file using -realSFS 1. And then do an optimization of the .sfs file which will estimate the Site frequency spectrum.


<classdiagram>
<classdiagram type="dir:LR">
  [sequence data]->GL[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]
  [sequence data]->GL[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]
[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]->realSFS[.sfs file]
[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]->realSFS[.sfs file]

Revision as of 15:08, 10 October 2012

realSFS

This method will estimate the site frequency spectrum, the method is described in Nielsen2012.

This is a 2 step procedure first generate a ".sfs" file using -realSFS 1. And then do an optimization of the .sfs file which will estimate the Site frequency spectrum.

<classdiagram type="dir:LR">

[sequence data]->GL[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]

[genotype likelihoods|SAMtools;GATK;SOAPsnp;Kim et.al]->realSFS[.sfs file] [.sfs file]->optimize[.ml file]

</classdiagram>


-realSFS 1
an sfs file will be generated.
-realSFS 2
snpcalling (not implemented, in this angsd)
-realSFS 4
genotypecalling (not implemented, int this angsd)