R/BIC.siland.R

In siland: Spatial Influence of Landscape

Defines functions BIC.siland

Documented in BIC.siland

BIC.siland<-function(object,...)
{
  N=length(object$fitted)
  loglik=object$loglik
  nparam=length(object$coefficients)

  if(object$family$family!="gaussian")
  BIC  =-2*loglik +nparam*log(N)
  else
  BIC=-2*loglik +(nparam+1)*log(N)

  cat(paste("BIC = ",signif(BIC,digits=6)))
  invisible(BIC)
}

Try the siland package in your browser