R/I_getVp.R

Defines functions .getVp

.getVp <- function(.fit, .obj, .lsp, .lpi)
{
  .Vp <- if(inherits(.obj$family, "general.family")){
    .fit$gcv.ubre <- as.numeric(.fit$REML) 
    .fit$outer.info <- NULL
    .fit$sp <- exp(.lsp)
    .fit$scale.estimated <- FALSE
    .fit$scale <- 1
    .fit$method <- "REML"
    .Vp <- gam.fit5.post.proc(.fit,.obj$Sl,.obj$L,.obj$lsp0,.obj$S,.obj$off)$Vb
    .Vp <- .Vp[.lpi[[1]], .lpi[[1]]]
  } else {
    .Vp <- .fit$Vb
  } 
  return( .Vp )
}

Try the qgam package in your browser

Any scripts or data that you put into this service are public.

qgam documentation built on Nov. 23, 2021, 1:07 a.m.