WLK.test {EWGoF} | R Documentation |
GoF tests based on the likelihood for the Weibull distribution
Description
Computes the Weibull GoF tests based on the three following statistics: the score, Wald and likelihood ratio GoF tests. These tests include the Weibull distribution in larger statistics and apply a parametric test to the additional parameter.
Usage
WLK.test(x, type = "GG1", funEstimate = "MLE", procedure = "S", nsim = 200, r=0)
Arguments
x |
a numeric vector of data values. |
type |
the type of the test statistic used:
|
funEstimate |
the method used to estimate the two Weibull parameters. "MLE" is the default used method based on the Maximum Likelihood Estimators, "LSE" for the Least Squares Estimators and "ME" for the Moment Estimators. |
procedure |
the procedure used as a default is the score "S". The procedure can be either "W" for the Wald test or "LR" for the test based on the likelihood ratio procedure. |
nsim |
an integer specifying the number of replicates used in Monte Carlo. |
r |
an integer specifying the number of right censored observations. |
Details
The tests are based on different generalized Weibull families: the tests (GG1) and (GG2) are based on the Generalized Gamma distribution, the tests (EW) are based on the Exponentiated Weibull, (PGW) on the Power Generalized Weibull, (MO) on the Marshall-Olkin distribution and (MW) are based on the Modified Weibull distribution. Each family can have nine versions depending on the procedure used (score, Wald or likelihood ratio statistic) and on the parameters estimation methods: maximum likelihood, moment or least squares method, except GG1 which has only three versions using the maximum likelihood estimators.
The tests statistics and
are a combination between two Wald test statistics: PGW with ME (
) and MW with MLE (
) after they are centered with their mean values (
and
) and normalized by their standard deviations (respectively
and
).
The expressions of the statistics and
are as follows:
All the previous tests can be applied to type II right censored samples (simple censoring). The censoring is introduced in the MLEs.
A second statistic similar to
, is combining Wald and likelihood ratio tests based on the Generalized Gamma distribution:
Value
An object of class htest.
Author(s)
Meryam KRIT
References
Krit M., Gaudoin O., Xie M. and Remy E., Simplified likelihood goodness-of-fit tests for the Weibull distribution, Communications in Statistics - Simulation and Computation.
Examples
x <- rlnorm(50,.3)
#Apply some likelihood based tests
WLK.test(x,type="GG1",funEstimate="MLE",procedure="W")
WLK.test(x,type="PGW",funEstimate="ME",procedure="S")
WLK.test(x,type="MO",funEstimate="LSE",procedure="LR")
#Apply G to censored sample at right r=10
a<- sort(x[1:40])
## Not run:
WLK.test(a,type="G",r=10)
## End(Not run)