R/gamobs.22.R

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

gamobs.22 <-
function (x, mu = 0, sig = 1, eta = 0, kappa = 1, alf = 2, bet = 2) 
{
    r <- (x - mu)/sig
    1/sig^2 * (1 + r^2 * (s.3.m.singly.censored(r, alf, bet) - 
        s.2.m.singly.censored(r, alf, bet)^2) + 2 * r * s.2.m.singly.censored(r, 
        alf, bet)) * dnorm(x, mean = eta, sd = kappa)
}