Comment ajuster les modèles mixtes pour les données de rapport?

Question

Je veux connaître la variance dans le rapport des arbres adultes aux semis, gaules aux semis chez les espèces. Comme la valeur de ratio va de 0 à 1, une façon de résoudre cela est d'adapter un LMM en transformant les données de réponse à la normale.

Si je ne transforme pas les données, y a-t-il un autre moyen de résoudre ce problème? À ma connaissance, le GLMM binomial peut contenir des données proportionnelles. Puis-je l'utiliser pour ajuster les données du ratio, à savoir, glmer (ratio ~ espèces + (1 | sampleplot), famille = binomial, poids = ??, données)? Comment définir la valeur des poids?

Answer 1

Vous devez probablement quelque chose le long de la ligne de:

glmer(AdultTrees/Seedlings~species + (1|sampleplot), 
    weights=Seedlings, 
    data=Data,family=binomial) 

ne suis pas sûr de l'espèce d'effet fixe cependant. Vous mentionnez que vous êtes intéressé par la variance du succès des semis parmi les espèces, ce qui suggère que vous voulez que cela soit ajusté comme un effet aléatoire à la place. Si toutes les espèces sont évaluées sur toutes les parcelles, cela ressemblera à ceci:

glmer(AdultTrees/Seedlings~ (1|species) + (1|sampleplot), 
    weights=Seedlings, 
    data=Data,family=binomial) 

Qui est une conception d'effets aléatoires croisés. Difficile à dire cependant sans regarder les données.