NgsAdmix
This will contain the program called NGSadmix, which is a very nice tool for finding admixture. It is based on genotype likelihoods or genotype probabilities. It is a fancy multithreaded c/c++ program.
Latest version is 32 from June 25 2013. It can be found [1]
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, and that we want to use 4 computing cores.
./ngsAdmix312 -likes input.gz -K 3 -P 4 -tol 0.00000001 -tolLike50 0.01 -outfiles myoutfiles -minMaf 0.05
Input Files
Input files are contains genotype likelihoods in beagle input file format [2]. 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.
- PREFIX.log
- PREFIX.fopt.gz
- 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.