Je vais décrire le modèle de chacun de vos appels lmer()
et leur différence, puis répondre à votre dernière question sur la sélection des effets aléatoires.
Chacun de vos trois modèles contient des effets fixes pour practice
, context
et l'interaction entre les deux. Les effets aléatoires diffèrent entre les modèles.
lmer(ERPindex ~ practice*context + (1|participants), data=base)
participants
participant
0
lmer(ERPindex ~ practice*context + (1+practice|participants), data=base)
Ce modèle, en plus d’une interception aléatoire, contient également une pente aléatoire dans practice
. Cela signifie que le rythme auquel les individus apprennent de la pratique est différent d'une personne à l'autre. Si une personne a un effet aléatoire positif, elle augmente plus vite avec la pratique que la moyenne, tandis qu'un effet aléatoire négatif indique qu'elle apprend moins vite avec la pratique que la moyenne, ou peut-être s'aggrave avec la pratique, en fonction de la variance de l'aléatoire. effet (en supposant que l'effet fixe de la pratique soit positif).
lmer(ERPindex ~ practice*context + (practice|participants) +
(practice|participants:context), data=base)
Ce modèle s'adapte à une pente aléatoire et intercepte dedans practice
(vous devez faire (practice-1|...)
pour supprimer l'interception), tout comme le modèle précédent, mais maintenant vous avez également ajouté une pente aléatoire et une interception dans le facteur participants:context
, qui est un nouveau facteur dont les niveaux Toutes les combinaisons des niveaux présents dans participants
et context
et les effets aléatoires correspondants sont partagées par des observations ayant la même valeur pour les deux participants
et context
. Pour adapter ce modèle , vous aurez besoin d'avoir plusieurs observations qui ont les mêmes valeurs pour les deux participants
etcontext
ou bien le modèle n'est pas estimable. Dans de nombreuses situations, les groupes créés par cette variable d'interaction sont très rares et donnent lieu à des modèles d'effets aléatoires très bruyants / difficiles à adapter. Vous devez donc faire attention lorsque vous utilisez un facteur d'interaction en tant que variable de regroupement.
Fondamentalement (lire: sans devenir trop compliqué), vous devez utiliser des effets aléatoires lorsque vous pensez que les variables de regroupement définissent des "poches" d’inhomogénéité dans l’ensemble de données ou que les individus partageant le même niveau de facteur de regroupement doivent être corrélés les uns aux autres ( les individus qui ne le sont pas ne devraient pas être corrélés) - les effets aléatoires y parviennent. Si vous pensez que des observations partageant les mêmes niveaux participants
et context
plus similaires que la somme des deux parties, il peut être approprié d’inclure l’effet aléatoire "interaction".
Edit: Comme @Henrik le mentionne dans les commentaires, les modèles que vous adaptez, par exemple:
lmer(ERPindex ~ practice*context + (1+practice|participants), data=base)
faire en sorte que la pente aléatoire et l'interception aléatoire soient corrélées l'une avec l'autre, et que cette corrélation soit estimée par le modèle. Pour contraindre le modèle de sorte que la pente aléatoire et l'interception aléatoire ne soient pas corrélées (et donc indépendantes, puisqu'elles sont normalement distribuées), vous devez plutôt adapter le modèle:
lmer(ERPindex ~ practice*context + (1|participants) + (practice-1|participants),
data=base)
Le choix entre ces deux options doit être fondé sur le fait que vous pensez, par exemple, que les participant
s ayant une ligne de base supérieure à la moyenne (c.-à-d. Une interception aléatoire positive) sont également susceptibles d'avoir un taux de changement supérieur à la moyenne (c.-à-d. Une pente aléatoire positive). Si tel est le cas, vous autoriseriez une corrélation entre les deux, sinon, vous les obligeriez à être indépendants. (Encore une fois, cet exemple suppose que la pente de l'effet fixe est positive).
lmer(ERPindex ~ practice*context + (1|participants) + (0 + practice|participants, data=base)
ou est-ce que je me trompe? (Unrelated: Excusez ma petite publication de votre message. Si vous êtes en désaccord avec cette précision, remettez-la en place.)x <-rnorm(1000); id <- rep(1:100,each=10); y <- rnorm(1000); g <- lmer(y ~ (1+x|id)); g2 <- lmer(y ~ (1|id) + (x-1|id)); attr(logLik(g),"df"); attr(logLik(g2),"df");
@ Macro a donné une bonne réponse ici, je veux juste ajouter un petit point. Si certaines personnes dans votre situation utilisent:
Je soupçonne qu'ils commettent une erreur. Considérez:
(practice|participants)
signifie qu’il existe une pente aléatoire (et une interception) pour l’effet depractice
pour chacunparticipant
, alors(practice|participants:context)
qu’il existe une pente aléatoire (et une interception) pour l’effet depractice
pour chaqueparticipant by context
combinaison . C'est bien, si c'est ce qu'ils veulent, mais je suppose qu'ils le souhaitent(practice:context|participants)
, ce qui signifie qu'il existe une pente aléatoire (et une interception) pour l'effet d'interaction depractice by context
chacunparticipant
.la source
Dans un modèle à effets aléatoires ou à effets mixtes, un effet aléatoire est utilisé lorsque vous souhaitez traiter l'effet observé comme s'il était issu d'une distribution de probabilité d'effets.
L'un des meilleurs exemples que je puisse donner concerne la modélisation des données d'essais cliniques issues d'un essai clinique multicentrique. Un effet de site est souvent modélisé comme un effet aléatoire. Cela est dû au fait que les quelque 20 sites réellement utilisés dans l’essai provenaient d’un groupe beaucoup plus vaste de sites potentiels. En pratique, la sélection n'a peut-être pas été faite au hasard, mais il peut toujours être utile de la traiter comme si elle l'était.
Bien que l'effet de site ait pu être modélisé comme un effet fixe, il serait difficile de généraliser les résultats à une population plus large si nous ne prenions pas en compte le fait que l'effet serait différent pour un ensemble sélectionné différent de 20 sites. Le traiter comme un effet aléatoire nous permet de l'expliquer de cette façon.
la source