The R Project SVN R-packages

Rev

Rev 4708 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed 


R version 2.4.0 beta (2006-09-25 r39510)
Copyright (C) 2006 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(mlmRev)
Loading required package: lme4
Loading required package: Matrix
Loading required package: lattice
> options(show.signif.stars = FALSE)
> cntr <- list()
> (fm01 <- lmer(langPOST ~ IQ.ver.cen+avg.IQ.ver.cen+(1|schoolNR), bdf))
Linear mixed-effects model fit by REML 
Formula: langPOST ~ IQ.ver.cen + avg.IQ.ver.cen + (1 | schoolNR) 
   Data: bdf 
   AIC   BIC logLik MLdeviance REMLdeviance
 15240 15263  -7616      15228        15232
Random effects:
 Groups   Name        Variance Std.Dev.
 schoolNR (Intercept)  7.8859  2.8082  
 Residual             42.1723  6.4940  
number of obs: 2287, groups: schoolNR, 131

Fixed effects:
               Estimate Std. Error t value
(Intercept)    40.74101    0.28660  142.16
IQ.ver.cen      2.41477    0.07168   33.69
avg.IQ.ver.cen  1.58916    0.31477    5.05

Correlation of Fixed Effects:
            (Intr) IQ.vr.
IQ.ver.cen   0.000       
avg.IQ.vr.c  0.077 -0.228
> (fm02 <- lmer(langPOST ~ IQ.ver.cen+avg.IQ.ver.cen+(IQ.ver.cen|schoolNR), bdf))
Linear mixed-effects model fit by REML 
Formula: langPOST ~ IQ.ver.cen + avg.IQ.ver.cen + (IQ.ver.cen | schoolNR) 
   Data: bdf 
   AIC   BIC logLik MLdeviance REMLdeviance
 15230 15264  -7609      15214        15218
Random effects:
 Groups   Name        Variance Std.Dev. Corr   
 schoolNR (Intercept)  8.07566 2.84177         
          IQ.ver.cen   0.20801 0.45608  -0.642 
 Residual             41.34967 6.43037         
number of obs: 2287, groups: schoolNR, 131

Fixed effects:
               Estimate Std. Error t value
(Intercept)    40.74958    0.28808  141.45
IQ.ver.cen      2.45983    0.08364   29.41
avg.IQ.ver.cen  1.40893    0.32374    4.35

Correlation of Fixed Effects:
            (Intr) IQ.vr.
IQ.ver.cen  -0.274       
avg.IQ.vr.c  0.029 -0.213
> (fm03 <- lmer(math ~ year+(1|childid)+(1|schoolid), egsingle))
Linear mixed-effects model fit by REML 
Formula: math ~ year + (1 | childid) + (1 | schoolid) 
   Data: egsingle 
   AIC   BIC logLik MLdeviance REMLdeviance
 16767 16795  -8380      16747        16759
Random effects:
 Groups   Name        Variance Std.Dev.
 childid  (Intercept) 0.66992  0.81849 
 schoolid (Intercept) 0.18693  0.43235 
 Residual             0.34700  0.58907 
number of obs: 7230, groups: childid, 1721; schoolid, 60

Fixed effects:
             Estimate Std. Error t value
(Intercept) -0.780482   0.061088  -12.78
year         0.746123   0.005396  138.26

Correlation of Fixed Effects:
     (Intr)
year -0.031
> (fm04 <- lmer(math ~ year+(1|childid)+(year|schoolid), egsingle))
Linear mixed-effects model fit by REML 
Formula: math ~ year + (1 | childid) + (year | schoolid) 
   Data: egsingle 
   AIC   BIC logLik MLdeviance REMLdeviance
 16494 16536  -8241      16472        16482
Random effects:
 Groups   Name        Variance Std.Dev. Corr  
 childid  (Intercept) 0.672304 0.81994        
 schoolid (Intercept) 0.166812 0.40843        
          year        0.011584 0.10763  0.443 
 Residual             0.324516 0.56966        
number of obs: 7230, groups: childid, 1721; schoolid, 60

Fixed effects:
            Estimate Std. Error t value
(Intercept) -0.77709    0.05823  -13.34
year         0.76296    0.01533   49.76

Correlation of Fixed Effects:
     (Intr)
year 0.358 
> (fm05 <- lmer(normexam ~ standLRT + sex + schgend + (1|school), Exam))
Linear mixed-effects model fit by REML 
Formula: normexam ~ standLRT + sex + schgend + (1 | school) 
   Data: Exam 
  AIC  BIC logLik MLdeviance REMLdeviance
 9360 9398  -4674       9326         9348
Random effects:
 Groups   Name        Variance Std.Dev.
 school   (Intercept) 0.085829 0.29297 
 Residual             0.562534 0.75002 
number of obs: 4059, groups: school, 65

Fixed effects:
              Estimate Std. Error t value
(Intercept)  -0.001049   0.055569   -0.02
standLRT      0.559754   0.012450   44.96
sexM         -0.167392   0.034100   -4.91
schgendboys   0.177691   0.113475    1.57
schgendgirls  0.158997   0.089403    1.78

Correlation of Fixed Effects:
            (Intr) stnLRT sexM   schgndb
standLRT    -0.014                      
sexM        -0.316  0.061               
schgendboys -0.395 -0.003 -0.145        
schgendgrls -0.622  0.009  0.197  0.245 
> (fm06 <- lmer(score ~ gcsecnt + (1|school) + (1|lea), Chem97))
Linear mixed-effects model fit by REML 
Formula: score ~ gcsecnt + (1 | school) + (1 | lea) 
   Data: Chem97 
    AIC    BIC logLik MLdeviance REMLdeviance
 141705 141738 -70848     141686       141697
Random effects:
 Groups   Name        Variance Std.Dev.
 school   (Intercept) 1.166197 1.07991 
 lea      (Intercept) 0.014766 0.12151 
 Residual             5.154203 2.27029 
number of obs: 31022, groups: school, 2410; lea, 131

Fixed effects:
            Estimate Std. Error t value
(Intercept)  5.63545    0.03123   180.4
gcsecnt      2.47256    0.01690   146.3

Correlation of Fixed Effects:
        (Intr)
gcsecnt 0.058 
> (fm07 <- lmer(height ~ age + I(age^2) + I(age^3) + I(age^4) + (age + I(age^2)|Subject),
+                  Oxboys))
Linear mixed-effects model fit by REML 
Formula: height ~ age + I(age^2) + I(age^3) + I(age^4) + (age + I(age^2) |      Subject) 
   Data: Oxboys 
   AIC   BIC logLik MLdeviance REMLdeviance
 669.5 707.5 -323.7      644.7        647.5
Random effects:
 Groups   Name        Variance Std.Dev. Corr        
 Subject  (Intercept) 62.74004 7.9209               
          age          2.81669 1.6783   0.607       
          I(age^2)     0.38576 0.6211   0.607 1.000 
 Residual              0.28112 0.5302               
number of obs: 234, groups: Subject, 26

Fixed effects:
            Estimate Std. Error t value
(Intercept) 149.0185     1.5547   95.85
age           6.1754     0.3609   17.11
I(age^2)      1.1310     0.3755    3.01
I(age^3)      0.4523     0.1847    2.45
I(age^4)     -0.3796     0.3414   -1.11

Correlation of Fixed Effects:
         (Intr) age    I(g^2) I(g^3)
age       0.554                     
I(age^2)  0.169  0.276              
I(age^3) -0.001 -0.382  0.027       
I(age^4)  0.024  0.018 -0.912 -0.021
Warning message:
nlminb returned message false convergence (8)
 in: `LMEoptimize<-`(`*tmp*`, value = list(maxIter = 200, tolerance = 1.49011611938477e-08,  
> (fm08 <- lmer(height ~ poly(age,4) + (age + I(age^2)|Subject), Oxboys))
Linear mixed-effects model fit by REML 
Formula: height ~ poly(age, 4) + (age + I(age^2) | Subject) 
   Data: Oxboys 
   AIC   BIC logLik MLdeviance REMLdeviance
 658.3 696.3 -318.2      644.5        636.3
Random effects:
 Groups   Name        Variance Std.Dev. Corr        
 Subject  (Intercept) 62.17477 7.88510              
          age          2.72854 1.65183  0.609       
          I(age^2)     0.35817 0.59847  0.609 1.000 
 Residual              0.28269 0.53169              
number of obs: 234, groups: Subject, 26

Fixed effects:
              Estimate Std. Error t value
(Intercept)   149.5197     1.5818   94.52
poly(age, 4)1  64.5424     3.2613   19.79
poly(age, 4)2   4.2034     0.8508    4.94
poly(age, 4)3   1.2862     0.5317    2.42
poly(age, 4)4  -0.5897     0.5317   -1.11

Correlation of Fixed Effects:
            (Intr) p(,4)1 p(,4)2 p(,4)3
poly(ag,4)1 0.623                      
poly(ag,4)2 0.493  0.770               
poly(ag,4)3 0.000  0.000  0.000        
poly(ag,4)4 0.000  0.000  0.000  0.000 
Warning message:
nlminb returned message function evaluation limit reached without convergence (9)
 in: `LMEoptimize<-`(`*tmp*`, value = list(maxIter = 200, tolerance = 1.49011611938477e-08,  
> (fm09 <- lmer(attain ~ verbal * sex+(1|primary)+(1|second), ScotsSec))
Linear mixed-effects model fit by REML 
Formula: attain ~ verbal * sex + (1 | primary) + (1 | second) 
   Data: ScotsSec 
   AIC   BIC logLik MLdeviance REMLdeviance
 14880 14917  -7434      14843        14868
Random effects:
 Groups   Name        Variance Std.Dev.
 primary  (Intercept) 0.275457 0.52484 
 second   (Intercept) 0.014747 0.12144 
 Residual             4.253113 2.06231 
number of obs: 3435, groups: primary, 148; second, 19

Fixed effects:
            Estimate Std. Error t value
(Intercept) 5.914713   0.076795   77.02
verbal      0.158356   0.003787   41.81
sexF        0.121553   0.072413    1.68
verbal:sexF 0.002593   0.005388    0.48

Correlation of Fixed Effects:
            (Intr) verbal sexF  
verbal       0.177              
sexF        -0.482 -0.178       
verbal:sexF -0.122 -0.680  0.161
> (fm10 <- lmer(attain ~ verbal + sex+(1|primary)+(sex|second), ScotsSec))
Linear mixed-effects model fit by REML 
Formula: attain ~ verbal + sex + (1 | primary) + (sex | second) 
   Data: ScotsSec 
   AIC   BIC logLik MLdeviance REMLdeviance
 14862 14905  -7424      14832        14848
Random effects:
 Groups   Name        Variance Std.Dev. Corr   
 primary  (Intercept) 0.268596 0.51826         
 second   (Intercept) 0.046389 0.21538         
          sexF        0.145830 0.38188  -0.798 
 Residual             4.218513 2.05390         
number of obs: 3435, groups: primary, 148; second, 19

Fixed effects:
            Estimate Std. Error t value
(Intercept) 5.927407   0.086603   68.44
verbal      0.159591   0.002771   57.59
sexF        0.098448   0.114517    0.86

Correlation of Fixed Effects:
       (Intr) verbal
verbal  0.115       
sexF   -0.624 -0.060
>