Displaying 1 result from an estimated 1 matches for "lls2".
Did you mean:
lls
2017 Mar 07
0
Potential clue for Bug 16975 - lme fixed sigma - inconsistent REML estimation
...val <- lapply(object, Initialize, data)
attr(val, "plen") <- unlist(lapply(val, function(el) length(coef(el))))
class(val) <- c("varSum","varComb", "varFunc")
val
}
logLik.varSum <-
function(object, ...)
{
lls <- lapply(object, logLik)
lls2 <- apply(as.data.frame(lapply(object, varWeights)), 1, function(x){
1/sqrt(sum((1/x)^2))
})
val <- sum(log(lls2))
attr(val, "df") <- sum(unlist(lapply(lls, attr, "df")))
class(val) <- "logLik"
val
}
####* The methods from here to the example...