R: diagramme des Interactions avec une constante et une variable catégorique pour un GLMM (lme4)
Je voudrais faire un diagramme des interactions pour représenter visuellement la différence ou la similitude dans les pentes de l'interaction d'une variable catégorielle (4 niveaux) et une assiette variable en continu à partir des résultats d'un modèle de régression.
with(GLMModel, interaction.plot(continuous.var, categorical.var, response.var))
N'est pas ce que je recherche. Elle produit d'un complot dans lequel les changements de pente pour chaque valeur de la variable continue. Je suis à la recherche de faire une intrigue avec des constantes et des pistes dans la suite de l'intrigue:
Des idées?
Je adapter à un modèle de la forme fit<-glmer(resp.var ~ cont.var*cat.var + (1|rand.eff) , data = sample.data , poisson)
Voici quelques exemples de données:
structure(list(cat.var = structure(c(4L, 4L, 1L, 4L, 1L, 2L,
1L, 1L, 1L, 1L, 4L, 1L, 1L, 3L, 2L, 4L, 1L, 1L, 1L, 2L, 1L, 2L,
2L, 1L, 3L, 1L, 1L, 2L, 4L, 1L, 2L, 1L, 1L, 4L, 1L, 3L, 1L, 3L,
3L, 4L, 3L, 4L, 1L, 3L, 3L, 1L, 2L, 3L, 4L, 3L, 4L, 2L, 1L, 1L,
4L, 1L, 1L, 1L, 1L, 1L, 1L, 4L, 1L, 4L, 4L, 3L, 3L, 1L, 3L, 3L,
3L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 4L, 1L, 3L, 4L, 1L, 1L, 4L,
1L, 3L, 1L, 1L, 3L, 2L, 4L, 1L, 4L, 1L, 4L, 4L, 4L, 4L, 2L, 4L,
4L, 1L, 2L, 1L, 4L, 3L, 1L, 1L, 3L, 2L, 4L, 4L, 1L, 4L, 1L, 3L,
2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 4L, 1L,
2L, 2L, 1L, 1L, 2L, 3L, 1L, 4L, 4L, 4L, 1L, 4L, 4L, 3L, 2L, 4L,
1L, 3L, 1L, 1L, 4L, 4L, 2L, 4L, 1L, 1L, 3L, 4L, 2L, 1L, 3L, 3L,
4L, 3L, 2L, 3L, 1L, 4L, 2L, 2L, 1L, 4L, 1L, 2L, 3L, 4L, 1L, 4L,
2L, 1L, 3L, 3L, 3L, 4L, 1L, 1L, 1L, 3L, 1L, 3L, 4L, 2L, 1L, 4L,
1L, 1L, 1L, 2L, 1L, 1L, 4L, 1L, 3L, 1L, 2L, 1L, 4L, 1L, 2L, 4L,
1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 3L, 1L, 3L, 4L, 1L, 4L, 3L,
3L, 3L, 4L, 1L, 3L, 1L, 1L, 4L, 4L, 4L, 4L, 2L, 1L, 1L, 3L, 2L,
1L, 4L, 4L, 2L, 4L, 2L, 4L, 1L, 3L, 4L, 1L, 1L, 2L, 3L, 2L, 4L,
1L, 1L, 3L, 4L, 2L, 2L, 3L, 4L, 1L, 2L, 3L, 1L, 2L, 4L, 1L, 4L,
2L, 4L, 3L, 4L, 2L, 1L, 1L, 1L, 1L, 1L, 4L, 4L, 1L, 4L, 4L, 1L,
4L, 2L, 1L, 1L, 1L, 1L, 3L, 1L, 1L, 3L, 3L, 2L, 2L, 1L, 1L, 4L,
1L, 4L, 3L, 1L, 2L, 1L, 4L, 2L, 4L, 4L, 1L, 2L, 1L, 1L, 1L, 4L,
1L, 4L, 1L, 2L, 1L, 3L, 1L, 3L, 3L, 1L, 1L, 4L, 3L, 1L, 4L, 1L,
2L, 4L, 1L, 1L, 3L, 3L, 2L, 4L, 4L, 1L, 1L, 2L, 2L, 1L, 2L, 4L,
3L, 4L, 4L, 4L, 4L, 1L, 3L, 1L, 2L, 2L, 2L, 4L, 2L, 3L, 4L, 1L,
3L, 2L, 2L, 1L, 1L, 1L, 3L, 1L, 2L, 2L, 1L, 1L, 3L, 2L, 1L, 1L,
1L, 1L, 2L, 1L, 1L, 1L, 4L, 4L, 4L, 3L, 3L, 2L, 1L, 3L, 2L, 1L,
1L, 1L, 4L, 1L, 1L, 2L, 3L, 1L, 1L, 2L, 4L, 3L, 2L, 4L, 3L, 2L,
1L, 3L, 1L, 3L, 1L, 4L, 3L, 1L, 4L, 4L, 2L, 4L, 1L, 1L, 2L, 4L,
4L, 2L, 3L, 4L, 4L, 3L, 1L, 4L, 1L, 2L, 4L, 1L, 1L, 4L, 1L, 1L,
1L, 1L, 1L, 3L, 4L, 1L, 4L, 4L, 2L, 2L, 2L, 2L, 3L, 4L, 4L, 1L,
1L, 4L, 2L, 3L, 3L, 1L, 1L, 1L, 1L, 3L, 1L, 1L, 1L, 3L, 4L, 2L,
3L, 1L, 1L, 1L, 4L, 1L, 1L, 4L, 4L, 4L, 1L, 1L, 1L, 1L), .Label = c("A",
"B", "C", "D"), class = "factor"), cont.var = c(-0.0682900527296927,
0.546320421837542, -0.273160210918771, -0.887770685486005, 0.136580105459385,
0.75119058002662, 0.546320421837542, -0.273160210918771, -0.682900527296927,
0.136580105459385, 0.75119058002662, 0.75119058002662, 0.75119058002662,
0.341450263648464, 0.75119058002662, 0.546320421837542, 0.546320421837542,
-0.478030369107849, -0.478030369107849, -0.682900527296927, -0.682900527296927,
0.546320421837542, -0.478030369107849, -0.0682900527296927, 0.136580105459385,
0.136580105459385, 0.75119058002662, -0.478030369107849, 0.75119058002662,
-0.887770685486005, 0.136580105459385, -0.478030369107849, 0.341450263648464,
-0.682900527296927, -0.478030369107849, 0.341450263648464, -0.478030369107849,
0.546320421837542, 0.75119058002662, -0.478030369107849, -0.273160210918771,
0.546320421837542, -0.682900527296927, 0.75119058002662, -0.478030369107849,
-0.887770685486005, 0.136580105459385, -0.887770685486005, -0.0682900527296927,
-0.478030369107849, 0.546320421837542, 0.75119058002662, 0.136580105459385,
-0.273160210918771, -0.273160210918771, 0.75119058002662, -0.682900527296927,
0.136580105459385, -0.273160210918771, -0.273160210918771, 0.136580105459385,
0.136580105459385, 0.341450263648464, 0.136580105459385, -0.273160210918771,
-0.273160210918771, -0.682900527296927, -0.887770685486005, -0.0682900527296927,
0.136580105459385, -0.0682900527296927, -0.273160210918771, -0.273160210918771,
0.341450263648464, 0.75119058002662, -0.682900527296927, -0.0682900527296927,
-0.273160210918771, -0.887770685486005, -0.0682900527296927,
0.75119058002662, 0.546320421837542, 0.75119058002662, 0.75119058002662,
-0.887770685486005, 0.341450263648464, 0.75119058002662, -0.887770685486005,
0.136580105459385, -0.273160210918771, 0.546320421837542, 0.546320421837542,
-0.682900527296927, 0.75119058002662, 0.136580105459385, -0.0682900527296927,
-0.478030369107849, 0.75119058002662, -0.478030369107849, 0.341450263648464,
0.136580105459385, -0.0682900527296927, -0.478030369107849, -0.0682900527296927,
-0.0682900527296927, 0.546320421837542, -0.273160210918771, 0.75119058002662,
0.341450263648464, 0.546320421837542, -0.478030369107849, 0.136580105459385,
-0.887770685486005, -0.273160210918771, -0.273160210918771, -0.478030369107849,
-0.478030369107849, 0.75119058002662, -0.682900527296927, -0.0682900527296927,
0.546320421837542, 0.75119058002662, 0.546320421837542, 0.136580105459385,
-0.478030369107849, 0.136580105459385, 0.546320421837542, -0.478030369107849,
-0.0682900527296927, -0.0682900527296927, 0.546320421837542,
-0.273160210918771, 0.136580105459385, -0.0682900527296927, 0.75119058002662,
-0.0682900527296927, 0.546320421837542, -0.887770685486005, -0.0682900527296927,
-0.682900527296927, -0.478030369107849, -0.478030369107849, -0.682900527296927,
0.75119058002662, 0.341450263648464, -0.0682900527296927, 0.341450263648464,
-0.0682900527296927, -0.887770685486005, -0.887770685486005,
-0.273160210918771, -0.0682900527296927, 0.546320421837542, -0.0682900527296927,
-0.0682900527296927, 0.75119058002662, -0.0682900527296927, -0.273160210918771,
-0.478030369107849, 0.546320421837542, 0.546320421837542, 0.546320421837542,
0.341450263648464, 0.136580105459385, -0.478030369107849, 0.136580105459385,
0.136580105459385, 0.136580105459385, -0.478030369107849, -0.273160210918771,
-0.273160210918771, -0.273160210918771, 0.341450263648464, -0.273160210918771,
-0.0682900527296927, 0.136580105459385, 0.546320421837542, -0.478030369107849,
-0.273160210918771, 0.546320421837542, 0.546320421837542, -0.273160210918771,
-0.0682900527296927, 0.341450263648464, 0.546320421837542, -0.0682900527296927,
0.136580105459385, -0.478030369107849, 0.75119058002662, -0.478030369107849,
-0.682900527296927, -0.478030369107849, 0.136580105459385, -0.273160210918771,
-0.0682900527296927, -0.887770685486005, -0.887770685486005,
0.546320421837542, -0.273160210918771, 0.546320421837542, -0.478030369107849,
0.546320421837542, -0.0682900527296927, 0.75119058002662, -0.273160210918771,
0.546320421837542, 0.341450263648464, -0.0682900527296927, -0.0682900527296927,
-0.0682900527296927, -0.887770685486005, 0.136580105459385, -0.273160210918771,
-0.478030369107849, 0.75119058002662, 0.341450263648464, 0.546320421837542,
-0.273160210918771, 0.546320421837542, 0.75119058002662, -0.273160210918771,
0.75119058002662, 0.546320421837542, -0.273160210918771, -0.273160210918771,
0.75119058002662, -0.273160210918771, -0.0682900527296927, 0.136580105459385,
-0.478030369107849, 0.75119058002662, 0.75119058002662, -0.887770685486005,
-0.887770685486005, 0.546320421837542, -0.682900527296927, -0.887770685486005,
0.136580105459385, 0.75119058002662, 0.75119058002662, -0.478030369107849,
0.136580105459385, 0.75119058002662, -0.273160210918771, -0.682900527296927,
-0.273160210918771, 0.136580105459385, 0.546320421837542, -0.682900527296927,
-0.478030369107849, 0.136580105459385, -0.682900527296927, -0.0682900527296927,
-0.478030369107849, 0.136580105459385, -0.887770685486005, -0.273160210918771,
-0.0682900527296927, -0.273160210918771, -0.887770685486005,
0.546320421837542, 0.546320421837542, -0.478030369107849, -0.273160210918771,
-0.0682900527296927, 0.136580105459385, -0.478030369107849, 0.75119058002662,
0.341450263648464, 0.136580105459385, 0.136580105459385, 0.75119058002662,
0.136580105459385, -0.0682900527296927, 0.546320421837542, -0.0682900527296927,
-0.887770685486005, 0.75119058002662, 0.75119058002662, 0.546320421837542,
-0.887770685486005, -0.0682900527296927, -0.682900527296927,
-0.682900527296927, 0.75119058002662, 0.75119058002662, -0.478030369107849,
0.546320421837542, -0.273160210918771, 0.75119058002662, -0.0682900527296927,
0.546320421837542, -0.0682900527296927, -0.273160210918771, 0.546320421837542,
0.75119058002662, -0.0682900527296927, 0.546320421837542, -0.682900527296927,
-0.273160210918771, -0.0682900527296927, -0.478030369107849,
-0.478030369107849, 0.136580105459385, -0.273160210918771, 0.136580105459385,
0.546320421837542, 0.75119058002662, -0.273160210918771, 0.341450263648464,
-0.273160210918771, 0.136580105459385, 0.546320421837542, 0.546320421837542,
0.136580105459385, 0.136580105459385, -0.682900527296927, 0.341450263648464,
0.341450263648464, -0.273160210918771, -0.682900527296927, -0.0682900527296927,
0.75119058002662, -0.887770685486005, -0.478030369107849, -0.273160210918771,
-0.478030369107849, -0.478030369107849, 0.136580105459385, -0.478030369107849,
0.136580105459385, -0.478030369107849, 0.136580105459385, -0.0682900527296927,
-0.273160210918771, 0.136580105459385, 0.341450263648464, -0.478030369107849,
0.75119058002662, 0.136580105459385, 0.341450263648464, 0.546320421837542,
-0.887770685486005, 0.75119058002662, 0.341450263648464, -0.0682900527296927,
-0.478030369107849, 0.546320421837542, 0.136580105459385, -0.682900527296927,
-0.0682900527296927, 0.341450263648464, -0.478030369107849, -0.0682900527296927,
-0.478030369107849, -0.0682900527296927, 0.341450263648464, -0.478030369107849,
-0.682900527296927, 0.75119058002662, -0.478030369107849, -0.682900527296927,
0.341450263648464, -0.887770685486005, -0.478030369107849, 0.546320421837542,
-0.887770685486005, -0.478030369107849, -0.478030369107849, 0.341450263648464,
0.75119058002662, -0.682900527296927, 0.75119058002662, 0.75119058002662,
0.341450263648464, -0.0682900527296927, 0.546320421837542, -0.0682900527296927,
0.136580105459385, 0.136580105459385, 0.136580105459385, 0.136580105459385,
0.546320421837542, 0.546320421837542, -0.0682900527296927, 0.75119058002662,
-0.0682900527296927, -0.0682900527296927, -0.682900527296927,
-0.273160210918771, -0.682900527296927, -0.478030369107849, 0.136580105459385,
0.75119058002662, 0.546320421837542, 0.341450263648464, -0.887770685486005,
-0.0682900527296927, 0.136580105459385, 0.75119058002662, -0.273160210918771,
-0.682900527296927, 0.136580105459385, -0.478030369107849, -0.273160210918771,
-0.273160210918771, 0.136580105459385, 0.341450263648464, -0.478030369107849,
-0.0682900527296927, -0.682900527296927, 0.75119058002662, -0.273160210918771,
-0.478030369107849, -0.0682900527296927, -0.0682900527296927,
-0.273160210918771, -0.0682900527296927, -0.478030369107849,
0.75119058002662, -0.0682900527296927, 0.136580105459385, 0.546320421837542,
0.546320421837542, -0.478030369107849, -0.273160210918771, 0.546320421837542,
-0.478030369107849, -0.682900527296927, 0.75119058002662, -0.0682900527296927,
-0.682900527296927, -0.682900527296927, 0.75119058002662, 0.341450263648464,
-0.478030369107849, 0.75119058002662, 0.136580105459385, -0.887770685486005,
0.341450263648464, 0.341450263648464, 0.546320421837542, -0.273160210918771,
0.136580105459385, 0.75119058002662, -0.0682900527296927, -0.682900527296927,
-0.478030369107849, -0.478030369107849, 0.75119058002662, 0.546320421837542,
-0.478030369107849, 0.546320421837542, 0.136580105459385, -0.887770685486005,
0.75119058002662, -0.0682900527296927, 0.75119058002662, 0.75119058002662,
-0.273160210918771, -0.682900527296927, 0.546320421837542, 0.546320421837542,
-0.887770685486005, 0.75119058002662, -0.273160210918771, 0.546320421837542,
-0.0682900527296927, 0.136580105459385, 0.341450263648464, -0.478030369107849,
0.136580105459385, 0.136580105459385, -0.273160210918771, 0.546320421837542,
-0.273160210918771, -0.273160210918771, -0.273160210918771, 0.75119058002662,
-0.887770685486005, -0.887770685486005, -0.0682900527296927,
-0.478030369107849, -0.0682900527296927, 0.75119058002662, -0.273160210918771,
0.136580105459385, -0.478030369107849, -0.273160210918771, 0.136580105459385,
0.75119058002662, 0.546320421837542, -0.478030369107849, -0.273160210918771,
-0.273160210918771, 0.136580105459385, -0.273160210918771, -0.0682900527296927,
0.75119058002662, 0.136580105459385), resp.var = c(2L, 1L, 0L,
1L, 0L, 0L, 0L, 0L, 0L, 1L, 3L, 1L, 0L, 1L, 0L, 1L, 2L, 0L, 1L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 2L,
1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 2L,
0L, 3L, 2L, 0L, 2L, 2L, 0L, 0L, 0L, 1L, 1L, 3L, 1L, 2L, 0L, 1L,
0L, 0L, 1L, 0L, 2L, 0L, 2L, 4L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 2L,
3L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 2L,
0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 2L, 0L, 1L, 0L, 4L, 1L, 0L,
1L, 1L, 0L, 0L, 0L, 1L, 3L, 0L, 2L, 0L, 0L, 2L, 1L, 0L, 0L, 2L,
0L, 0L, 0L, 2L, 0L, 0L, 3L, 0L, 0L, 2L, 1L, 1L, 0L, 0L, 3L, 1L,
1L, 2L, 0L, 2L, 0L, 2L, 2L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 2L, 2L, 1L, 0L, 0L, 1L,
0L, 0L, 0L, 0L, 6L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 0L,
1L, 0L, 0L, 1L, 3L, 1L, 0L, 2L, 3L, 0L, 0L, 1L, 0L, 0L, 1L, 1L,
0L, 0L, 0L, 0L, 1L, 2L, 1L, 1L, 0L, 0L, 2L, 0L, 2L, 0L, 0L, 1L,
1L, 0L, 0L, 2L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L,
0L, 1L, 0L, 2L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L,
0L, 3L, 0L, 0L, 3L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L, 2L, 1L, 1L, 0L, 2L, 2L, 0L, 2L, 1L, 0L, 2L, 0L, 0L, 0L, 0L,
3L, 0L, 2L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 2L, 0L, 1L, 1L, 0L, 1L,
0L, 3L, 1L, 3L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 2L, 0L,
2L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, 2L, 0L, 3L, 0L, 0L, 0L,
0L, 1L, 0L, 0L, 3L, 1L, 1L, 2L, 0L, 0L, 3L, 0L, 0L, 0L, 1L, 1L,
0L, 1L, 3L, 0L, 2L, 0L, 0L, 1L, 3L, 1L, 0L, 0L, 4L, 3L, 0L, 2L,
0L, 0L, 0L, 3L, 0L, 0L, 2L, 3L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L,
0L, 0L, 0L, 3L, 3L, 2L, 0L, 0L, 2L, 0L, 0L, 0L, 0L, 2L, 0L, 0L,
0L, 0L, 0L, 1L, 0L, 2L, 0L, 0L, 1L, 0L, 0L, 1L, 2L, 0L, 1L, 0L,
2L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 3L, 1L, 0L, 0L, 0L, 0L, 0L,
1L, 2L, 0L, 2L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L,
0L, 0L, 3L, 2L, 2L, 0L, 1L, 0L, 5L, 0L, 4L, 2L, 0L, 3L, 0L, 0L,
1L, 1L, 0L, 0L, 0L, 2L, 0L, 1L, 0L, 3L, 0L, 2L, 0L, 0L, 0L, 2L,
0L), rand.eff = c(37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L,
37L, 37L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L, 40L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L,
43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L, 43L)), .Names = c("cat.var",
"cont.var", "resp.var", "rand.eff"), row.names = c(NA, 500L), class = "data.frame")
- pas trop dur en combinant
predict
avecggplot
oulattice::xyplot
. Reproductible exemple s'il vous plaît?
Vous devez vous connecter pour publier un commentaire.
Voici une réponse de toutes sortes (en passant, vous avez eu, il manque des guillemets dans votre bloc de données ci-dessus, qui ont du être corrigés manuellement ...)
L'ajustement du modèle:
(Notez que c'est un peu bizarre de spécification de modèle -- les forces de toutes les catégories ont la même valeur en
cont.var==0
. Vouliez-vous direcont.var*cat.var
?Rapide et sale régressions linéaires:
Maintenant avec une loi de Poisson GLM (mais ne tenant pas compte de l'effet aléatoire), et en montrant les points de données:
Le prochain bit exige le développement (r-forge) version de
lme4
, qui a unpredict
méthode:Configurer le bloc de données à des fins de prédiction:
Prévoir au niveau de la population (
REform=NA
), sur le prédicteur linéaire (logit) échelle (c'est la seule façon, vous obtiendrez des lignes droites sur la parcelle)Maintenant, sur la réponse de l'échelle:
*
au lieu de:
dans la spécification du modèle.install.packages("lme4",repos="http://lme4.r-forge.r-project.org/repos")
-- le construire sur le principal r-forge référentiel est brisé à l'instant (vous pouvez toujours re-installer à partir de CRAN si nécessaire).La
jtools
paquet (CRAN lien) peut faire le tracé de ce genre de modèle assez simple. Je suis le développeur de ce package.Nous l'ajustement du modèle comme Ben l'a fait dans sa réponse:
Et avec
jtools
nous suffit d'utiliser lainteract_plot
fonction comme ceci:Le résultat:
Par défaut des parcelles de la réponse de l'échelle, mais vous pouvez l'avoir tracée sur l'échelle linéaire avec la
outcome.scale = "link"
argument (la valeur par défaut est"response"
).La effets package de support pour les
lme4
modèles, et devrait être en mesure de faire ce que vous voulez.Il est également livré avec deux légèrement dépassées papiers (vous pouvez y penser comme des vignettes).