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].
La question du CV J'essaie de donner (a) une représentation mathématique détaillée et concise d'un modèle à effets mixtes. J'utilise le lme4package dans R. Quelle est la représentation mathématique correcte pour mon modèle? Les données, la question scientifique et le code R Mon ensemble de données comprend des espèces dans …
J'expérimente l'algorithme de la machine de renforcement de gradient via le caretpackage en R. À l'aide d'un petit ensemble de données d'admission à l'université, j'ai exécuté le code suivant: library(caret) ### Load admissions dataset. ### mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") ### Create yes/no levels for admission. ### mydata$admit_factor[mydata$admit==0] <- "no" mydata$admit_factor[mydata$admit==1] <- …
Les tests de permutation (également appelés test de randomisation, test de re-randomisation ou test exact) sont très utiles et s'avèrent utiles lorsque l'hypothèse de distribution normale requise par exemple t-testn'est pas remplie et lorsque la transformation des valeurs par classement des un test non paramétrique comme Mann-Whitney-U-testcela entraînerait la perte …
Je voudrais faire correspondre les sorties de lmer (vraiment glmer) avec un exemple binomial jouet. J'ai lu les vignettes et je crois comprendre ce qui se passe. Mais apparemment non. Après être resté coincé, j'ai corrigé la "vérité" en termes d'effets aléatoires et je suis allé après l'estimation des effets …
J'ai un ensemble de données dans lequel la variable que j'aimerais utiliser comme effet aléatoire n'a qu'une seule observation pour certains niveaux. Sur la base des réponses aux questions précédentes, j'ai compris que, en principe, cela peut convenir. Puis-je adapter un modèle mixte avec des sujets qui n'ont qu'une seule …
J'examine actuellement certains travaux et suis tombé sur les suivants, ce qui me semble faux. Deux modèles mixtes sont montés (en R) à l'aide de lmer. Les modèles ne sont pas imbriqués et sont comparés par des tests de rapport de vraisemblance. En bref, voici un exemple reproductible de ce …
J'ai un ensemble de données de mesures répétées déséquilibrées à analyser, et j'ai lu que la façon dont la plupart des progiciels statistiques traitent cela avec l'ANOVA (c'est-à-dire la somme des carrés de type III) est erronée. Par conséquent, je voudrais utiliser un modèle à effets mixtes pour analyser ces …
En spécifiant un modèle à effets mixtes croisés, j'essaie d'inclure des interactions. Cependant, j'obtiens le message d'erreur suivant: Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) Le modèle présente …
J'ai une question sur mon utilisation d'un modèle mixte / lmer. Le modèle de base est le suivant: lmer(DV ~ group * condition + (1|pptid), data= df) Le groupe et la condition sont deux facteurs: le groupe a deux niveaux (groupe A, groupe B) et la condition a trois niveaux …
J'ai réalisé une expérience où j'ai élevé différentes familles provenant de deux populations sources différentes. Chaque famille a reçu l'un des deux traitements. Après l'expérience, j'ai mesuré plusieurs traits sur chaque individu. Pour tester un effet du traitement ou de la source ainsi que leur interaction, j'ai utilisé un modèle …
Supposons que j'ai NNN participants, chacun donnant une réponse YYY 20 fois, 10 dans une condition et 10 dans une autre. J'adapte un modèle linéaire à effets mixtes comparant YYY dans chaque condition. Voici un exemple reproductible simulant cette situation en utilisant le lme4package dans R: library(lme4) fml <- "~ …
Supposons que j'ai une variable de réponse yijyijy_{ij} qui a été mesurée à partir du jjj e frère de la iii e famille. De plus, certaines données comportementales xijxijx_{ij} ont été collectées en même temps auprès de chaque sujet. J'essaie d'analyser la situation avec le modèle linéaire à effets mixtes …
Quelqu'un peut-il me dire comment R peut estimer le point de rupture dans un modèle linéaire par morceaux (en tant que paramètre fixe ou aléatoire), alors que j'ai également besoin d'estimer d'autres effets aléatoires? J'ai inclus un exemple de jouet ci-dessous qui correspond à une régression de bâton de hockey …
Douglas Bates déclare que les modèles suivants sont équivalents «si la matrice de variance-covariance pour les effets aléatoires à valeur vectorielle a une forme spéciale, appelée symétrie composée» ( diapositive 91 dans cette présentation ): m1 <- lmer(y ~ factor + (0 + factor|group), data) m2 <- lmer(y ~ factor …
J'utilise un modèle de régression logistique sous la forme: lmer(response~1+(1|site), family=binomial, REML = FALSE) Normalement, je calculerais l'ICC à partir des variances d'interception et résiduelles, mais le résumé du modèle n'inclut pas la variance résiduelle. Comment puis-je calculer cela?
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.