startbernweibull: Rough parameter estimate for the Bernoulli-Weibull...

In qmap: Statistical Transformations for Post-Processing Climate Model Output

Description

Estimates rough starting values for the Bernoulli-Weibull distribution using the method of moments for the shape and the scale parameters. The probability of non-zero events is estimated as the fraction of values that are larger than zero.

Usage

startbernweibull(x)

Arguments

x	numeric vector.

Value

A list containing:

prob	probability of non-zero event.
scale	scale parameter of the weibull distribution.
shape	shape parameter of the weibull distribution.

Note

In this package startbernweibull is intended to be used in conjunction with fitQmapDIST (and mledist) with parameter distr="bernweibull".

Author(s)

Lukas Gudmundsson

See Also

fitQmapDIST, bernweibull, fitdist

Examples

gg <- rbernweibull(n=300, prob=0.2, scale=1, shape=1)
startbernweibull(gg)
mledist(gg,"bernweibull",startbernweibull(gg))

Example output

Loading required package: fitdistrplus
Loading required package: MASS
Loading required package: survival
Loading required package: npsurv
Loading required package: lsei
$prob
[1] 0.1966667

$shape
[1] 1.092038

$scale
[1] 1.008519

$estimate
     prob     shape     scale 
0.1966526 1.1012666 0.9982891 

$convergence
[1] 0

$value
[1] 204.9969

$hessian
               prob         shape     scale
prob   1.899153e+03 -7.105427e-09   0.00000
shape -7.105427e-09  8.957628e+01 -25.66546
scale  0.000000e+00 -2.566546e+01  71.79244

$optim.function
[1] "optim"

$optim.method
[1] "Nelder-Mead"

$fix.arg
NULL

$fix.arg.fun
NULL

$weights
NULL

$counts
function gradient 
      62       NA 

$optim.message
NULL

$loglik
[1] -204.9969

qmap documentation built on May 1, 2019, 7:31 p.m.