J'ai cherché cela en ligne pendant des heures, mais aucun des messages en ligne n'est ce que je recherche. Ma question est très facile à implémenter dans la procédure mixte SAS Proc mais je ne sais pas comment le faire dans les packages lme et / ou lmer. Supposons que j'ai un modèle, , où est fixe mais et sont aléatoires. Mon code R est
f1 = lme(y ~ factor(a), data = mydata,
random = list(factor(b) = ~ 1, factor(a):factor(b) = ~ 1))
Erreur: inattendue =
dans:
f1 = lme(y ~ factor(a), data = mydata,
random = list(factor(a) =
Quelqu'un pourrait-il me dire comment spécifier ces effets aléatoires dans lme? Merci d'avance
dput
obtenir le code nécessaire pour recréer vos données. D'après le commentaire que vous avez laissé, le résultat eststructure(list(method = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"), day = structure(c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L, 1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L), .Label = c("1", "2", "3", "4"), class = "factor"), level = c(142.3, 144, 134.9, 146.3, 148.6, 156.5, 152, 151.4, 142.9, 147.4, 125.9, 127.6, 135.5, 138.9, 142.9, 142.3)), .Names = c("method", "day", "level"), row.names = c(NA, -16L), class = "data.frame")
Réponses:
Essayez ceci, c'est une façon standard de faire un tracé divisé. La notation
/
signifie que la méthode est imbriquée dans day.la source
Cela aiderait beaucoup si vous fournissiez un data.frame. Maintenant, il n'est pas clair ce qu'est un facteur de regroupement. Je juge que c'est . Ensuite, en notation, votre modèle doit être écrit comme suit:β
lme
la source