Go Back   SEQanswers > Bioinformatics > Bioinformatics

Thread Tools
Old 11-03-2020, 11:45 AM   #1
Location: Germany

Join Date: Dec 2010
Posts: 29
Default Constrained ordination with SNP data

I have a genome-wide SNP dataset from RADseq containing 100 individuals from 4 groups. I converted the data so that it's in the form of a matrix (columns are SNPs, rows are individuals).
I'd like to use a constrained ordination (like a discriminant function) analysis in which I use the individuals from 2 of the 4 total groups to infer the axis that best separates these 2 groups across all SNPs. I then want to predict, along this axis, the position of the individuals of the other two groups.
I can easily convert the current SNP matrix into a numeric count matrix, such as by recoding, at each SNP, homozygous genotypes for allele A as 0, homozygous genotypes for the alternative allele B as 2, and heterozygous genotypes as 1.

My question is: what type of constrained ordination is appropriate for this type of data and goal? I'd like to do this in R.

Last edited by Marius; 11-03-2020 at 12:02 PM.
Marius is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 11:36 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO