https://github.com/cran/IQCC
Raw File
Tip revision: 53c14e2d35deeae90e359b99cad0db06d573eb53 authored by Emanuel P. Barbosa on 10 March 2010, 00:00:00 UTC
version 0.5
Tip revision: 53c14e2
alpha.risk.R
alpha.risk <- function(n)
{
    D1 <- function(n)
    {
        d1 <- max(0, d2(n) - 3 * d3(n))
        return(d1)
    }
    D2 <- function(n)
    {
        D2 <- d2(n) + 3 * d3(n)
        return(D2)
    }
    risco <- function(n)
    {
        risco <- 1 - (ptukey(D2(n), n, Inf) - ptukey(D1(n), n, Inf))
        return(risco)
    }
    risk <- rep(0, length(n))
    for(i in 1:length(n))
        risk[i] <- risco(n[i])
    return(risk)
}
back to top