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 1167, Sun Jan 15 21:13:54 2006 UTC revision 1168, Mon Jan 16 14:40:39 2006 UTC
# Line 609  Line 609 
609            if (useScale)            if (useScale)
610                sc <- .Call("mer_sigma", x, REML, PACKAGE = "Matrix")                sc <- .Call("mer_sigma", x, REML, PACKAGE = "Matrix")
611            sc2 <- sc * sc            sc2 <- sc * sc
612            ans <- lapply(x@Omega, function(el) {            cnames <- x@cnames
613                el <- as(sc2 * solve(el), "dpoMatrix")            ans <- x@Omega
614              for (i in seq(a=ans)) {
615                  el <- as(sc2 * solve(ans[[i]]), "dpoMatrix")
616                  el@Dimnames <- list(cnames[[i]], cnames[[i]])
617                el@factors$correlation <- as(el, "correlation")                el@factors$correlation <- as(el, "correlation")
618                el                ans[[i]] <- el
619            })            }
620            attr(ans, "sc") <- sc            attr(ans, "sc") <- sc
621            ans            ans
622        })        })

Legend:
Removed from v.1167  
changed lines
  Added in v.1168

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