Output Genotype file format used in PCA analyses

300 views
Skip to first unread message

jha...@nevada.unr.edu

unread,
Apr 25, 2016, 5:18:51 PM4/25/16
to Stacks
Hi All,

I have been searching the Stacks tutorials and I am unable to find information on how to export genotype files in a format used in principal components analyses.  Is stacks able to do this?  I have sequenced what is believed to be 3 populations and want to run a few PCAs on the data.  I tried genotypes function but it doesn't give out the formatted needed for these analyses.

ind1_id genotype_locus1 genotype_locus2 genotype_locus3 genotype_locus4 genotype_locus5 ...
ind2_id genotype_locus1 genotype_locus2 genotype_locus3 genotype_locus4 genotype_locus5 ...
ind3_id genotype_locus1 genotype_locus2 genotype_locus3 genotype_locus4 genotype_locus5 ...
ind4_id genotype_locus1 genotype_locus2 genotype_locus3 genotype_locus4 genotype_locus5 ...

where the first column is the ID and subsequent columns are the loci.  0=homozygous 1=heterzygous 2=homozygous

Any help would be greatly appreciated.

-Josh Hallas

Louis Plough

unread,
Apr 25, 2016, 10:37:21 PM4/25/16
to stacks...@googlegroups.com, Stacks
I would read the manual/info on the 'populations' module . You can export to any number of standard pop gen genotype formats. e.g. Genepop is a good one for downstream analyses like PCA. 

LP
--
Stacks website: http://catchenlab.life.illinois.edu/stacks/
---
You received this message because you are subscribed to the Google Groups "Stacks" group.
To unsubscribe from this group and stop receiving emails from it, send an email to stacks-users...@googlegroups.com.
Visit this group at https://groups.google.com/group/stacks-users.
For more options, visit https://groups.google.com/d/optout.

Claire Daguin Thiébaut

unread,
Apr 26, 2016, 2:43:14 AM4/26/16
to stacks...@googlegroups.com
Hi all,

you can also use the vcf file output from "populations". Then you can use the StackR package (https://github.com/thierrygosselin/stackr) to transform your vcf into a genind file with the vcf2genind function. The genind file will be use as an input for producing the PCA with Adegenet.

Cheers

Claire Daguin Thiebaut

Joshua Hallas

unread,
Apr 26, 2016, 10:46:02 AM4/26/16
to stacks...@googlegroups.com
Thanks for the help Claire. 

-josh

-----------------------------------------------------------------
Joshua M. Hallas
PhD. Student, Parchman Lab
Ecology, Evolution, and Conservation Biology
University of Nevada, Reno
1664 N. Virginia St.
Reno, Nevada 89557
https://jmhallasresearch.wordpress.com/
------------------------------------------------------------------

You received this message because you are subscribed to a topic in the Google Groups "Stacks" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/stacks-users/gp8dzpCMp7E/unsubscribe.
To unsubscribe from this group and all its topics, send an email to stacks-users...@googlegroups.com.
Reply all
Reply to author
Forward
0 new messages