R/summary.hyperbFit.R

In HyperbolicDist: The Hyperbolic Distribution

### Calculate summary for hyperbFit object
### Only applies when hessian is asked for
###
### DJS 11/08/06
summary.hyperbFit <- function(object, ...)
{
    if (!inherits(object,"hyperbFit")){
        stop("Object must belong to class hyperbFit")
    }
    if (!is.null(object$hessian)){
        sds <- sqrt(diag(solve(object$hessian)))
        sds[2] <- object$Theta[2]*sds[2]
        sds[3] <- object$Theta[3]*sds[3]
        names(sds) <- c("pi","zeta","delta","mu")
        object$sds <- sds
    }
    class(object) <- "summary.hyperbFit"
    return(object)
} ## End of summary.hyperbFit

### Print summary
print.summary.hyperbFit <- function(x,
                                    digits = max(3, getOption("digits") - 3),
                                    ...)
{
    if (!inherits(x,"summary.hyperbFit")){
        stop("Object must belong to class summary.hyperbFit")
    }
    cat("\nData:     ", x$obsName, "\n")
    cat("Parameter estimates:\n")
    if (is.null(x$sds)){
        print.default(format(x$Theta, digits = digits),
                      print.gap = 2, quote = FALSE)
    }else{
        ans <- format(rbind(x$Theta, x$sds), digits = digits)
        ans[1, ] <- sapply(ans[1, ], function(obs) paste("", obs))
        ans[2, ] <- sapply(ans[2, ],
                           function(obs) paste("(", obs, ")", sep = ""))
        dn <- dimnames(ans)
        dn[[1]] <- rep("", 2)
        dn[[2]] <-
            paste(substring("      ", 1,
                            (nchar(ans[2, ]) - nchar(dn[[2]]))%/%2), dn[[2]])
        dn[[2]] <- paste(dn[[2]],
                         substring("      ", 1,
                                   (nchar(ans[2, ]) - nchar(dn[[2]]))%/%2))
        dimnames(ans) <- dn
        print.default(ans, print.gap = 2, quote = FALSE)
    }
    cat("Likelihood:        ",x$maxLik,"\n")
    cat("Method:            ",x$method,"\n")
    cat("Convergence code:  ",x$conv,"\n")
    cat("Iterations:        ",x$iter,"\n")
    invisible(x)
}