Question: How to do AMOVA on GBS data
0
gravatar for jaafari.omid
8 weeks ago by
jaafari.omid10
jaafari.omid10 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 (https://rdrr.io/cran/poppr/man/poppr.amova.html) 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 • 174 views
ADD COMMENTlink written 8 weeks ago by jaafari.omid10

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 8 weeks ago by jean.elbers550

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 8 weeks ago by jaafari.omid10

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 8 weeks ago by jean.elbers550
Please log in to add an answer.

Help
Access

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