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...
Pour un exemple simple, supposons qu'il existe deux modèles de régression linéaire Modèle 1 a trois prédicteurs, x1a, x2betx2c Le modèle 2 a trois prédicteurs du modèle 1 et deux prédicteurs supplémentaires x2aetx2b Il existe une équation de régression de la population où la variance de la...
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'essaie de calculer les prédictions d'effet aléatoire à partir d'un modèle mixte linéaire à la main, et en utilisant la notation fournie par Wood dans les modèles additifs généralisés: une introduction avec R (pg 294 / pg 307 du pdf), je suis confus sur ce que chaque paramètre représente. Voici un...
La formule que l'on doit spécifier pour la formation d'un modèle à plusieurs niveaux (à l'aide lmerde la lme4 Rbibliothèque) m'obtient toujours. J'ai lu d'innombrables manuels et tutoriels, mais je ne l'ai jamais bien compris. Voici donc un exemple de ce tutoriel que j'aimerais voir formulé dans...
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'essaie de reproduire plusieurs tests d'interaction entre les deux lmet lmersur des mesures répétées (2x2x2). La raison pour laquelle je veux comparer les deux méthodes est que le GLM de SPSS pour les mesures répétées donne exactement les mêmes résultats que l' lmapproche présentée ici, donc à la...
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...
Comment puis-je obtenir des effets aléatoires groupés pour lmer après une imputation multiple? J'utilise des souris pour imputer plusieurs trames de données. Et lme4 pour un modèle mixte avec interception aléatoire et pente aléatoire. La mise en commun de lmer se passe bien, sauf qu'elle ne...
É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é....
J'effectue des tests post-hoc sur un modèle linéaire à effets mixtes dans R( lme4package). J'utilise multcomppackage ( glht()fonction) pour effectuer les tests post-hoc. Mon plan expérimental est des mesures répétées, avec un effet de bloc aléatoire. Les modèles sont spécifiés comme: mymod <-...
J'analyse les données oculométriques d'une expérience conçue. Une version simplifiée de mes données ressemble à ceci (vous pouvez obtenir les données dput () ici ), head(lookDATA) participant fixationImage fixationCount 1 9 Automobile 81 2 9 Bird 63 3 9 Chair 82 4 9 Dog 64 5 9 Face 90 6 9 Plant 75...
Dans la lmerfonction au sein lme4de Ril y a un appel à la construction d' une matrice de modèle d'effets aléatoires, , comme expliqué ici , pages 7 - 9.ZZZ Le calcul de implique les produits de KhatriRao et / ou Kronecker de deux matrices, et . J i X iZZZJiJiJ_iXiXiX_i La matrice est une bouchée:...
En contrepartie de ce poste , j'ai travaillé sur la simulation de données avec des variables continues, se prêtant à des interceptions et des pentes corrélées. Bien qu'il existe d'excellents articles sur ce sujet sur le site et en dehors du site , j'ai eu du mal à trouver un exemple de bout en bout...
Donc, j'ai 16 essais dans lesquels j'essaie d'authentifier une personne à partir d'un trait biométrique en utilisant Hamming Distance. Mon seuil est fixé à 3,5. Mes données sont ci-dessous et seul l'essai 1 est un vrai positif: Trial Hamming Distance 1 0.34 2 0.37 3 0.34 4 0.29 5 0.55 6 0.47 7 0.47...
Je monte un glmermodèle dans le lme4package R. Je cherche une table anova avec une valeur de p montrée ici, mais je ne trouve aucun paquet qui lui convient. Est-il possible de le faire en R? Le modèle que je monte est de la forme: model1<-glmer(dmn~period*teethTreated+(1|fullName),...
J'ai du mal à obtenir des résultats équivalents entre un aovmodèle de mesures répétées inter-intra et un lmermodèle mixte. Mes données et mon script se présentent comme suit data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST...
Dans un modèle mixte linéaire généralisé logistique (famille = binôme), je ne sais pas comment interpréter la variance à effets aléatoires: Random effects: Groups Name Variance Std.Dev. HOSPITAL (Intercept) 0.4295 0.6554 Number of obs: 2275, groups: HOSPITAL, 14 Comment interpréter ce résultat...
Je simule des essais de Bernoulli avec un aléatoire entre les groupes, puis j'adapte le modèle correspondant avec le paquet:logitθ ∼ N( logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number of...