[adegenet-forum] amova on genind object (using poppr)
Vikram Chhatre
crypticlineage at gmail.com
Tue Apr 7 23:45:32 CEST 2015
I am trying to do amova on individual loci within a genind object created
from a structure formatted data set.
Steps I have followed so far:
mygenind <- import2genind(.....)
myhier <- read.table('myhier.txt', header=T)
head(myhier)
Pop Reg
1 1 1
2 1 1
3 2 1
4 2 1
5 3 2
6 3 2
7 4 2
8 4 2
9 4 2
mygenclone <- as.genclone(mygenind)
sethierarchy(mygenclone) <- myhier
mygenclone_seploc <- seploc(mygenclone, truenames=TRUE, res.type('genind',
'matrix'))
amova_res <- poppr.amova(mygenclone_seploc, ~Reg/Pop, within=TRUE)
Error: mygenclone_seploc must be a genind object.
------------
Has anyone seen this error before? If anyone knows a workaround, I would
appreciate hearing about it. The statistic works fine for 'over all loci',
but I need to obtain per locus estimates.
Thanks for any feedback.
Vikram
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/adegenet-forum/attachments/20150407/5ca3991e/attachment.html>
More information about the adegenet-forum
mailing list