https://github.com/cran/IQCC
Raw File
Tip revision: ef4202fc61e884dce81758e0f31c6f1ae429e523 authored by Emanuel P. Barbosa on 10 March 2010, 00:00:00 UTC
version 0.6
Tip revision: ef4202f
d3.R
d3 <- function(n)
{
    d2 <- d2(n)
    e <- vector()
    for(i in 1:length(n))
    {
        int <- integrate(function(w) { w * (1 - ptukey(w, n[i], Inf)) }, 0, Inf)
        e <- append(e, sqrt(2 * int[[1]] - (d2[1]) ^ 2))
    }
    return(e)
}
back to top