[R] "effects" package with "lme4"

[Axel Urbiz]

Dear John,

Thank you for your response. 

My apologies as I’m only recently getting exposed to mixed-model. My understanding, is that the model specified below also has random intercepts and slope, as they vary by Subject. `coef(fm1)` shows this. I was looking to plot the fitted splines by Subject.

Sorry if my interpretation is incorrect. 

Best,
Axel.



> On May 15, 2020, at 5:51 PM, Fox, John <jfox using mcmaster.ca> wrote:
> 
> Dear Axel,
> 
> There only one fixed effect in the model, ns(Days, 3), so I don't know what you expected.
> 
> Best,
> John
> --------------------------------------
> John Fox, Professor Emeritus
> McMaster University
> Hamilton, Ontario, Canada
> Web: socialsciences.mcmaster.ca/jfox/
> 
> 
> 
>> -----Original Message-----
>> From: Axel Urbiz <axel.urbiz using gmail.com>
>> Sent: Friday, May 15, 2020 5:33 PM
>> To: Fox, John <jfox using mcmaster.ca>; R-help using r-project.org
>> Subject: "effects" package with "lme4"
>> 
>> Hello John and others,
>> 
>> I’d appreciate your help as I’m trying to plot the effect of predictor
>> “Days” on Reaction by Subject. I’m only getting one plot in the example
>> below.
>> 
>> ### Start example
>> 
>> library(lme4)
>> library(splines)
>> data("sleepstudy")
>> 
>> fm1 <- lmer(Reaction ~ ns(Days, 3) + (ns(Days, 3) | Subject), sleepstudy)
>> coef(fm1)
>> plot(allEffects(fm1))
>> 
>> ### End example
>> 
>> Thanks,
>> Axel.