Displaying 3 results from an estimated 3 matches for "basehr".
Did you mean:
baehr
2006 May 15
1
anova statistics in lmer
...ws 5: 27-30)
> data("HR",package="SASmixed")
> library(lme4)
Loading required package: Matrix
Loading required package: lattice
Attaching package: 'lattice'
The following object(s) are masked from package:Matrix :
qqmath
> (fm1<-lmer(HR~baseHR+Time*Drug+(1|Patient),HR))
Linear mixed-effects model fit by REML
Formula: HR ~ baseHR + Time * Drug + (1 | Patient)
Data: HR
AIC BIC logLik MLdeviance REMLdeviance
788.6769 810.9768 -386.3384 791.8952 772.6769
Random effects:
Groups Name Variance Std.Dev....
2012 Jul 20
1
Extracting standard errors for adjusted fixed effect sizes in lmer
...ckage = "SASmixed")
> str(HR)
'data.frame': 120 obs. of 5 variables:
$ Patient: Factor w/ 24 levels "201","202","203",..: 1 1 1 1 1 2 2 2 2 2
...
$ Drug : Factor w/ 3 levels "a","b","p": 3 3 3 3 3 2 2 2 2 2 ...
$ baseHR : num 92 92 92 92 92 54 54 54 54 54 ...
$ HR : num 76 84 88 96 84 58 60 60 60 64 ...
$ Time : num 0.0167 0.0833 0.25 0.5 1 ...
> fm1 <- lmer(HR ~ baseHR + Time + Drug + (1 | Patient), HR)
> fixef(fm1) ##Extract estimates of fixed effects
(Intercept) baseHR Time...
2003 May 02
2
stepAIC/lme (1.6.2)
Based on the stepAIC help, I have assumed that it only was for lm, aov, and
glm models. I gather from the following correspondence that it also works
with lme models.
Thomas Lumley 07:40 a.m. 28/04/03 -0700 4 Re: [R] stepAIC/lme problem
(1.7.0 only)
Prof Brian Ripley 04:19 p.m. 28/04/03 +0100 6 Re: [R] stepAIC/lme problem
(1.7.0 only)
Prof Brian Ripley 06:09 p.m. 29/04/03 +0100 6 Re: [R]