mgpd (version 1.99)

bgpd_maxlik: internal

Description

internal use only

Usage

bgpd_maxlik(param, dat, model = "log", fixed = FALSE, psi = 3, ...)

Arguments

param

dat

model

fixed

psi

...

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, model = "log", fixed = FALSE, psi = 3, 
    ...) 
{
    mlmax = 1e+15
    models = c("log", "psilog", "philog", "neglog", "psineglog", 
        "phineglog", "bilog", "negbilog", "mix", "ct", "taj", 
        "smith")
    if (!(model %in% models)) 
        stop(paste("'", model, "' is not a valid model.", sep = ""))
    else {
        if (model == "log") 
            ml = ml_log(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "smith") 
            ml = ml_smith(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "psilog") 
            ml = ml_psilog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "philog") 
            ml = ml_philog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "neglog") 
            ml = ml_neglog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "psineglog") 
            ml = ml_psineglog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "phineglog") 
            ml = ml_phineglog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "bilog") 
            ml = ml_bilog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "negbilog") 
            ml = ml_negbilog(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "mix") 
            ml = ml_mix(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "ct") 
            ml = ml_ct(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        if (model == "taj") 
            ml = ml_taj(param = param, dat = dat, mlmax = mlmax, 
                fixed = fixed)
        ml
    }
  }