Question: How to do AMOVA on GBS data
gravatar for jaafari.omid
7 months ago by
jaafari.omid50 wrote:

Hello dears all, Actually I have a vcf file including five populations. This vcf file comes from GBS data and the stacks pipeline. I was going to calculate between and within population diversity based on AMOVA in R ( but I got an error which I couldn’t how can I solve that. Here is the command and the error which I have on my face. I will be so grateful if anybody can help me, either with command or any other protocols. Regards, Omid

poppr.amova(omid.genlight, hier = NULL, clonecorrect = FALSE, within = TRUE, dist = omid.D.pop.dist, squared = TRUE, freq = TRUE, correction = "quasieuclid", sep = "-", filter = FALSE, threshold = 0, algorithm = "farthest_neighbor", threads = 1L, missing = "loci", cutoff = 0.05, quiet = FALSE, method = c("ade4", "pegas"), nperm = 0)

Error in poppr.amova(omid.genlight, hier = NULL, clonecorrect = FALSE,: A population hierarchy must be specified

gbs snp R genome • 342 views
ADD COMMENTlink written 7 months ago by jaafari.omid50

I don't know anything about this R package or the specific function, but it looks like you need to supply an R object to hier. You would need to read the function's documentation. Perhaps the author(s) have a vignette or example data that you should try to see if it works with the example data.

ADD REPLYlink written 7 months ago by jean.elbers1.2k

Yes I should define an object but when I read it I couldn't understand what exactly I should to to have this object.

ADD REPLYlink written 7 months ago by jaafari.omid50

Presumably, you need to specify the hierarchy (structure) of your populations. Granted I haven't read the documentation for the R package or specific function.

ADD REPLYlink written 7 months ago by jean.elbers1.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1035 users visited in the last hour