https://github.com/cran/IQCC
Tip revision: ef4202fc61e884dce81758e0f31c6f1ae429e523 authored by Emanuel P. Barbosa on 10 March 2010, 00:00:00 UTC
version 0.6
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)
}