SCM

SCM Repository

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

Diff of /pkg/R/lmer.R

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

revision 752, Sun May 29 02:14:58 2005 UTC revision 754, Mon May 30 18:40:45 2005 UTC
# Line 693  Line 693 
693                               control = list(trace = getOption("verbose"),                               control = list(trace = getOption("verbose"),
694                               iter.max = controlvals$msMaxIter))                               iter.max = controlvals$msMaxIter))
695                    optpars <- optimRes$par                    optpars <- optimRes$par
696                      if (optimRes$convergence != 0)
697                          warning("nlminb failed to converge")
698                } else {                } else {
699                    optimRes <-                    optimRes <-
700                        optim(fn =                        optim(c(fixef(obj),
701                              c(fixef(obj), .Call("lmer_coef", obj, TRUE,                                .Call("lmer_coef", obj, TRUE, PACKAGE = "Matrix")),
                                                 PACKAGE = "Matrix")),  
702                              devLaplace,                              devLaplace,
703                              method = "BFGS", hessian = TRUE,                              method = "BFGS", hessian = TRUE,
704                              control = list(trace = getOption("verbose"),                              control = list(trace = getOption("verbose"),
705                              reltol = controlvals$msTol,                              reltol = controlvals$msTol,
706                              maxit = controlvals$msMaxIter))                              maxit = controlvals$msMaxIter))
707                    optpars <- optimRes$par                    optpars <- optimRes$par
                   Hessian <- optimRes$hessian  
               }  
708                if (optimRes$convergence != 0)                if (optimRes$convergence != 0)
709                    warning("optim failed to converge")                    warning("optim failed to converge")
710                      Hessian <- optimRes$hessian
711                  }
712    
713                ##fixef(obj) <- optimRes$par[seq(length = responseIndex - 1)]                ##fixef(obj) <- optimRes$par[seq(length = responseIndex - 1)]
714                if (getOption("verbose")) {                if (getOption("verbose")) {

Legend:
Removed from v.752  
changed lines
  Added in v.754

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