Register for an account

X

Enter your name and email address below.

Your email address is used to log in and will not be shared or sold. Read our privacy policy.

X

Website access code

Enter your access code into the form field below.

If you are a Zinio, Nook, Kindle, Apple, or Google Play subscriber, you can enter your website access code to gain subscriber access. Your website access code is located in the upper right corner of the Table of Contents page of your digital edition.

Health

An informative ADMIXTURE plot (perhaps?)

Gene ExpressionBy Razib KhanJanuary 14, 2013 12:48 AM

Newsletter

Sign up for our email newsletter for the latest science news

In my earlier posts where I gave a short intro to using Plink I distributed a data set termed PHLYO. One thing I did not mention is that I've also been running it on Admixture. But here's an important point: I ran the data set 10 times from K = 2 to K = 15. Why? Because the algorithm produces somewhat different results on each run (if you use a different seed, which you should), and I wanted to not be biased by one particular result. Additionally, I also turned on cross-validation error, which gives me a better sense of which K's to trust. But after I select the K which I want to visualize which replicate run will I then use to generate the bar plots? I won't pick any specific one. Rather, I'll merge them together with an off-the-shelf algorithm. Additionally, I also want to sort the individuals by their modal population cluster. This sounds rather convoluted, and it is somewhat. I have a pipeline that I use, but it's not too user friendly. One of my projects is to clean it up, document it, and publish it online. Though if you have your own pipeline all ready to go, please post it in the comments with a link! The general steps are as follows for me: 1) Convert Admixture Q files into Structure format, transform family identifications to numeric values, and generate a file with family identification and numeral pairs 2) Merge the results across runs using Clumpp 3) Sort the individual results within populations 4) The use Distruct to produce an output file Before I show you the resultant bar plot, here are the cross-validation results with standard deviation ticks:

cv1.jpg

cv3.jpg

I'm going to select K = 10.

CORE.10.png

2 Free Articles Left

Want it all? Get unlimited access when you subscribe.

Subscribe

Already a subscriber? Register or Log In

Want unlimited access?

Subscribe today and save 70%

Subscribe

Already a subscriber? Register or Log In