se {ParetoPosStable}

Approximated standard errors of Pareto Positive Stable (PPS) parameter estimates

Description

It approximates the stantard errors of PPS parameter estimates by bootstrapping.

Usage

se(PPSfit, k = 2000, parallel = TRUE, ncores = 2, ...)

Arguments

PPSfit

a PPSfit Object, typically from PPS.fit().

k

the number of steps in the bootstrapping procedure.

parallel

A logical argument specifying if parallelization is allowed in the bootstrapping procedure.

ncores

is the number of cores that we use if parallel is TRUE

...

other arguments.

Details

The function simulates k samples from the model given in the PPSfit argument, fits them with the same method of estimation and uses the parameter estimates to approximate the standard errors.

Value

A list with the standard errors.

References

Sarabia, J.M and Prieto, F. (2009). The Pareto-positive stable distribution: A new descriptive model for city size data, Physica A: Statistical Mechanics and its Applications, 388(19), 4179-4191.

See Also

PPS.fit

Examples

x <- rPPS(50, 1.2, 100, 2.3)
fit <- PPS.fit(x)
coef(fit)
se(fit, k = 50)