Nothing
###Summary Method for skewhypFit Objects #################################
summary.skewhypFit <- function(object,...){
if (!inherits(object,"skewhypFit"))
stop("Object must belong to class skewhypFit")
if (!is.null(object$hessian)){
sds <- sqrt(diag(solve(object$hessian)))
sds[2] <- object$param[2]*sds[2]
sds[4] <- object$param[4]*sds[4]
names(sds) <- c("mu","delta","beta","nu")
object$sds <- sds
}
class(object) <- "summary.skewhypFit"
return(object)
}
### Print Summary #####################################################
print.summary.skewhypFit <-
function(x, digits = max(3, getOption("digits") - 3), ...)
{
if (!inherits(x,"summary.skewhypFit")){
stop("Object must belong to class summary.skewhypFit")
}
cat("\nData: ", x$obsName, "\n")
cat("Parameter estimates:\n")
if (is.null(x$sds)){
print.default(format(x$param, digits = digits),
print.gap = 2, quote = FALSE)
}else{
ans <- format(rbind(x$param, 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)
}
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.