Message-ID: <F3396991-2B77-4A3C-830C-976DCD8EEF50@gmail.com>
Date: 2016-11-21T19:29:15Z
From: Zhian Kamvar
Subject: Bug in ade4 randtest.amova and solution
Hello,
I am the maintainer for the poppr package, which depends on adegenet (and thus ade4). One of my users noticed a bug in the output of randtest.amova: https://groups.google.com/forum/#!topic/poppr/D1gpqgQM2F0
Basically, the alternate hypothesis for the randomization tests would be incorrect at every fourth hierarchical level. I contacted the maintainer and this bug has been fixed on the ade4 github page: https://github.com/sdray/ade4/pull/4
If you do run into this problem, there is a workaround by passing the result to as.krandtest()
n <- 3 # (set this to the number of hierarchical levels)
with(randomization_result, as.krandtest(sim, obs, alter = c("less", rep("greater", n)), call = call, names = names)
Cheers,
Zhian N. Kamvar