If I am to run an ANOVA as such;
d3.totalcover.aov <- aov((Ecklonia + Phyllospora.comosa) ~ Location*Habitat, data=d3)
Where values of Ecklonia and Phyllospora are % covers of canopy forming kelp species (6 datapoints of each species within each Location*Habitat interaction). Is this a valid argument to make to observe total canopy cover. Or am I inadvertently affecting the mean total canopy cover?
I.e do I have to combine the total % coverage, defining it as Y in an ANOVA?
I tried:
total.cover <- (Ecklonia + Phyllospora, d3)
But wasn't a valid argument.
Cheers
P.S sorry if poorly worded, stats is my second language... 