Nothing
check_if_order_implemented <- function(family, order) {
if(order > 2L) {
stop("order-", order, " Laplace approximation not yet implemented",
call. = FALSE)
}
if(order > 1L & !is_canonical(family)) {
stop("order-", order, " Laplace approximation not yet implemented for non-canonical link functions")
}
}
find_lfun_Laplace <- function(modfr, devfun_laplace_1, order) {
if(order == 1L ) {
lfun <- function(x) { -devfun_laplace_1(x) / 2 }
}
else {
check_if_order_implemented(modfr$family, order)
if(order == 2L) {
lfun <- function(x){
epsilon_1 <- find_epsilon_1(x, modfr, devfun_laplace_1)
-devfun_laplace_1(x) / 2 + epsilon_1
}
}
}
return(lfun)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.