NgsAdmix

From software
Jump to navigation Jump to search

This will contain the program called NGSadmix, which is a very nice tool for finding admixture. It is based on genotype likelihoods. It is a fancy multithreaded c/c++ program.


Latest version is 32 from June 25 2013. It can be found [1]. Older versions can be found here: [2]

Installation

wget popgen.dk/software/NGSadmix/ngsadmix32.cpp 
g++ ngsadmix32.cpp -O3 -lpthread -lz -o NGSadmix

Run example

Assuming we have an input file called input.gz and we assume 3 ancestral populations (-K 3), and that we want to use 4 computing cores (-P 4). The prefix of the output files is myoutfiles (-o)

./ngsAdmix -likes input.gz -K 3 -P 4 -tol 0.00000001 -tolLike50 0.01 -o myoutfiles -minMaf 0.05 

Input Files

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

Options

./NGSadmix
Arguments:
	-likes Beagle likelihood filename
	-K Number of ancestral populations
Optional:
	-fname Ancestral population frequencies
	-qname Admixture proportions
	-outfiles Prefix for output files
	-printInfo print ID and mean maf for the SNPs that were analysed
Setup:
	-seed Seed for initial guess in EM
	-P Number of threads
	-method If 0 no acceleration of EM algorithm
	-misTol Tolerance for considering site as missing
Stop chriteria:
	-tolLike50 Loglikelihood difference in 50 iterations
	-tol Tolerance for convergence
	-dymBound Use dymamic boundaries (1: yes (default) 0: no)
	-maxiter Maximum number of EM iterations
Filtering
	-minMaf Minimum minor allele frequency
	-minLrt Minimum likelihood ratio value for maf>0
	-minInd Minumum number of informative individuals

Output Files

Program outputs 3 files.

  1. PREFIX.log
  2. PREFIX.fopt.gz
  3. PREFIX.qopt
  • The log file contains log information of the run. Commandline used for running the program, what the likelihood is every 50 iterations, and finally how long it took to do the run.
  • The fopt.gz file is an compressed file, which contains an estimate of the frequency for each site for all populations.
  • The qopt file contains the admixture proportions for all individuals.

Examples of the output files are found below.


Log file

Contents of the file log file

	-> Dumping file: tskSim/tsk6GL.beagle.s1.log
	-> Dumping file: tskSim/tsk6GL.beagle.s1.filter
Input: lname=tskSim/tsk6GL.beagle nPop=3, fname=(null) qname=(null) outfiles=tskSim/tsk6GL.beagle.s1
Setup: seed=1 nThreads=10 method=1
Convergence: maxIter=2000 tol=0.000000 tolLike50=0.010000 dymBound=0
Filters: misTol=0.050000 minMaf=0.000000 minLrt=0.000000 minInd=0
Input file has dim: nsites=100000 nind=75
Input file has dim (AFTER filtering): nsites=100000 nind=75
iter[start] like is=9299805.984931
iter[50] like is=-6531138.892608 thres=0.002800
iter[100] like is=-6528710.773349 thres=0.001289
iter[150] like is=-6528405.896951 thres=0.001211
iter[200] like is=-6528306.803820 thres=0.000420
iter[250] like is=-6528277.160993 thres=0.000546
iter[300] like is=-6528271.925055 thres=0.000033
iter[350] like is=-6528271.177692 thres=0.000008
iter[400] like is=-6528270.876315 thres=0.000005
iter[450] like is=-6528270.772894 thres=0.000140
iter[500] like is=-6528270.747721 thres=0.000002
iter[550] like is=-6528270.740654 thres=0.000002
Convergence achived because log likelihooditer difference for 50 iteraction is less than 0.010000
best like=-6528270.740654 after 550 iterations
	-> Dumping file: tskSim/tsk6GL.beagle.s1.qopt
	-> Dumping file: tskSim/tsk6GL.beagle.s1.fopt.gz
	[ALL done] cpu-time used =  671.82 sec
	[ALL done] walltime used =  114.00 sec

fopt file

Contents of the fopt file

...
0.75331646167520038837 0.51190946588401886608 0.50134051056701267601 
0.99999999900000002828 0.80165850924934911603 0.97470665326916294813 
0.99999999900000002828 0.89560828888972687789 0.88062641752218895341 
0.99999999900000002828 0.99999999900000002828 0.86109994249930577048 
0.70560445653074521655 0.78994686954000448154 0.93076614062025020413 
0.99999999900000002828 0.88878537780630872955 0.92662857068149151463 
0.05322676762098016434 0.22871739860812340117 0.17394852600322696645 
0.00000000100000000000 0.27428885137150410545 0.19029599645013275944 
0.57086006389212373691 0.42232596591112880891 0.74080063581586474974 
0.77359733910003525281 0.47380864146016693494 0.72073560889718923939 
0.49946404159405927148 0.21684946347150244050 0.15201985942558055021 
0.41802171086717271331 0.55490556205954566504 0.85691127728452165524 
0.77095213528720529794 0.60074618451005279418 0.70219544996184157792 
0.26517850405564091787 0.48500265408436060710 0.85432254709914456914 
0.80055081986260245852 0.74423201242010783574 0.87110476762969968334 
0.30563054476851375663 0.05233529475348827620 0.25911912824038613179 
0.51084997710733415222 0.62263692178557350498 0.50738250264097506381 
0.64790272562679740442 0.91230541484222271720 0.73015721390331478347 
0.07124629651164265942 0.37896482494356753534 0.29218012479334326548 
0.00000000100000000000 0.26969100790961914038 0.28395781874856029781 
0.97074775756045073027 0.79093498372643300520 0.64006920058897498471 
0.64661948716978157048 0.84130009558421925409 0.76730057769159087933 
0.86990900887920663553 0.79410745692063922085 0.69416721874359499367 
0.34956069940263900797 0.27773038429396151860 0.25923476721423144298 
0.77739744690560164120 0.51272232330145017798 0.53888718200036844763 
0.35431569298041332150 0.20022780744715171219 0.43176580786072032980 
0.91858160919413811563 0.99999999900000002828 0.93584179237779097082 
0.90339823126358831384 0.94729687041528465308 0.84358671720630329371 
0.87068129661127857677 0.65267891763324525911 0.59315740612546075106 
0.24102496839012735319 0.42777100607917967201 0.39594098602469629533 
0.99999999900000002828 0.99999999900000002828 0.78549330115836857313 
0.15386277372522660922 0.18035502891341426146 0.26583557049163752950 
0.22456748943597096280 0.25110807159057474403 0.17244618960511531869 
0.74816053649164548922 0.54769319158907958656 0.44532166240679449398 
0.76350303696805599252 0.86547244122202959815 0.94111974586621383043 
0.40940400475566068872 0.67767095908245833513 0.40793761498610620064 
0.85389765162910868934 0.78901563183853873351 0.93614065916219291186 
0.54108661985898742763 0.61895909938546000983 0.88522763262549941654 
0.99051495581855464323 0.78855843624128341141 0.77646441702623147929 
0.51133721761171413434 0.74521610846562824637 0.32689774480116673416 
0.66618479413060949224 0.67891474309775079465 0.80762116232856140385 
0.81793598261160704865 0.77752326447671193943 0.95349025244041396565 
0.82120324647844433752 0.99999999900000002828 0.89800731971059466474 
...

There is currently no way to know the position of the lines of the fopt file if some sites have been flltered from the analysis (-minMaf, minInd, minLRT etc)

qopt file

Contents of the qopt file # cat tsk48GL.beagle.gz.s1.qopt

0.00254460532103031574 0.00108987228478324210 0.99636552239418640919 
0.00000015905647541105 0.00000000100000000000 0.99999983994352459327 
0.00034770382567266174 0.02639209238328452459 0.97326020379104283275 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00000467398081877176 0.00000000100000000000 0.99999532501918120264 
0.00000000907496942853 0.00585150933779484805 0.99414848158723567728 
0.00515826525767644137 0.01138897436535154552 0.98345276037697204607 
0.03914841746468285949 0.00000000100000000000 0.96085158153531713410 
0.00000000100000000000 0.00629199375758324100 0.99370800524241675866 
0.00771173022930659625 0.00000154720357311662 0.99228672256712036059 
0.00000000100000000000 0.00075135345721917719 0.99924864554278081119 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000005468413042120 0.00087279924180633879 0.99912714607406327705 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00712941313019542066 0.00118955677574110528 0.99168103009406338710 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00165385222968000606 0.99834614677032007535 
0.00000000100000000000 0.00006297763597355473 0.99993702136402651259 
0.00519087111391381209 0.00000000100000000000 0.99480912788608621966 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00202872783596746379 0.00000000100000000000 0.99797127116403261393 
0.00876424336999809782 0.00949457841911990376 0.98174117821088191516 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.01820430093358888640 0.00000694033297829119 0.98178875873343274261 
0.00351013812443964728 0.00000020340562512923 0.99648965846993520223 
0.00771897550085272680 0.00605259705033356268 0.98622842744881378252 
0.00600595292580561029 0.00000000100000000000 0.99399404607419439284 
0.01454910070242997067 0.00543457657939076105 0.98001632271817917808 
0.02567862615486414535 0.00160921436783232220 0.97271215947730349516 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00000000100000000000 0.00001041560507852223 0.99998958339492149960 
0.00000000100000000000 0.01383432553657116572 0.98616567346342876021 
0.00343840097404925389 0.00000000100000000000 0.99656159802595079000 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00051244065751142103 0.00404846039501185508 0.99543909894747661937 
0.02003953974792894652 0.00000004934009128878 0.97996041091197982897 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00000000100000000000 0.00000000100000000000 0.99999999799999994554 
0.02176809890633762956 0.00000000100000000000 0.97823190009366245423 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.01563096189267457192 0.00970868396771427770 0.97466035413961116252 
0.00000000100000000000 0.00000000100000000000 0.99999999800000005656 
0.00002540964943070735 0.00000000100000000000 0.99997458935056915408 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99501476026684787524 0.00000000100000000000 0.00498523873315206718 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99520671498720802983 0.00479241730266987201 0.00000086771012207898 
0.95884374919730619435 0.00000000100000000000 0.04115624980269377842 
0.99002104218586972628 0.00000000100000000000 0.00997895681413022567 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99999999770925251941 0.00000000129074746013 0.00000000100000000000 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.98980053177767901573 0.00000005577971952226 0.01019941244260143612 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99999785004878083416 0.00000000100000000000 0.00000214895121910354 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99220030909132039820 0.00000000100000000000 0.00779968990867968733 
0.99999996788621803301 0.00000000100000000000 0.00000003111378189772 
0.99736783433174225344 0.00255940950853666971 0.00007275615972113173 
0.99998096423035520708 0.00000000574461213317 0.00001903002503262207 
0.99711097909957713270 0.00288887008493822353 0.00000015081548462101 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99769262012085335734 0.00000000100000000000 0.00230737887914652393 
0.99999820787375570674 0.00000000433914936351 0.00000178778709493472 
0.98047422489554170166 0.00012980111977614777 0.01939597398468214523 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.98208006049140339488 0.00000000100000000000 0.01791993850859651197 
0.97530298545159921364 0.00000000100000000000 0.02469701354840085974 
0.99657542812406740840 0.00000000100000000000 0.00342457087593254226 
0.99954556420189066834 0.00045443479810919004 0.00000000100000000000 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99531584565237773976 0.00410740812985130408 0.00057674621777084644 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.99878572597704817770 0.00000000100000000000 0.00121427302295177490 
0.98571687209123504125 0.00400077401169816448 0.01028235389706666329 
0.99027397554762419674 0.00840892511494516215 0.00131709933743062008 
0.99999993504923445631 0.00000000100000000000 0.00000006395076564386 
0.95946639819101930957 0.00000000100000000000 0.04053360080898076034 
0.99999999800000005656 0.00000000100000000000 0.00000000100000000000 
0.98414939425022363029 0.01585059024074651421 0.00000001550902978739 
0.99999999622245250297 0.00000000277754757396 0.00000000100000000000 
0.99525652466242930938 0.00000001683386219288 0.00474345850370842034 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999999799999994554 0.00000000100000000000 0.00000000100000000000 
0.99999965447943561792 0.00000000100000000000 0.00000034452056438734 
0.99864814059528783652 0.00135185840471215468 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000001076370464123 0.99999998823629543399 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000100000000000 0.99999999800000005656 0.00000000100000000000 
0.00000000100000000000 0.99999986659623718577 0.00000013240376283687 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000099999999999 0.99999999799999994554 0.00000000099999999999 
0.00000000100000000000 0.99632783404679736705 0.00367216495320256799 
0.00000000100000000000 0.99999999799999994554 0.00000000100000000000 
0.00000000099999999999 0.99999999800000005656 0.00000000099999999999 
0.35919621347731411909 0.32381633362411937904 0.31698745289856661289 
0.31048363757756514136 0.30902410742704566893 0.38049225499538924522 
0.36341140678787386964 0.33678307361394943520 0.29980551959817652863 
0.34550713774447228133 0.34037087985425079628 0.31412198240127681137 
0.34705579219215104692 0.35218792485566730033 0.30075628295218165276 
0.33646039412306782967 0.32632754139618752598 0.33721206448074481088 
0.31881401220765009930 0.34885621407165418040 0.33232977372069577582 
0.34999374672052624424 0.33030931848049555066 0.31969693479897826061 
0.33152251818028721786 0.32339147992992234304 0.34508600188979043910 
0.31959998197389311025 0.33152491237148390413 0.34887510565462298562 
0.34724548642936803322 0.31809475756470984020 0.33465975600592196004 
0.33378069767858009609 0.33223636639277298599 0.33398293592864686241 
0.32023090400419051971 0.33179989332826043125 0.34796920266754916007 
0.35205158009776410521 0.33547091017851976558 0.31247750972371612921 
0.34291063455495451873 0.31853488093100223999 0.33855448451404313026 
0.31929132670383747472 0.32755905579808902717 0.35314961749807355362 
0.34114474726121107873 0.34607065583774476725 0.31278459690104404300 
0.33725705347681012025 0.32910919226619778089 0.33363375425699209886 
0.33918213722968154622 0.32278745806952213737 0.33803040470079642743 
0.33788659799509024317 0.34692305448657090317 0.31519034751833896468 
0.35876135180876589370 0.33843260979944000955 0.30280603839179404124 
0.34721570614318736370 0.34395335873604998556 0.30883093512076259524 
0.34165097731337079612 0.32814110943000784903 0.33020791325662146587 
0.33922542743931027864 0.32639619830977489867 0.33437837425091476717 
0.34461619391735059947 0.33133174331942943924 0.32405206276321996128 
0.34277551565686120716 0.32746953398981676342 0.32975495035332202942 
0.33842982221926010133 0.31224638933762871584 0.34932378844311123833 
0.34443810815667752490 0.32640113997211872565 0.32916075187120380496 
0.31723258569943768581 0.34955203711397470068 0.33321537718658750249 
0.35394053250677920408 0.33291498389624818444 0.31314448359697255597 
0.33504517457864940733 0.34188143503173562543 0.32307339038961496724 
0.33240938202788244960 0.34671459781042585080 0.32087602016169164409 
0.31745792352948248860 0.33722730677636020280 0.34531476969415725309 
0.33098224522913716195 0.33312298285105168549 0.33589477191981131909 
0.34090909280056919117 0.32423671881295645925 0.33485418838647434958 
0.32985465610121944557 0.32124851771265583444 0.34889682618612483100 
0.33525528582568764335 0.31967441393853385234 0.34507030023577844879 
0.33823045943274382408 0.33932114218381809190 0.32244839838343819505 
0.34374166546335593875 0.33527470302709477812 0.32098363150954922762 
0.32177399566214615056 0.34277626859597382092 0.33544973574188002852 
0.34915111840878915173 0.33072079898488659921 0.32012808260632419355 
0.31132788816691708833 0.32844185942225745389 0.36023025241082540227 
0.33067206673512555826 0.34601992411426535368 0.32330800915060908807 
0.31337643746173032833 0.33835721859074846529 0.34826634394752131740 
0.32762993090356395953 0.34856645453438306337 0.32380361456205303261 
0.33558678075595765877 0.34449062515269568419 0.31992259409134682357 
0.33433652456352996873 0.32868556951924504661 0.33697790591722515119 
0.32115036446030281736 0.35050069566489522321 0.32834893987480190392 
0.32524569843140932468 0.33953480032298033464 0.33521950124561045170 
0.33520046917246110185 0.31124301814705779279 0.35355651268048110536 
0.51565151014669796670 0.00027180960956305278 0.48407668024373901039 
0.51978922685130035664 0.01333903580405943964 0.46687173734464021413 
0.48123878312258933088 0.00648941795451128591 0.51227179892289931296 
0.48941833241028537271 0.00512373007237581363 0.50545793751733880672 
0.48421136927686320162 0.00600153379448644612 0.50978709692865020742 
0.53246468447754891073 0.00000000100000000000 0.46753531452245111755 
0.50637710620505416159 0.01564455874020675985 0.47797833505473913407 
0.49416813414210103428 0.00000000100000000000 0.50583186485789899400 
0.51328206693115174808 0.00000000100000000000 0.48671793206884833571 
0.50420356848059588728 0.00779539942445491366 0.48800103209494921641 
0.51589943710654184716 0.00000000100000000000 0.48410056189345807010 
0.46643393286795947761 0.00024627960390510270 0.53331978752813535838 
0.50134326603627110686 0.00000000100000000000 0.49865673296372897694 
0.52516062216154979492 0.00887494007947397384 0.46596443775897622430 
0.50553300231497877437 0.00610541400596737328 0.48836158367905380118 
0.48505848053244243756 0.00412236953776635561 0.51081914992979127188 
0.50419106430093152404 0.00671707921410998055 0.48909185648495850929 
0.51266037905765671212 0.00565931340437971983 0.48168030753796364785 
0.50479638826213368841 0.00082364200405335279 0.49437996973381287402 
0.48963785250324892706 0.00000000100000000000 0.51036214649675115673 
0.49861342640726780129 0.00000000100000000000 0.50138657259273211597 
0.49321745088202589846 0.00000000100000000000 0.50678254811797418533 
0.52297921048641760056 0.00000000100000000000 0.47702078851358242773 
0.51351947193443381323 0.00000000100000000000 0.48648052706556610403 
0.49861600587139209839 0.01143470350387426789 0.48994929062473369097 
0.47497824395255133778 0.00413641430709298184 0.52088534174035572288 
0.50602874958787047444 0.00000013752429825494 0.49397111288783129845 
0.51347175918678078510 0.00477133273041653854 0.48175690808280269284 
0.50359809216181616875 0.00000002299679746021 0.49640188484138642044 
0.52201190781479689385 0.00000000100000000000 0.47798809118520296790 
0.52427554763933403859 0.01637369304678280152 0.45935075931388308357 
0.50464335890649447691 0.01062810063722730188 0.48472854045627822295 
0.48795095623978190780 0.00032508303858300066 0.51172396072163517378 
0.49273360783177866384 0.03185613233234574349 0.47541025983587564818 
0.49075081269029041664 0.00043182816413278401 0.50881735914557668643 
0.51236233643387329995 0.01050799870797843559 0.47712966485814828355 
0.51939186110717183720 0.00638063180499700081 0.47422750708783106832 
0.49685157861691658931 0.00000000100000000000 0.50314842038308338346 
0.50376251978896124939 0.00609062514993390959 0.49014685506110500235 
0.50469879197514677660 0.00000000100000000000 0.49530120702485330719 
0.48806858812981018803 0.00000000100000000000 0.51193141087018978475 
0.49345173654735252633 0.00767168036095551131 0.49887658309169191639 
0.51926063211476558568 0.00000000100000000000 0.48073936688523438709 
0.49182360714466144547 0.00000000100000000000 0.50817639185533869384 
0.50012065040991493525 0.00101172020552988784 0.49886762938455525562 
0.49490771372946151807 0.00000000100000000000 0.50509228527053839919 
0.50981594186492362741 0.01168450085559137597 0.47849955727948501050 
0.48459184220397827358 0.00000007440008454733 0.51540808339593724430 
0.51153925961371649045 0.00045999176804108893 0.48800074861824249695 
0.49380129779182529992 0.00214174101547949525 0.50405696119269527422 
0.10504303642339951619 0.45848347542219436423 0.43647348815440606407 
0.09383999674587484296 0.44580529318052469767 0.46035471007360045936 
0.11801124345951279071 0.44619343422410290279 0.43579532231638429263 
0.10150817897299509174 0.44474184109029252232 0.45374997993671234431 
0.14144944553914898244 0.47426718065022838156 0.38428337381062249722 
0.08656596263718574491 0.47201374694852676894 0.44142029041428754166 
0.10422682420288104099 0.45665008652196642513 0.43912308927515242285 
0.07422281507005458467 0.46668026430253822801 0.45909692062740725671 
0.11152984148911383733 0.44326164444242566187 0.44520851406846068121 
0.12101900721666984662 0.45534926548479054409 0.42363172729853953991 
0.19287147372937366030 0.40220634979635128126 0.40492217647427497518 
0.19868166550667537562 0.39952077624337684059 0.40179755824994778379 
0.20144056442189406386 0.40552701281654912613 0.39303242276155692103 
0.17400131741109717276 0.41572345587205422612 0.41027522671684846234 
0.19363830614785534912 0.39941552029693161430 0.40694617355521295332 
0.20932370419936904837 0.41063785306931777086 0.38003844273131326403 
0.21496306930156286463 0.41077627378883840858 0.37426065690959875454 
0.20887311245081657818 0.39219787302656328176 0.39892901452262014006 
0.18789467459437667052 0.42880445734573224836 0.38330086805989094234 
0.21467435158258502126 0.41396326091136687042 0.37136238750604805281 
0.30215275924600598634 0.35114326369103593395 0.34670397706295807971 
0.27985580964526363124 0.36766711333486662427 0.35247707701986974449 
0.29214764907998119758 0.34353124024041165052 0.36432111067960715189 
0.28098186396660507214 0.35436535705487937076 0.36465277897851555711 
0.29909659519210785028 0.34708664349540557792 0.35381676131248662731 
0.29960230758566036569 0.34764467237891033546 0.35275302003542929885 
0.28690707484319816212 0.36958476358894237768 0.34350816156785934918 
0.31218824558522878521 0.35988855578362860532 0.32792319863114272049 
0.29371283648699086921 0.34536893102077848017 0.36091823249223065062 
0.32028624797598659324 0.35059182523172049972 0.32912192679229296255 
0.39315538655109805166 0.30778919233772789044 0.29905542111117405790 
0.39625700997625840083 0.29350948690034872612 0.31023350312339292856 
0.40087160410050781678 0.31851581382017457589 0.28061258207931755182 
0.40117357253398744366 0.30569836130272198815 0.29312806616329067921 
0.40013703551439627759 0.28691859513594913933 0.31294436934965452757 
0.39131222513930874474 0.30759794867682349606 0.30108982618386764818 
0.40826221599444090238 0.30658973748486684219 0.28514804652069231095 
0.41420080477834714250 0.28227625784283560950 0.30352293737881719249 
0.39119930707342420728 0.32102763805993583812 0.28777305486664006562 
0.37635520411942069430 0.29805329179310008358 0.32559150408747933314 
0.51400585200303006150 0.26100245041580294458 0.22499169758116702167 
0.50336119658518030384 0.25110166586697690860 0.24553713754784287082 
0.47299237773462793344 0.26084178003823194070 0.26616584222714018138 
0.49359314224598493936 0.26013978211456978418 0.24626707563944530421 
0.52795469779405246324 0.26499345968140075591 0.20705184252454675309 
0.48219467330650939152 0.25987283477635270135 0.25793249191713785162 
0.47626160019217189667 0.25351817092177358903 0.27022022888605456981 
0.51617477226059282902 0.23162353057460718930 0.25220169716479995392 
0.49698887507445854705 0.24557159475841641716 0.25743953016712495252 
0.52733914260860248469 0.25309832534801629533 0.21956253204338116447 
0.56749881833694781896 0.19172441472755546998 0.24077676693549673881 
0.59339160859286765870 0.19241414198845174788 0.21419424941868048240 
0.62308540846251914136 0.18054125203843729430 0.19637333949904353658 
0.59485531592769125275 0.20909554531024135415 0.19604913876206744860 
0.61310545246842529377 0.20645329445333451823 0.18044125307824007698 
0.60102956519838679483 0.21237444166376903687 0.18659599313784405727 
0.59278179178128642679 0.20826418834431797977 0.19895401987439562119 
0.60456224253100432353 0.20686687908046738626 0.18857087838852840123 
0.59417710257213784963 0.21264514488765640099 0.19317775254020574938 
0.59059286756608764257 0.21451811369415349495 0.19488901873975889023 
0.69484036887292865980 0.14634823390637874407 0.15881139722069256837 
0.69945423984127830241 0.16333221995631252987 0.13721354020240922322 
0.69115689116107958956 0.14927316115273414621 0.15956994768618620872 
0.68851717088680941536 0.14201541767923545057 0.16946741143395496754 
0.69288781352263861812 0.14270021794166909412 0.16441196853569234326 
0.68819873910998985433 0.16242980538224471854 0.14937145550776548264 
0.68619763716276405141 0.14370194479775053042 0.17010041803948539041 
0.68596343194490616568 0.16051691534743553480 0.15351965270765843830 
0.70684340251150390433 0.16654037983665334610 0.12661621765184280508 
0.70657158115262697073 0.14984891346689468983 0.14357950538047842270 
0.79161214498168253062 0.10430887542937690438 0.10407897958894059276 
0.79477141808375573184 0.10274451187208989700 0.10248407004415439892 
0.80425538032447896342 0.10720945367236509038 0.08853516600315590457 
0.79445836435866723502 0.11481368508653701233 0.09072795055479568327 
0.80626524450581027459 0.08599284906042292675 0.10774190643376663212 
0.77991736902186048486 0.08777798585427237787 0.13230464512386716502 
0.77897241390666871474 0.11419808069913564563 0.10682950539419577840 
0.80225596727756287585 0.10739115862914316857 0.09035287409329402497 
0.81035643868218754093 0.11405964018980654928 0.07558392112800596530 
0.80474324803558927588 0.09992219310105134034 0.09533455886335934215 
0.89147290804053958002 0.05818869713285088757 0.05033839482660958098 
0.87135519951168793895 0.04885203404408157424 0.07979276644423052844 
0.90273220877706750187 0.05642671780738096193 0.04084107341555152232 
0.90299890240805003039 0.05982401615206547896 0.03717708143988454617 
0.88622329583732417646 0.03227381365259313073 0.08150289051008267893 
0.89149278212958615875 0.03556871666107842139 0.07293850120933542680 
0.90540444756330573650 0.06637446770308205735 0.02822108473361228942 
0.89581315874618450135 0.06675457610008654619 0.03743226515372900798 
0.86941364504212315101 0.03330392614486758773 0.09728242881300920575 
0.88098981477392690476 0.04673780362475228600 0.07227238160132080924 

Plot results

Use R

admix<-t(as.matrix(read.table("tsk48GL.beagle.gz.s1.qopt")))
barplot(admix,col=1:3,space=0,border=NA,xlab="Individuals",ylab="admixture")

citation

log

  • v32 june 25-2013; modified code such that it now compiles on OSX
  • v31 june 24-2013; First public version.