R/9_cpc.R

#' @export
#'
cpc<- function(r,f){
  n <- sum(anova(f)[,1])+1
  p <- n-df.residual(r)
  #print(n)
  #print(p)
  return((deviance(r)/(deviance(f)/df.residual(f)))-(n-2*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.