Y a-t-il une raison intuitive pour que les effets aléatoires soient réduits à leur valeur attendue dans le modèle mixte linéaire
Structure de données > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2...
J'ai un GLMM du formulaire: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Lorsque j'utilise drop1(model, test="Chi"), j'obtiens des résultats différents de ceux que j'utilise à Anova(model, type="III")partir du package de voiture ou...
J'essaie d'intégrer un modèle mixte multivarié (c'est-à-dire à réponses multiples) R. Mis à part les packages ASReml-ret SabreR(qui nécessitent un logiciel externe), il semble que cela ne soit possible qu'en MCMCglmm. Dans l' article qui accompagne le MCMCglmmpaquet (p. 6), Jarrod Hadfield décrit...
Étant donné trois variables, yet xqui sont positives continues et zqui sont catégoriques, j'ai deux modèles candidats donnés par: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) et fit.fe <- lm( y ~ 1 + x ) J'espère comparer ces modèles pour déterminer quel modèle est le plus approprié....
Est-ce que quelqu'un sait comment calculer (ou extraire) l'effet de levier et les distances de Cook pour un merobjet de classe (obtenu via le lme4package)? Je voudrais les représenter pour une analyse des
Je comprends l'effet d'un effet aléatoire catégorique sur un modèle à effets mixtes en ce qu'il effectue un regroupement partiel des observations par niveau dans l'effet aléatoire, en supposant effectivement que les observations ne sont pas indépendantes elles-mêmes, mais seulement leurs pools...
Je travaille sur un ensemble de données. Après avoir utilisé certaines techniques d'identification de modèle, je suis sorti avec un modèle ARIMA (0,2,1). J'ai utilisé la detectIOfonction dans le package TSAen R pour détecter une valeur aberrante innovante (IO) à la 48e observation de mon ensemble...
On m'a toujours enseigné que les effets aléatoires n'influencent que la variance (erreur) et que les effets fixes n'influencent que la moyenne. Mais j'ai trouvé un exemple où les effets aléatoires influencent également la moyenne - l'estimation du coefficient: require(nlme) set.seed(128) n <-...
Il existe plusieurs discussions sur R-sig-ME sur l'obtention d'intervalles de confiance pour les prédictions utilisant lme4et nlmedans R. Par exemple ici et ici en 2010, y compris certains commentaires de Dougals Bates, l'un des auteurs des deux packages. J'hésite à le citer mot pour mot, de peur...
Cette question / sujet a été soulevé lors d'une discussion avec un collègue et je cherchais des opinions à ce sujet: Je modélise certaines données en utilisant une régression logistique à effets aléatoires, plus précisément une régression logistique à interception aléatoire. Pour les effets fixes,...
Dans le modèle mixte, nous supposons que les effets aléatoires (paramètres) sont des variables aléatoires qui suivent des distributions normales. Elle ressemble beaucoup à la méthode bayésienne, dans laquelle tous les paramètres sont supposés être aléatoires. Le modèle à effet aléatoire est-il donc...
J'analyse les résultats d'une expérience de temps de réaction chez R. J'ai effectué une ANOVA à mesures répétées (1 facteur intra-sujet avec 2 niveaux et 1 facteur inter-sujet avec 2 niveaux). J'ai exécuté un modèle mixte linéaire similaire et je voulais résumer les résultats de lmer sous la forme...
Quelle est la structure de variance-covariance par défaut pour les effets aléatoires dans glmerou lmerdans le lme4package? Comment spécifier une autre structure variance-covariance pour les effets aléatoires dans le code? Je n'ai trouvé aucune information à ce sujet dans la
J'analyse certaines données comportementales en utilisant lme4in R, principalement en suivant les excellents tutoriels de Bodo Winter , mais je ne comprends pas si je gère correctement les interactions. Pire encore, personne d'autre impliqué dans cette recherche n'utilise des modèles mixtes, donc...
J'ai le modèle suivant m_plotéquipé d' lme4::lmereffets aléatoires croisés pour les participants ( lfdn) et les éléments ( content): Random effects: Groups Name Variance Std.Dev. Corr lfdn (Intercept) 172.173 13.121 role1 62.351 7.896 0.03 inference1 24.640 4.964 0.08 -0.30 inference2 52.366 7.236...
Avant de poser cette question, j'ai fait une recherche sur notre site et j'ai trouvé beaucoup de questions similaires (comme ici , ici et ici ). Mais je pense que ces questions connexes n'ont pas été bien répondues ou discutées, je voudrais donc soulever à nouveau cette question. Je pense qu'il...
J'espère sincèrement avoir formulé cette question de manière à ce qu'elle puisse être résolue définitivement - sinon, faites-le moi savoir et je vais réessayer! Je devrais également, je suppose, noter que j'utiliserai R pour ces analyses. J'ai plusieurs mesures plant performance (Ys)dont je...
Je suis confus par les déclarations sur une page Web de l'UCLA sur la régression logistique à effets mixtes. Ils montrent un tableau des coefficients d'effets fixes de l'ajustement d'un tel modèle et le premier paragraphe ci-dessous semble interpréter les coefficients exactement comme une...
Envisagez une expérience avec plusieurs participants humains, chacun mesuré plusieurs fois dans deux conditions. Un modèle d'effets mixtes peut être formulé (en utilisant la syntaxe lme4 ) comme: fit = lmer( formula = measure ~ (1|participant) + condition ) Supposons maintenant que je souhaite...