FastNgsAdmixOld: Difference between revisions

From software
Jump to navigation Jump to search
No edit summary
No edit summary
Line 19: Line 19:
<pre>
<pre>
YRI, CEU, Sardinian, WestGreenland, CHB, JPT, Sindhi, Papuan, Chukchi, Karitiana
YRI, CEU, Sardinian, WestGreenland, CHB, JPT, Sindhi, Papuan, Chukchi, Karitiana
<\pre>
</pre>


Preset reference populations can be altered.
Preset reference populations can be altered.

Revision as of 21:13, 21 September 2014

This page contains information about the program called FastNGSadmix, which is a very fast tool for finding admixture proportions from NGS data of a single individual. It is based on genotype likelihoods. The program is written in R.


Installation

wget popgen.dk/??
OR
Use shiny: popgen.dk/???

Run example

First download some example test files which has been generated on basis of data from the 1000 genomes project.

wget popgen.dk/??????/popfreq.RData
wget popgen.dk/??????/beagle.gz

tralalala

YRI, CEU, Sardinian, WestGreenland, CHB, JPT, Sindhi, Papuan, Chukchi, Karitiana

Preset reference populations can be altered.

HOME$ Rscript NGS_Admix_rscript.r beagle_GL='indi_genotypelikelihood.beagle'  popfreq_RDATA='pop_freq_all.RData' outputfile='test'

Input Files

Input files are contains genotype likelihoods in genotype likelihood beagle input file format [1]. We recommend [ANGSD] for easy transformation of Next-generation sequencing data to beagle format.

Example of a beagle genotype likelihood input file for 3 individuals.

marker       allele1  allele2   Ind0      Ind0    Ind0
1_14000023      1       0       0.941    0.058    0.000
1_14000072      2       3       0.709    0.177    0.112
1_14000113      0       2       0.855    0.106    0.037
1_14000202      2       0       0.835    0.104    0.060
...