Description Usage Arguments Details Value References See Also Examples
Bias correction of second order of the maximum likelihood estimators of the parameters of the beta regression model.
1 | bias(fit)
|
fit |
Fit beta regression models for rates and proportions via maximum likelihood using a parametrization with mean (depending through a link function on the covariates) and precision parameter (called phi). |
The parameters of the beta regression model are estimated by the maximum likelihood method (see Ferrari and Cribari-Neto, 2004). These estimators are generally biased in models that use link function. This bias is not a serious problem when the sample size is large, however, when the sample is small, this bias can be large compared with the standard-error estimator.
Simas et al (2010) defines formulas general for second-order biases of the beta regression model with constant or variable-precision accuracy.
bias() returns a matrix with corrected coefficients.
Ferrari, S.L.P., and Cribari-Neto, F. (2004). Beta Regression for Modeling Rates and Proportions. Journal of Applied Statistics, 31(7), 799-815.
Simas, A.B., Barreto-Souza, W., and Rocha, A.V. (2010). Improved Estimators for a General Class of Beta Regression Models. Computational Statistics and Data Analysis, 54(2), 348-366.
betareg
1 2 3 4 5 6 7 8 9 10 11 | require(betareg)
## Section 4 from Ferrari and Cribari-Neto (2004)
data("GasolineYield", package = "betareg")
bbt <- betareg(yield ~ batch + temp, data = GasolineYield)
bias(bbt)
## Section 3 from online supplements to Simas et al. (2010)
## mean model as in gy above
## precision model with regressor temp
bbt2 <- betareg(yield ~ batch + temp | temp, data = GasolineYield)
bias(bbt2)
|
Loading required package: betareg
[,1] [,2]
(Intercept) "(Intercept)" "-6.15747245892992"
batch1 "batch1" "1.7261925685144"
batch2 "batch2" "1.32245170818145"
batch3 "batch3" "1.57183129174859"
batch4 "batch4" "1.05830812846942"
batch5 "batch5" "1.13255313746645"
batch6 "batch6" "1.03926017086979"
batch7 "batch7" "0.542829115983543"
batch8 "batch8" "0.495389514347454"
batch9 "batch9" "0.386716566211749"
temp "temp" "0.0109684071678406"
(phi) "(phi)" "497.836056240778"
[,1] [,2]
(Intercept) "(Intercept)" "-5.92776986572297"
batch1 "batch1" "1.60213076457744"
batch2 "batch2" "1.30015030359806"
batch3 "batch3" "1.56803555614975"
batch4 "batch4" "1.03024697021589"
batch5 "batch5" "1.15464613537861"
batch6 "batch6" "1.01990210932617"
batch7 "batch7" "0.622505608783371"
batch8 "batch8" "0.564175434468378"
batch9 "batch9" "0.359759380595258"
temp "temp" "0.0103717445271745"
(Intercept) "(phi)_(Intercept)" "1.05408390253782"
temp "(phi)_temp" "0.0158182304828595"
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.