Puis-je obtenir de l'aide pour terminer cette tentative (en cours) de tentative de repérage sur les équivalents ANOVA et REGRESSION? J'ai essayé de concilier les concepts, la nomenclature et la syntaxe de ces deux méthodologies. Il existe de nombreux messages sur ce site concernant leur similitude,...
J'ai des données de patients traités avec 2 types de traitements différents pendant la chirurgie. J'ai besoin d'analyser son effet sur la fréquence cardiaque. La mesure de la fréquence cardiaque est prise toutes les 15 minutes. Étant donné que la durée de la chirurgie peut être différente pour...
J'ai effectué une évaluation informatisée des différentes méthodes d'ajustement d'un type particulier de modèle utilisé dans les sciences paléo. J'avais un ensemble d'entraînement de grande taille et j'ai donc mis au hasard (échantillonnage aléatoire stratifié) un ensemble de tests de côté. J'ai...
Existe-t-il un tel ensemble qui prévoit une estimation du modèle binomial à effets mixtes négatifs gonflés zéro dans R? J'entends par là: Zero-inflation où vous pouvez spécifier le modèle binomial pour zéro inflation, comme dans la fonction zeroinfl dans le package pscl: zeroinfl (y ~ X | Z, dist =...
Dites que vous êtes dans la bibliothèque de votre département de statistique et que vous tombez sur un livre avec l'image suivante en première page. Vous penserez probablement qu'il s'agit d'un livre sur les choses de régression linéaire. Quelle serait l'image qui vous ferait penser aux modèles...
Disons que nous avons un modèle mod <- Y ~ X*Condition + (X*Condition|subject) # Y = logit variable # X = continuous variable # Condition = values A and B, dummy coded; the design is repeated # so all participants go through both Conditions # subject = random effects for different subjects...
J'ai du mal à comprendre la sortie de mon lmer()modèle. Il s'agit d'un modèle simple d'une variable de résultat (prise en charge) avec différentes interceptions d'état / effets aléatoires d'état: mlm1 <- lmer(Support ~ (1 | State)) Les résultats de summary(mlm1)sont: Linear mixed model fit by...
J'ai une question concernant la technique d'amorçage appropriée à utiliser avec des données où un clustering puissant est présent. J'ai été chargé d'évaluer un modèle prédictif multivarié à effets mixtes sur les données de sinistres d'assurance en notant le modèle de référence actuel sur les...
Étant donné une variable prédite (P), un effet aléatoire (R) et un effet fixe (F), on pourrait adapter deux * modèles à effets mixtes ( syntaxe lme4 ): m1 = lmer( P ~ (1|R) + F ) m2 = lmer( P ~ (1+F|R) + F) Si je comprends bien, le deuxième modèle est celui qui permet à l'effet fixe de varier selon...
Formules de modèle en R telles que y ~ x + a*b + c:d sont basées sur la notation dite de Wilkinson : Wilkinson et Rogers 1973, Description symbolique des modèles factoriels pour l'analyse de la variance . Cet article n'a pas discuté des notations pour les modèles mixtes (qui n'existaient peut-être...
Bonjour gourous statistiques et assistants de programmation R, Je m'intéresse à la modélisation des captures d'animaux en fonction des conditions environnementales et du jour de l'année. Dans le cadre d'une autre étude, j'ai dénombré des captures sur environ 160 jours sur trois ans. Chaque jour,...
Diverses descriptions sur la sélection des modèles sur les effets aléatoires des modèles mixtes linéaires indiquent l'utilisation de REML. Je connais la différence entre REML et ML à un certain niveau, mais je ne comprends pas pourquoi REML devrait être utilisé parce que ML est biaisé. Par exemple,...
J'ai un ensemble de données composé d'une série de décomptes mensuels de «bâtons cassés» provenant d'une poignée de sites. J'essaie d'obtenir une seule estimation récapitulative à partir de deux techniques différentes: Technique 1: Ajustez un «bâton cassé» avec un GLM de Poisson avec une variable...
J'ai utilisé le MCMCglmmpackage récemment. Je suis confus par ce que l'on appelle dans la documentation structure R et structure G. Ceux-ci semblent se rapporter aux effets aléatoires - en particulier en spécifiant les paramètres de la distribution antérieure sur eux, mais la discussion dans la...
J'ai une expérience que je vais essayer de résumer ici. Imaginez que je lance trois pierres blanches devant vous et vous demande de juger de leur position. J'enregistre une variété de propriétés des pierres et votre réponse. Je le fais sur un certain nombre de sujets. Je génère deux modèles. L'un...
Dans un modèle à effets mixtes, la recommandation consiste à utiliser un effet fixe pour estimer un paramètre si tous les niveaux possibles sont inclus (par exemple, hommes et femmes). Il est en outre recommandé d'utiliser un effet aléatoire pour tenir compte d'une variable si les niveaux inclus ne...
Je connais assez bien les modèles à effets mixtes (MEM), mais un collègue m'a récemment demandé comment ils se comparent aux modèles de croissance latente (LGM). J'ai fait un peu de recherche sur Google, et il semble que LGM soit une variante de la modélisation des équations structurelles qui est...
L'un des problèmes que j'ai toujours eu avec les modèles mixtes est de trouver des visualisations de données - du type qui pourraient se retrouver sur un papier ou une affiche - une fois que l'on a les résultats. En ce moment, je travaille sur un modèle d'effets mixtes de Poisson avec une formule...
J'essaie d'utiliser lme4::glmer()pour adapter un modèle mixte généralisé binomial (GLMM) avec une variable dépendante qui n'est pas binaire, mais une variable continue entre zéro et un. On peut considérer cette variable comme une probabilité; en fait , il est probable que rapporté par des sujets...
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...