Help making ADMIXTURE output digestible
1
2
Entering edit mode
8.7 years ago
devenvyas ▴ 740

I am new to the program ADMIXTURE.

I am running it from K=1 to K=20 on a set of 2033 individuals and 276,899 SNPs post-LD-pruning. The sample order is based on the default order in the dataset I downloaded, so population members are not necessarily back-to-back and related populations are not necessarily listed near each other.

Ideally, in the morning I will have a set of 40 *.Q files (20 with CV and bootstrap and 20 without), so in preparation I am trying to figure out how to graph them with the bars in a meaningful order. I was wondering if anyone had suggestions on how to this (or general suggestions about managing ADMIXTURE output).

Thanks!

admixture • 6.0k views
ADD COMMENT
0
Entering edit mode

I have the exact same question! Here is a similar one: Can't figure out order of output rows in ADMIXTURE

ADD REPLY
0
Entering edit mode
7.3 years ago
beausoleilmo ▴ 580

I found this website explaining how to do the full analysis. Basically, you have to look in the .fam file. You can import this file in R like this: fam = read.table("~/path-to-the-file/my_analysis.fam"). After that I was reusing the first column and rename the rows of your .Q file: tbl = read.table("~/path-to-the-file/my_analysis.Q"). Use this dataset to print your barplot!

ADD COMMENT
0
Entering edit mode

Hello i was looking for the same question. I found multiple unanswered posts over this. Thank you for your answer but i am unable to reach the link provided by you 'this website'. you may please guide again. Thank you

ADD REPLY
1
Entering edit mode

Hi, I am also using ADMIXTURE and found this marvelous tutorial that includes an Rscript that will plot the admixture output. The only thing you need to provide is a two column file including in the first column your sample IDs (in the same order that are found in the .fam file) and second column the population. The script will sort the samples by population and plot them in the columns and show different K in the rows. Hope it helps!

ADD REPLY
1
Entering edit mode

You can also check out pong (README here), which creates interactive plots with the clusters and colors automatically aligned across K, groups redundant outputs (which is useful when you have multiple *.Q files per _K_), and highlights multimodality in the data.

ADD REPLY
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Thank you. I will check it out

ADD REPLY

Login before adding your answer.

Traffic: 2931 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6