Hi Poornima,

As far as i can see you have:

1. Missing values in your dataset (NAs), which you should delete (or complete) them before you run the distances

2. You do not provide a structures dataset, which should be a DataFrame

please check the following code

########## R - code #############

samp

@poornima: AMOVA is analyse on molecular variance

@Kostas:

There isn't any missing value, I just proposed a tiny part on a 200x200matrix to give an aspect of it ;)

I did not provide structure because from ?amova:

"samples

a data frame with haplotypes (or genotypes) as rows, populations as columns and abundance as entries"

Then I should give abondance (then summing/getting frequency according to populations) that mean I only have 3 populations left without any structure to gather them …

hi Benjamin,

I apologise for any unintentional mistake with your name, due to lack of time, in my previous note.

I am not sure I understand your reasoning for the structures dataframe but it does not really matter, as it should run without it.

The only way I can reproduce your error is when the dataset is not in a "data.frame" object, but in a matrix

So, please provide us with the output of the command

class(YourDataSet)

Well, I had this error :p but I quickly figured that out with the help of the exemple:

> class(frec)

[1] "data.frame"

But I sure it is the way of doing my object that isn't good… I start from an object name "Betula" looking like that:

al270 al273 al121 al122 al123 al132 al134 al137 al146 al163.1 al217 al229

B-40-2x-T 0 0 0 1 1 0 0 0 1 0 0 0

B-41-2x-T 1 0 0 1 1 0 0 0 1 0 1 0

B-43-2x-T 0 0 1 1 1 0 0 0 1 0 0 0

B-47-2x-T 1 0 0 1 0 0 1 0 0 0 0 0

B-52-2x-T 1 0 0 1 1 0 0 0 0 0 0 0

B-72-3x-T 1 0 0 1 0 0 0 0 0 0 0 0

B-73-2x-A 1 0 0 1 1 0 0 0 1 0 1 0

B-75-4x-T 0 0 1 1 0 0 0 0 1 0 0 0

B-83-4x-A 0 0 1 1 1 1 0 0 1 0 0 0

(of course it is just a tiny random part of it)

Then I tried to prepare my sample object…

#prep data

apply(Betula[code$Ploidy==2,],2,sum)->num

nrow(Betula[code$Ploidy==2,])->dem

round(100*num/dem)->P2

apply(Betula[code$Ploidy==3,],2,sum)->num

nrow(Betula[code$Ploidy==3,])->dem

round(100*num/dem)->P3

apply(Betula[code$Ploidy==4,],2,sum)->num

nrow(Betula[code$Ploidy==4,])->dem

round(100*num/dem)->P4

frec

OK. Lets do it step by step

1. you need a data.frame to be you "samples" object. You are using the data.frame named "frec" which has a (?) number of rows (without row.names) and 3 columns named "V1","V2", "V3", with numerical values.

2. Then you are using a second object named "gel" to produce the distances (distE).

Does the "gel" object have the same number of rows and columns with "frec" data.frame??? If not the error is "Error in amova(frec, distE) : Non convenient samples"

Also you standardize it for presence-absence "binary=TRUE". Does this give you a second error like: "In is.euclid(distances) : Zero distance(s)" ????

Well… The "gel" object is either betula or t(betula) in function of the computer out put. In the present case, they are just identical.

number of row of "frec" is equal to number of row of betula (or gel) id est the number of marker…

Then distE is the euclidian distance between sample, so number of row= number of column = number of sample

I do not have any "In is.euclid(distances) : Zero distance(s)" like error…

hi,

vegdist requires a community matrix and not a data.frame.

Therefore the number of rows will not be the right one.

check the results of the following code:

################ R-code ######################

gel