AMOVA coefficients

[Gabriella Scatà]

Hi everyone,

I just tested the AMOVa this week with the pegas implementation (as I understood that ade4 has a bug?).

This is what I obtain:

Analysis of Molecular Variance

Call: pegas::amova(formula = hier, data = hierdf, nperm = nperm, is.squared = FALSE)

                SSD        MSD  df
clusters 0.03572749 0.03572749   1
pop      0.05107500 0.01276875   4
Error    2.23568319 0.01330764 168
Total    2.32248567 0.01342477 173

Variance components:
              sigma2 P.value
clusters  2.9689e-04  0.0000
pop      -1.8593e-05  0.8808
Error     1.3308e-02        

Phi-statistics:
clusters.in.GLOBAL (Phi_CT)      pop.in.GLOBAL (Phi_ST)    pop.in.clusters (Phi_SC)
                0.021852940                 0.020484360                -0.001399156

Variance coefficients:
       a        b        c
28.98276 29.03448 77.33333

How should I interpret the coefficients?

I understand that they represent the slope/effect of each treatment (in this case group/population/individuals) on the genotype...but how would you interpret them along with the variances?

I would really appreciate some feedback

Thank you!

Best

Gabriella