mgpd (version 1.99)

ml_log: internal

Description

internal use only

Usage

ml_log(param, dat, mlmax = 1e+15, fixed = FALSE, ...)

Arguments

param

dat

mlmax

fixed

...

Value

• internal use only

Details

internal use only

References

internal use only

See Also

internal use only

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (param, dat, mlmax = 1e+15, fixed = FALSE, ...) 
{
    loglik = mlmax
    lik = NULL
    x = dat[, 1]
    y = dat[, 2]
    if (fixed) 
        param[1] = 0
    lik = try(dbgpd(x, y, model = "log", mar1 = param[1:3], mar2 = param[4:6], 
        dep = param[7]))
    if (!is.null(lik)) {
        loglik = -sum(log(lik))
        if (min(1 + param[3] * (x - param[1])/param[2]) < 0) 
            loglik = mlmax
        if (min(1 + param[6] * (y - param[4])/param[5]) < 0) 
            loglik = mlmax
    }
    loglik
  }