R/9_SBCp.R

#' @export
SBCp<- function(model){
  n <- sum(anova(model)[,1])+1
  p <- n-df.residual(model)
  #print(n)
  #print(p)
  return(n*log(deviance(model))-n*log(n)+(log(n))*p)

}

Try the ALSM package in your browser

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

ALSM documentation built on May 2, 2019, 10:19 a.m.