R/predict.humpfit.R

In jarioksa/natto: An Extreme 'vegan' Package of Experimental Code

#' @param newdata Values of \code{mass} used in \code{predict}. The
#'     original data values are used if missing.
#' @importFrom stats fitted coef
#' @rdname humpfit
#' @export
`predict.humpfit` <-
    function(object, newdata = NULL, ...)
{
    if (is.null(newdata))
        return(fitted(object))
    else {
        p <- coef(object)
        x <- unlist(newdata)
        x <- ifelse(x < p[1], x/p[1], p[1]*p[1]/x/x)
        fv <- p[3]*log(1 + p[2]*x/p[3])
    }
    fv
}