Take the 2-minute tour ×
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It's 100% free, no registration required.

I make habitat suitability models for animal species. The purpose of my research is to investigate the accuracy of different models.

I clearly have a nested design:

  • accuracy_measure -> response variable
  • 2 model types (model_type) -> fixed effect
  • 230 species (species) -> random effect
  • 10 replicates/species (replicate) -> random effect
  • 10 subreplicate/replicate -> observation

So I have: $10\times 10\times 230$ observations/model, identified as speciesID-replicate-subreplicate (species_ID ranging 1:230, replicate 1:10 and subreplicate 1:10)

One could think about such mixed-effect model:

my.model <- lme(fixed = accuracy_measure~model_type, random = ~1|species/replicate)

I do not expand here into model simplification, yet... but here are my questions:

  1. my replicates & subreplicates are paired in the sense that they come from the same split of the data. As an example for species X, the 20 observations of replicate1 of model A and B (10 for A and 10 for B) are linked by a same data split which is likely to influence my accuracy measure. In the same way the 2 observations replicate1-subreplicate2 of model A and B (1 for A and 1 for B) are also linked. Is there a way to introduce such pairing in a mixed effect model?

  2. several continuous covariates, attributes of species (number of points for modelization, size in km2 of the range) may influence the measures of accuracy and I may be interested in investigating those effects. How could I include a covariate in such model? How should I strucutre it given that there are species covariates (high order of nesting)?

I hope that my questions are relevant!

share|improve this question

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.