I am trying to run an ANCOVA in R with abiotic factours as multiple independent variables and ecoregion as a covariate.
This is the code I had originally ran but I'm not sure if it accounts for ecoregion as a covariate.
aov(Quantity ~ temperature + dissolvedoxygen + specificconductance + discharge + ecoregion, data data)
I had originally though to use * instead of + between the different abiotic factours but that only looks at interactions if I'm not mistaken.
What would be the appropriate code to use with multiple independent variables and one covariate for an ANCOVA?
