SCM

SCM Repository

[matrix] Diff of /pkg/tests/lmer.R
ViewVC logotype

Diff of /pkg/tests/lmer.R

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1161, Fri Jan 13 23:46:29 2006 UTC revision 1162, Fri Jan 13 23:47:38 2006 UTC
# Line 1  Line 1 
1  library(Matrix)  library(Matrix)
2  options(show.signif.stars = FALSE)  options(show.signif.stars = FALSE)
3    
4  (fm1 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos),  (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
5               OrchardSprays))  (fm2 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy, method = "ML"))
 (fm2 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos),  
              OrchardSprays, method = "ML"))  
6    
7  ## should produce a warning but fit by REML  ## should produce a warning but fit by REML
8  (fm1b <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos),  (fm1b <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy, method = "AGQ"))
               OrchardSprays, method = "AGQ"))  
9    
10  ## generalized linear mixed model  ## generalized linear mixed model
11  (fm3 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos),  (fm3 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos),
# Line 32  Line 29 
29  ## fit and summarize the model  ## fit and summarize the model
30  fit.1 <- lmer (y ~ 1 + (1 | group.id))  fit.1 <- lmer (y ~ 1 + (1 | group.id))
31  coef (fit.1)# failed in Matrix 0.99-6 -- FIXME: should get a show() method  coef (fit.1)# failed in Matrix 0.99-6 -- FIXME: should get a show() method
32  summary(fit.1)  try(summary(fit.1)) # FIXME: add a summary method
   
33    
34  ## Many family = binomial cases  ## Many family = binomial cases
35  if (isTRUE(try(data(Contraception, package = 'mlmRev')) ==  if (isTRUE(try(data(Contraception, package = 'mlmRev')) ==

Legend:
Removed from v.1161  
changed lines
  Added in v.1162

root@r-forge.r-project.org
ViewVC Help
Powered by ViewVC 1.0.0  
Thanks to:
Vienna University of Economics and Business Powered By FusionForge