Nothing
## deviance ##
deviance.nrm <-
function(object, ...)
{
### number of parameters
nme <- length(object$erg_distr$mean_est) - 1
#object$ctrl$sigmaest
nva <- object$ctrl$sigmaest * (length(object$erg_distr$sig_est) -1)
npar <- ncol(object$reshOBJ$Qmat) + nme + nva - length(object$ctrl$Clist)
structure(2*object$last_mstep$value, df=npar)
}
## logLik ##
logLik.nrm <-
function(object,...)
{
### number of parameters
nme <- length(object$erg_distr$mean_est) - 1
nva <- object$ctrl$sigmaest * (length(object$erg_distr$sig_est) -1)
npar <- ncol(object$reshOBJ$Qmat) + nme + nva - length(object$ctrl$Clist)
# number of observations
nobs <- sum(sapply(object$reshOBJ$d,nrow))
return(structure((-1)*object$last_mstep$value, df=npar, nobs=nobs))
}
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.