R/aggregation.misclass.R

In peperr: Parallelised Estimation of Prediction Error

aggregation.misclass <- function (full.data = NULL, response, x, model, cplx = NULL,
    type = c("apparent", "noinf"), fullsample.attr = NULL, ...)
{
    data <- as.data.frame(x)
    data$response <- response
    if (class(model)[1] == "penfit") {
        probs <- predict(model, data = data, penalized = x, ...)
    }  else {
        if(class(model)[1] == "glm") {
            probs <- predict(model, newdata = data, penalized = x, , type = "response", ...)
        } else {
        probs <- predict(model, data = data, type = "response",
            ...)
            }
    }
    type <- match.arg(type)
    if (type == "apparent") {
        mr <- sum(abs(round(probs) - response))/length(response)
    }
    if (type == "noinf") {
        mr <- mean(abs((matrix(response, length(response), length(response),
            byrow = TRUE) - round(probs))))
    }
    mr
}