J'ai parcouru Google et ce site et je suis toujours confus à propos de la fonction lmer dans la bibliothèque lme4. J'ai quelques données recueillies auprès de différents services psychiatriques, qui ont une structure à plusieurs niveaux. Pour simplifier, je vais choisir deux variables de niveau 2...
L'occurrence pas si rare lorsqu'il s'agit de modèles mixtes maximaux complexes (estimation de tous les effets aléatoires possibles pour des données et un modèle donnés) est parfaite (+1 ou -1) ou une corrélation presque parfaite entre certains effets aléatoires. Aux fins de la discussion, observons...
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...
Supposons que, dans une étude de 15 sujets, la variable de réponse (res) est modélisée avec deux variables explicatives, une (niveau) est catégorique avec 5 niveaux et l'autre (temps de réponse: RT) est continue. Avec lmer dans le package lme4 de R, j'ai: fm1 <- lmer(res ~ level * RT + (level-1...
Cette question peut sembler stupide, mais ... est-il exact que les effets aléatoires ne peuvent s'appliquer qu'aux variables catégorielles (comme l'identifiant individuel, l'identifiant de la population, ...), par exemple, disons que est une variable catégorielle:xixix_i yiyiy_i ~βxiβxi\beta_{x_i}...
J'ai le type de données suivant. J'ai évalué 10 individus chacun répété 10 fois. J'ai une matrice de relations 10x10 (relation entre toutes les combinaisons des individus). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5,...
Mon collègue et moi ajustons une gamme de modèles d'effets mixtes linéaires et non linéaires dans R. On nous demande d'effectuer une validation croisée sur les modèles ajustés afin que l'on puisse vérifier que les effets observés sont relativement généralisables. C'est normalement une tâche...
Dans R, comment spécifier le modèle lmer sans effet fixe global? Par exemple, si je dis quelque chose comme lmer(y ~ (1 | group) + (0 + x | group), data = my_df) le modèle équipé sera yje j= a +αje+βjeXje jyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Comment puis-je adapter le modèle...
J'essaie de faire face à une analyse du temps sur événement en utilisant des résultats binaires répétés. Supposons que le temps écoulé soit mesuré en jours, mais pour le moment, nous discrétisons le temps en semaines. Je veux approximer un estimateur de Kaplan-Meier (mais tenir compte des...
Les greffes suivantes sont extraites de cet article . Je suis novice dans le bootstrap et j'essaie d'implémenter le bootstrap paramétrique, semi-paramétrique et non paramétrique pour le modèle mixte linéaire avec le R bootpackage. Code R Voici mon Rcode: library(SASmixed) library(lme4)...
J'essaie d'analyser l'effet de l'année sur logInd variable pour un groupe particulier d'individus (j'ai 3 groupes). Le modèle le plus simple: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data = mydata)...
J'essaie de produire un modèle mixte linéaire, le code R est le suivant. lme (Average.payoff ~ Game + Type + Others.Type + Game: Type + Game: Others.Type + Type: Others.Type, random = ~ 1 | Subjects, method = "REML", data = Subjectsm1) -> lme1 Le terme de réponse Average.payoff est continu...
Dans une étude longitudinale, les résultats des unités sont mesurés à plusieurs reprises aux points temporels avec un total de occasions de mesure fixes (fixe = les mesures sur les unités sont prises en même temps). i t mYitYitY_{it}iiitttmmm Les unités sont assignées au hasard soit à un...
J'ai une question de base concernant les approches de la moyenne des modèles à l'aide de critères informatiques pour pondérer les modèles au sein d'un ensemble candidat. La plupart des sources que j'ai lues sur la moyenne du modèle préconisent de faire la moyenne des estimations des coefficients...
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),...
Dois-je inclure des effets aléatoires dans un modèle même s'ils ne sont pas statistiquement significatifs? J'ai un plan expérimental de mesures répétées, dans lequel chaque individu expérimente trois traitements différents dans un ordre aléatoire. Je voudrais contrôler les effets de l'individu et...
Exemples: J'ai une phrase dans la description de poste: "Java senior engineer in UK". Je veux utiliser un modèle d'apprentissage profond pour le prédire en 2 catégories: English et IT jobs. Si j'utilise un modèle de classification traditionnel, il ne peut prédire qu'une seule étiquette avec...
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...
La distinction entre les variables strictement aléatoires (qui devraient être modélisées en tant que telles) et les variables non aléatoires qui, selon certains, pourraient être modélisées comme aléatoires s'il s'agit d'un modèle hiérarchique / à plusieurs niveaux, est floue pour moi. Bates et...
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...