R/summary.nplr.R

In fredcommo/nplr: N-Parameter Logistic Regression

summary.nplr <- function(object, ...){
    
    pars <- unlist(getPar(object))
    pars[-1] <- format(pars[-1], digits = 6, scientific = TRUE)
    
    goodness <- getGoodness(object)
    gof <- format(goodness$gof, digits = 6, scientific = TRUE)
    wgof <- format(goodness$wgof, digits = 6, scientific = TRUE)
    errors <- t(getStdErr(object))
    auc <- format(getAUC(object), digits = 6)
    infl <- as.numeric(getInflexion(object))
    inflpt <- cbind.data.frame(
        xInfl = format(infl[1], digits = 6),
        yInfl = format(infl[2], digits = 6)
    )
    estim <- getEstimates(object, .5)
    interv <- sprintf("[%s | %s]",
                      format(estim$x.025, digits = 6, scientific = TRUE),
                      format(estim$x.975, digits = 6, scientific = TRUE)
    )
    ic50 <- format(estim$x, digits = 6, scientific = TRUE)
    IC50 <- cbind(IC50 = ic50, "[95%]" = interv)
    nplrVersion <- as.character(packageVersion("nplr"))
    nplrDate <- as.character(packageDescription("nplr")["Date"])
    Rversion <- as.character(version["version.string"])
    
    out <- cbind.data.frame(t(pars),
                            GOF = gof,
                            weightedGOF = wgof,
                            errors,
                            auc,
                            inflpt,
                            "Log10(IC50)" = log10(estim$x),
                            IC50,
                            "date (Y-m-d)" = format(Sys.Date(), "%Y-%m-%d"),
                            "nplr version" = sprintf("%s (%s)",nplrVersion, nplrDate),
                            "R version" = gsub("R version ", "", Rversion)
    )
    
    out <- data.frame(summary = t(out))
    colnames(out) <- 'value'
    print(out)
}