Je suis sur le point de plonger dans l'apprentissage R et mon projet d'apprentissage consistera à appliquer une régression à effets mixtes ou aléatoires à un ensemble de données afin de développer une équation prédictive. Je partage l'inquiétude de l'auteur dans ce post Comment choisir la bibliothèque nlme ou lme4 R pour les modèles d'effets mixtes? en me demandant si NLME ou LME4 est le meilleur package pour me familiariser. Une question plus fondamentale est: quelle est la différence entre la modélisation à effets mixtes linéaire et non linéaire?
Pour le fond, j'ai appliqué la modélisation ME dans ma recherche sur la SEP (dans MATLAB, pas R), donc je connais la façon dont les variables fixes vs aléatoires sont traitées. Mais je ne sais pas si le travail que j'ai fait était considéré comme ME linéaire ou non linéaire. Est-ce simplement la forme fonctionnelle de l'équation utilisée, ou autre chose?
Réponses:
Il existe plusieurs distinctions entre les modèles de régression linéaire et non linéaire, mais la principale mathématique est que les modèles linéaires sont linéaires dans les paramètres, tandis que les modèles non linéaires sont non linéaires dans les paramètres. Pinheiro et Bates (2000, pp. 284-285), auteurs du
nlme
package R, ont décrit avec élégance les considérations les plus importantes dans la sélection des modèles:Il existe également de grandes différences entre les packages nlme et lme4 qui vont au-delà du problème de linéarité. Par exemple, en utilisant nlme, vous pouvez ajuster des modèles linéaires ou non linéaires et, pour l'un ou l'autre type, spécifier les structures de variance et de corrélation pour les erreurs intra-groupe (par exemple, autorégressives); lme4 ne peut pas faire ça. De plus, les effets aléatoires peuvent être fixes ou croisés dans les deux packages, mais il est beaucoup plus facile (et plus efficace en termes de calcul) de spécifier et de modéliser les effets aléatoires croisés dans lme4.
Je conseillerais d'abord de considérer a) si vous aurez besoin d'un modèle non linéaire, et b) si vous devrez spécifier la variance intra-groupe ou les structures de corrélation. Si l'une de ces réponses est oui, alors vous devez utiliser nlme (étant donné que vous vous en tenez à R). Si vous travaillez beaucoup avec des modèles linéaires qui ont des effets aléatoires croisés, ou des combinaisons compliquées d'effets aléatoires imbriqués et croisés, alors lme4 est probablement un meilleur choix. Vous devrez peut-être apprendre à utiliser les deux packages. J'ai d'abord appris lme4, puis j'ai réalisé que je devais utiliser nlme car je travaille presque toujours avec des structures d'erreur autorégressives. Cependant, je préfère toujours lme4 lorsque j'analyse les données d'expériences avec des facteurs croisés. La bonne nouvelle est qu'une grande partie de ce que j'ai appris sur lme4 a bien été transférée à nlme. D'une manière ou d'une autre,
Les références
Pinheiro, JC et Bates, DM (2000). Modèles à effets mixtes en S et S-PLUS . New York: Springer-Verlag.
la source
Pour la partie linéaire-non linéaire, voir: article CrossValidated sur le sujet , en particulier la réponse classée deuxième par Charlie. Je ne pense pas qu'il y ait de changements dans le traitement des effets mixtes.
la source