R/fix0102.R

In EnvStats: Package for Environmental Statistics, Including US EPA Guidance

fix0102 <-
function (d = -1, alpha = 2, beta = 2, mu = 0, sigma = 1, eta = 0, 
    kappa = 1) 
{
    r <- (d - mu)/sigma
    rt <- (d - eta)/kappa
    g20 <- g2.m.singly.censored(r, 0, alpha, beta)
    g1.5 <- g1.m.singly.censored(r, 0.5, alpha, beta)
    s.r <- g1.5/g20
    exp.psi <- pnorm(rt) * s.r - fix0107(fix0106, d, eta + 10 * 
        kappa, 1000, alpha, beta, mu, sigma, eta, kappa) * (2 * 
        alpha + 1)
    exp.psi
}