ml_phineglog: internal

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/ml_phineglog.R

Description

internal use only

Usage

1
ml_phineglog(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 = "phineglog", 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
  }


mgpd documentation built on May 19, 2017, 11:06 p.m.
Search within the mgpd package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.