1 files changed, 47 insertions, 0 deletions
diff --git a/R/kinreport.R b/R/kinreport.R
new file mode 100644
index 0000000..4156803
--- /dev/null
+++ b/R/kinreport.R
@@ -0,0 +1,47 @@
+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()
+}