lme4 et nlme sont des packages R utilisés pour ajuster des modèles d'effets mixtes linéaires, linéaires généralisés et non linéaires. Pour les questions générales sur les modèles mixtes, utilisez la balise [mixed-model].
Permettez-moi d'abord de donner quelques informations; Je résumerai mes questions à la fin. La distribution bêta, paramétrée par sa moyenne et , a , où est la fonction de variance.μμ\muϕϕ\phiVar(Y)=V(μ)/(ϕ+1)Var(Y)=V(μ)/(ϕ+1)\operatorname{Var}(Y) = \operatorname{V}(\mu)/(\phi+1)V(μ)=μ(1−μ)V(μ)=μ(1−μ)\operatorname{V}(\mu) = \mu(1-\mu) Dans une régression bêta (par exemple, en utilisant le package betareg dans R), la régression suppose …
Quelle est la difference entre (1|DNA.concentration/mouse.id)et (DNA.concentration|mouse.id)? Que signifient les symboles |et la /signification dans la syntaxe de l'effet aléatoire?
Je suis curieux de savoir comment le package lmerTest dans R, en particulier la fonction "rand", gère les tests d'effets aléatoires. Prenons l'exemple du pdf lmerTest sur CRAN qui utilise l'ensemble de données "carottes" intégré: #import lme4 package and lmerTest package library(lmerTest) #lmer model with correlation between intercept and slopes …
J'ai un problème avec l'interprétation des interactions à 2 et 3 voies dans lmer. Mon DV est la hauteur qui est une variable continue. Tous les IV sont des variables catégorielles. Le premier facteur est animal, rat ou lion. Le deuxième facteur est le sexe, masculin ou féminin. Le troisième …
Mon projet actuel peut m'obliger à construire un modèle pour prédire le comportement d'un certain groupe de personnes. l'ensemble de données de formation ne contient que 6 variables (id est uniquement à des fins d'identification): id, age, income, gender, job category, monthly spend dans laquelle se monthly spendtrouve la variable …
Ma conception est la suivante. yyy est la réponse de Bernoulli x1x1x_1 est une variable continue x2x2x_2 est une variable catégorielle (facteur) à deux niveaux L'expérience est complètement dans les sujets. Autrement dit, chaque sujet reçoit chaque combinaison dex1x1x_1 et x2x2x_2. Il s'agit d'une configuration de régression logistique à mesures …
J'utilise le lme4package dans R pour faire de la modélisation logistique à effets mixtes. J'ai cru comprendre que la somme de chaque effet aléatoire devait être nulle. Lorsque je crée des modèles mixtes linéaires de jouets à l'aide de lmer, les effets aléatoires sont généralement < confirmant ma conviction que …
J'ai des problèmes à trouver une solution concernant la façon d'exécuter un test post-hoc (Tukey HSD) après une ANOVA à mesures répétées à 2 facteurs (tous deux intra-sujets) en R. Pour l'ANOVA, j'ai utilisé la fonction aov: summary(aov(dv ~ x1 * x2 + Error(subject/(x1*x2)), data=df1)) Après avoir lu les réponses …
J'ai 2 facteurs Aet B(5 × 3) et une covariable Xdans une conception intra-sujet. Voici comment je spécifie mon modèle global: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Mon interprétation est que je regarde un graphique y~x, où la pente change en raison de la covariable, et la ligne se déplace …
Mise à jour : Puisque je sais maintenant que mon problème est appelé séparation quasi-complète, j'ai mis à jour la question pour refléter cela (merci à Aaron). J'ai un ensemble de données provenant d'une expérience dans laquelle 29 participants humains (facteur code) ont travaillé sur un ensemble d'essais et la …
J'essaie d'obtenir des prédictions pour les observations d'un objet lme. C'est censé être assez simple. Pourtant, comme je reçois différents types d'erreurs pour différents essais, il me semble que je manque quelque chose. Mon modèle est le suivant: model <- lme(log(child_mortality) ~ as.factor(cluster)*time + my.new.time.one.transition.low.and.middle + ttd + maternal_educ+ log(IHME_id_gdppc) …
Supposons que A est un facteur fixe et B est un facteur aléatoire. A et B sont croisés. Un autre facteur aléatoire C est imbriqué dans B. Compte tenu des interactions AB et AC, comment dois-je adapter le modèle en utilisant la fonction lme dans le package nlme? Merci.
J'ai lu que la distribution de l'estimateur à effets aléatoires dans lme4 est très asymétrique et pour cette raison, les erreurs standard ne sont pas signalées. Je me demande si quelqu'un peut fournir une référence pour cela? J'ai accès au livre de Bates et Pinherio, mais pas Raudenbush et Bryk. …
Je travaille sur un ensemble de données afin d'évaluer l'impact du séchage sur les activités microbiennes des sédiments. L'objectif est de déterminer si l'impact du séchage varie selon les types de sédiments et / ou la profondeur dans les sédiments. Le plan expérimental est le suivant: Le premier facteur sédiment …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.