R/print.summary.dbglm.r

In dbstats: Distance-Based Statistics

 #############################
 #### print.summary.dbglm ####
 #############################


print.summary.dbglm <-function (x, digits = 2, ...){

    # print the call
   x$call[[1]]<-as.name("dbglm")
   cat("\nCall:  ", paste(deparse(x$call), sep = "\n", collapse = "\n"),
        "\n", sep = "")

    # Deviance Residuals
    cat("\nDeviance Residuals:\n")
    print(summary(as.numeric(format(round(x$deviance.resid,digits=4)))))
    
    # dispersion
    if (x$family$family %in% c("poisson","binomial"))
     cat(gettextf("\n(Dispersion parameter for %s family taken to be %i)",
            x$family$family,x$dispersion),"\n")
    else
     cat(gettextf("\n(Dispersion parameter for %s family taken to be %f)",
          x$family$family,x$dispersion),"\n")
    
    # deviance: null or estimated model
    cat(gettextf("\n    Null deviance: %s  on %i degrees of freedom",
          format(round(x$null.deviance,digits)),x$df.null))
    cat(gettextf("\nResidual deviance: %s  on %i degrees of freedom",
          format(round(x$deviance,digits=digits)),x$df.residual))
    
    
    # Number of Fisher Scoring iterations and Convergence criterio
    cat(gettextf("\n\nNumber of Fisher Scoring iterations: %i",x$iter),"\n")
    cat(gettextf("Convergence criterion: %s",x$convcrit),"\n\n")
 
   # AIC
   cat("AIC:\t   ", format(round(x$aic, digits)), "\n")
   cat("BIC:\t   ", format(round(x$bic, digits)), "\n")
   cat("GCV:\t   ", format(round(x$gcv, digits)),"\n\n")
 
 }