internal

Share:

Description

internal use only

Usage

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

Arguments

param
dat
mlmax
fixed
...

Details

internal use only

Value

internal use only

Note

internal use only

Author(s)

P. Rakonczai

References

internal use only

See Also

internal use only

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
##---- 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 = "psilog", mar1 = param[1:3], 
        mar2 = param[4:6], dep = param[7], asy = param[8], p = param[9]), 
        silent = TRUE)
    if (!is.null(lik)) {
        if (is.null(attr(lik, "class"))) {
            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
  }

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.