~ubuntu-branches/ubuntu/natty/lme4/natty

« back to all changes in this revision

Viewing changes to tests/lmer-1.Rout.save

  • Committer: Bazaar Package Importer
  • Author(s): Dirk Eddelbuettel
  • Date: 2010-11-11 11:00:44 UTC
  • mfrom: (1.1.23 upstream) (2.2.9 sid)
  • Revision ID: james.westby@ubuntu.com-20101111110044-dt302242aebr14fm
Tags: 0.999375-37-1
New upstream release

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
1
 
2
 
R version 2.12.0 Patched (2010-10-19 r53372)
 
2
R version 2.12.0 (2010-10-15)
3
3
Copyright (C) 2010 The R Foundation for Statistical Computing
4
4
ISBN 3-900051-07-0
5
 
Platform: i686-pc-linux-gnu (32-bit)
 
5
Platform: x86_64-pc-linux-gnu (64-bit)
6
6
 
7
7
R is free software and comes with ABSOLUTELY NO WARRANTY.
8
8
You are welcome to redistribute it under certain conditions.
28
28
 1756 1775 -871.8     1752    1744
29
29
Random effects:
30
30
 Groups   Name        Variance Std.Dev. Corr  
31
 
 Subject  (Intercept) 612.092  24.7405        
 
31
 Subject  (Intercept) 612.090  24.7405        
32
32
          Days         35.072   5.9221  0.066 
33
33
 Residual             654.941  25.5918        
34
34
Number of obs: 180, groups: Subject, 18
153
153
5    1.4381076
154
154
6   -1.7594785
155
155
7   -0.2546842
156
 
8    0.6977648
 
156
8    0.6977649
157
157
9    1.2121675
158
 
10   4.9434830
 
158
10   4.9434829
159
159
 
160
160
> (sf1 <- summary(fit.1)) # show() is as without summary()
161
161
Linear mixed model fit by REML 
446
446
 101.5 107.5 -44.77    87.59   89.54
447
447
Random effects:
448
448
 Groups   Name        Variance   Std.Dev. Corr   
449
 
 group    (Intercept) 140.641521 11.85924        
 
449
 group    (Intercept) 140.646168 11.85943        
450
450
          x             0.014427  0.12011 -1.000 
451
 
 Residual               2.722977  1.65014        
 
451
 Residual               2.722972  1.65014        
452
452
Number of obs: 20, groups: group, 2
453
453
 
454
454
Fixed effects:
455
455
            Estimate Std. Error t value
456
 
(Intercept) 15.23544    8.40776   1.812
 
456
(Intercept) 15.23553    8.40752   1.812
457
457
x            0.18559    0.08508   2.181
458
458
 
459
459
Correlation of Fixed Effects:
471
471
[1] 0 4
472
472
 
473
473
Slot "ST":
474
 
            [,1]
475
 
[1,]  7.18678678
476
 
[2,]  0.00000000
477
 
[3,] -0.01012811
 
474
              [,1]
 
475
[1,]  7.186912e+00
 
476
[2,]  5.962606e-08
 
477
[3,] -1.012811e-02
478
478
 
479
479
Slot "call":
480
480
lmer(formula = y ~ x + (x | group))
481
481
 
482
482
Slot "deviance":
483
 
[1] 87.58735
 
483
[1] 87.58728
484
484
 
485
485
Slot "dims":
486
486
   nt     n     p     q     s    np   LMM  REML  fTyp  lTyp  vTyp  nest useSc 
490
490
 
491
491
Slot "fixef":
492
492
                  [,1]
493
 
(Intercept) 15.2354383
494
 
x            0.1855865
 
493
(Intercept) 15.2355284
 
494
x            0.1855855
495
495
 
496
496
Slot "nc":
497
497
[1] 2
498
498
 
499
499
Slot "ranef":
500
 
          [,1]
501
 
[1,] -8.350427
502
 
[2,]  8.350427
503
 
[3,]  0.084574
504
 
[4,] -0.084574
 
500
            [,1]
 
501
[1,] -8.35085457
 
502
[2,]  8.35085457
 
503
[3,]  0.08457838
 
504
[4,] -0.08457838
505
505
 
506
506
Slot "sigma":
507
507
         [,1]
508
 
[1,] 1.650145
 
508
[1,] 1.650143
509
509
 
510
510
511
511
> if (FALSE) {  # mcmcsamp for  glmer  not yet available
632
632
m0  3 1034.17 1044.61 -514.08                         
633
633
m1  4  585.49  599.41 -288.75 450.67      1  < 2.2e-16
634
634
 
635
> ## Check that quasi families throw an error
 
636
> try(gm1 <- lmer(cbind(incidence, size - incidence) ~ period + (1|herd),
 
637
+                 data = cbpp, family = quasibinomial))
 
638
Error in glmer(formula = cbind(incidence, size - incidence) ~ period +  : 
 
639
  "quasi" families cannot be used in glmer
 
640
> try(gm1 <- lmer(incidence ~ period + (1|herd),
 
641
+                 data = cbpp, family = quasipoisson))
 
642
Error in glmer(formula = incidence ~ period + (1 | herd), data = cbpp,  : 
 
643
  "quasi" families cannot be used in glmer
 
644
> try(gm1 <- lmer(incidence ~ period + (1|herd),
 
645
+                 data = cbpp, family = quasi))
 
646
Error in glmer(formula = incidence ~ period + (1 | herd), data = cbpp,  : 
 
647
  "quasi" families cannot be used in glmer
 
648
635
649
> if(.unsafe.BLAS) rm(identical)
636
650
637
651
> cat('Time elapsed: ', proc.time(),'\n') # for ``statistical reasons''
638
 
Time elapsed:  6.609 0.185 7.371 0 0 
 
652
Time elapsed:  9.55 0.08 9.636 0 0 
639
653