kinreport <- function(kinobject, file = NA, vcov = FALSE, endpoint.digits = 1)
{
	if (!is.na(file)) {
		sink(file, split=TRUE)
	}

	cat("Parent compound: ", kinobject$parent, "\n")
        if (!is.null(kinobject$label)) cat("Label position:\t\t", kinobject$label, "\n")
	cat("Study type:      ", kinobject$type, "\n")
	cat("System:          ", kinobject$system, "\n")
	if (!is.null(kinobject$source)) {
          cat("Source:          ", kinobject$source, "\n")
        }
	cat("\n")
	fit.names <- names(kinobject$fits)
	for (kinmodel in fit.names)
	{
                m <- kinobject$fits[[kinmodel]]
                if (!(class(m) == "try-error")) {
                    cat("\n\n---\n")
                    cat("Nonlinear least squares fit of the", kinmodel, "model\n\n")
                    cat("Parameter estimation:\t")
                    s <- summary(m)
                    df <- s$df[2]
                    p <- 1 - pt(s$parameters[,3], df = df)
                    parms <- cbind(s$parameters[,c(1,2,3)], "Pr(>t)" = p)
                    cat("\n")
                    print(parms, digits=3)
                    cat("\n")
                    if(vcov)
                    {
                        cat("Variance-covariance matrix:\n")
                        print(vcov(m))
                        cat("\n")
                    }
                    cat("Chi2 error estimation:\t", 
                            round(100 * kinobject$results$stats[kinmodel, "err.min"], digits=2), 
                            " %\n", sep="")
                    cat("\n")
                }
	}
	cat("\n\n---\n")
	cat("Endpoint estimates\n\n")
	print(round(kinobject$results$results, digits=endpoint.digits))

	if (!is.na(file)) sink()
}