XploRe Help : estWeibull

Keywords - Function groups - @ A B C D E F G H I J K L M N O P Q R S T U V W X Y Z

Library:	stats
See also:	pdfWeibull cdfWeibull

Quantlet:		estWeibull
Description:		computes estimators of parameters (beta, tau) of Weibull distribution.

Reference(s):

K. Burnecki, A. Misiorek, R.Weron (2004) "Loss distibution", in “Statistical Tools for Finance and Insurance”, eds. P. Cizek, W. Härdle, R. Weron, Springer.

Usage:	`{db,de,vb,ve,w2b,w2e,a2b,a2e,beta,tau,betamin,taumin,likconv} = estWeibull(analyzedsample,method,numberofiterations)`
Input:
	`analyzedsample`	n x 1 vector, sample for which Weibull distribution parameters will be estimated.
	`method`	integer, estimation method selection flag (0 - all methods, 1 - ML, 2 - A2).
	`numberofiterations`	integer, number of iterations for minimization procedures.
Output:
	`db`	scalar, Value of Kolmogorov statistic for ML parameters.
	`de`	scalar, Value of Kolmogorov statistic for A2 parameters.
	`vb`	scalar, Value of Kuiper statistic for ML parameters.
	`ve`	scalar, Value of Kuiper statistic for A2 parameters.
	`w2b`	scalar, Value of Cramer-von Mises statistic for ML parameters.
	`w2e`	scalar, Value of Cramer-von Mises statistic for A2 parameters.
	`a2b`	scalar, Value of Anderson-Darling statistic for ML parameters.
	`a2e`	scalar, Value of Anderson-Darling statistic for A2 parameters.
	`beta`	scalar, Parameter beta from ML estimation.
	`tau`	scalar, Parameter tau from ML estimation.
	`betamin`	scalar, Parameter beta from A2 estimation.
	`taumin`	scalar, Parameter tau from A2 estimation.
	`likconv`	scalar, 1 - if ML estimation converged, 0 - otherwise.

Example:

library("xplore") library("stats") library("distribs") library("nummath") randomize(101) sample=rndWeibull(1000,1,2,5) estWeibull(sample,0,3000)

Result:

Contents of _tmp.db
[1,]  0.41295
Contents of _tmp.de
[1,]  0.013435
Contents of _tmp.vb
[1,]  0.63337
Contents of _tmp.ve
[1,]  0.023236
Contents of _tmp.w2b
[1,]   57.763
Contents of _tmp.w2e
[1,]  0.018431
Contents of _tmp.a2b
[1,]   276.49
Contents of _tmp.a2e
[1,]  0.16705
Contents of _tmp.beta
[1,]   2.0232
Contents of _tmp.tau
[1,]    5.085
Contents of _tmp.betamin
[1,]   2.0307
Contents of _tmp.taumin
[1,]   5.0745
Contents of _tmp.likconv
[1,]        1

Author: A. Misiorek, 20041110 license MD*Tech