Question: Replicates in Canonical Correspondence Analysis or PERMANOVA
gravatar for lelle
4.2 years ago by
lelle830 wrote:

My problem is the statistical analysis of community data with gradients of environmental data.

I have a table of species occurrences (from metabarcoding) per sample site and measurements of environmental parameters (like pH, nutrient concentration, ...). We did all our samples in technical triplicates. That is, we took each sample three times and ran it trough our complete lab process.

Now these replicates are obviously not statistically independent from each other so I have to take that into account when doing statistical test on the data.

I want to test if the environmental data we measured explain the variance of the metabarcoding data. For this I have so far considered Canonical Correspondence Analysis and PERMANOVA. If I understand what I read correctly, in a linear model I would have to include a "random effect" of the replicate. As far as I can see I can not include such an effect in the before mentioned tests.

I found in a FAQ of the vegan package the statement that this is not possible (at least in vegan) and as a work around they suggest to use partial CCA or include a "strata argument" in the PERMANOVA. I think this is not possible to use in my case, because this would imply that the effect of the replicate number is consistent over the sample sites.

The only other idea I have dome up with so fat is to average over the replicates, but I would rather find a different solution.

Any ideas or pointers would be greatly appreciated.

R metabarcoding vegan • 1.1k views
ADD COMMENTlink modified 4.2 years ago • written 4.2 years ago by lelle830

I have a feeling you'll have better luck with this question on cross validated.

ADD REPLYlink written 4.2 years ago by Devon Ryan97k

Thanks I will give it a try

ADD REPLYlink written 4.2 years ago by lelle830
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: 2278 users visited in the last hour