[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