Package: brglm2 1.1.0

brglm2: Bias Reduction in Generalized Linear Models

Estimation and inference from generalized linear models based on various methods for bias reduction and maximum penalized likelihood with powers of the Jeffreys prior as penalty. The 'brglmFit()' fitting method can achieve reduction of estimation bias by solving either the mean bias-reducing adjusted score equations in Firth (1993) <doi:10.1093/biomet/80.1.27> and Kosmidis and Firth (2009) <doi:10.1093/biomet/asp055>, or the median bias-reducing adjusted score equations in Kenne et al. (2017) <doi:10.1093/biomet/asx046>, or through the direct subtraction of an estimate of the bias of the maximum likelihood estimator from the maximum likelihood estimates as in Cordeiro and McCullagh (1991) <https://www.jstor.org/stable/2345592>. See Kosmidis et al (2020) <doi:10.1007/s11222-019-09860-6> for more details. Estimation in all cases takes place via a quasi Fisher scoring algorithm, and S3 methods for the construction of of confidence intervals for the reduced-bias estimates are provided. In the special case of generalized linear models for binomial and multinomial responses (both ordinal and nominal), the adjusted score approaches to mean and media bias reduction have been found to return estimates with improved frequentist properties, that are also always finite, even in cases where the maximum likelihood estimates are infinite (e.g. complete and quasi-complete separation; see Kosmidis and Firth, 2020 <doi:10.1093/biomet/asaa052>, for a proof for mean bias reduction in logistic regression). The 'mdyplFit()' fitting method fits logistic regression models using maximum Diaconis-Ylvisaker prior penalized likelihood, which also guarantees finite estimates. High-dimensionality corrections under proportional asymptotics can be applied to the resulting objects; see Sterzinger and Kosmidis (2024) <doi:10.48550/arXiv.2311.07419> for details.

Authors:Ioannis Kosmidis [aut, cre], Euloge Clovis Kenne Pagui [aut], Federico Boiocchi [ctb], Philipp Sterzinger [ctb], Kjell Konis [ctb], Nicola Sartori [ctb]

brglm2_1.1.0.tar.gz
brglm2_1.1.0.zip(r-4.7)brglm2_1.1.0.zip(r-4.6)brglm2_1.1.0.zip(r-4.5)
brglm2_1.1.0.tgz(r-4.6-x86_64)brglm2_1.1.0.tgz(r-4.6-arm64)brglm2_1.1.0.tgz(r-4.5-x86_64)brglm2_1.1.0.tgz(r-4.5-arm64)
brglm2_1.1.0.tar.gz(r-4.7-arm64)brglm2_1.1.0.tar.gz(r-4.7-x86_64)brglm2_1.1.0.tar.gz(r-4.6-arm64)brglm2_1.1.0.tar.gz(r-4.6-x86_64)
brglm2_1.1.0.tgz(r-4.6-emscripten)
manual.pdf |manual.html
DESCRIPTION |NEWS
card.svg |card.png
brglm2/json (API)

# Install 'brglm2' in R:
install.packages('brglm2', repos = c('https://ikosmidis.r-universe.dev', 'https://cloud.r-project.org'))

Bug tracker:https://github.com/ikosmidis/brglm2/issues

Datasets:
  • aids - Effects of AZT in slowing the development of AIDS symptoms
  • alligators - Alligator food choice data
  • coalition - Coalition data
  • endometrial - Histology grade and risk factors for 79 cases of endometrial cancer
  • enzymes - Liver enzyme data
  • hepatitis - Post-transfusion hepatitis: impact of non-A, non-B hepatitis surrogate tests
  • lizards - Habitat preferences of lizards
  • MultipleFeatures - Multiple features data
  • stemcell - Opinion on stem cell research and religious fundamentalism

On CRAN:

Conda:

adjusted-score-equationsalgorithmsbias-reducing-adjustmentsbias-reductionestimationglmhigh-dimensional-inferencelogistic-regressionmultinomial-regressionnegative-binomial-regressionnominal-responsesodds-ratioordinal-regressionordinal-responsesregressionregression-algorithmsrelative-riskstatistics

12.50 score 35 stars 15 packages 274 scripts 12k downloads 15 mentions 21 exports 8 dependencies

Last updated from:aac8a817d2. Checks:13 OK. Indexed: yes.

TargetResultTimeFilesSyslog
linux-devel-arm64OK159
linux-devel-x86_64OK178
source / vignettesOK200
linux-release-arm64OK168
linux-release-x86_64OK167
macos-release-arm64OK94
macos-release-x86_64OK223
macos-oldrel-arm64OK90
macos-oldrel-x86_64OK218
windows-develOK130
windows-releaseOK109
windows-oldrelOK125
wasm-releaseOK113

Exports:braclbrglm_controlbrglm_fitbrglmControlbrglmFitbrmultinombrnbcheck_infinite_estimatesdetect_separationexpomdypl_controlmdypl_fitmdyplControlmdyplFitmisordinal_superiorityplrtestse0se1sloesolve_se

Dependencies:enrichwithlatticeMASSMatrixnleqslvnnetnumDerivstatmod

Bias reduction in generalized linear models
The brglm2 package | This vignette | Other resources | Generalized linear models | Model | Score functions and information matrix | The expected information matrix about $\beta$ and $\phi$ is$$i =\left[\begin{array}{cc}i_{\beta\beta} & 0_p \0_p^\top & i_{\phi\phi}\end{array}\right] | Maximum likelihood estimation | Mean bias-reducing adjusted score functions | Median bias-reducing adjusted score functions | Mixed adjustments | Maximum penalized likelihood with powers of Jeffreys prior as penalty | Fitting algorithm in brglmFit | Input | Output | Iteration | Notes | Contributions to this vignette | Citation | References

Last update: 2025-08-29
Started: 2017-03-21

Adjacent category logit models using brglm2
bracl | Citation | Opinion on stem cell research and religious fundamentalism | Maximum likelihood estimation | Mean and median bias reduction | Prediction | Relevant resources | References

Last update: 2024-09-12
Started: 2019-02-06

Multinomial logistic regression using brglm2
brmultinom | Alligator data | Maximum likelihood estimation | Mean and median bias reduction | Infinite estimates and multinomial logistic regression | Relevant resources | Citation | References

Last update: 2024-09-12
Started: 2017-07-01

Negative binomial regression using brglm2
brnb | Ames salmonella data | Maximum likelihood estimation | Bias reduction | Asymptotic mean-bias correction | Mean-bias reducing adjusted score equations | Median-bias reducing adjusted score equations | Mixed bias reducing adjusted score equations | Relevant resources | Citation | References

Last update: 2024-09-12
Started: 2021-07-18

Estimating the exponential of regression parameters using brglm2
The expo() method | AIDS and AZT use | Infinite odds ratio estimates | brglmFit objects | References

Last update: 2023-02-07
Started: 2023-02-07

Readme and manuals

Help Manual

Help pageTopics
Effects of AZT in slowing the development of AIDS symptomsaids
Alligator food choice dataalligators
Bias reduction for adjacent category logit models for ordinal responses using the Poisson trick.bracl
brglm2: Bias Reduction in Generalized Linear Modelsbrglm2-package brglm2
Defunct Functions in package 'brglm2'brglm2-defunct check_infinite_estimates detect_separation
Auxiliary function for 'glm()' fitting using the 'brglmFit()' method.brglmControl brglm_control
Fitting function for 'glm()' for reduced-bias estimation and inferencebrglmFit brglm_fit
Bias reduction for multinomial response models using the Poisson trick.brmultinom
Bias reduction for negative binomial regression modelsbrnb
Coalition datacoalition
Extract model coefficients from '"brglmFit"' objectscoef.brglmFit
Extract estimates from '"brglmFit_expo"' objectscoef.brglmFit_expo
Extract model coefficients from '"brnb"' objectscoef.brnb
Method for computing confidence intervals for one or more regression parameters in a '"brglmFit"' objectconfint.brglmFit
Method for computing confidence intervals for one or more regression parameters in a '"brmultinom"' objectconfint.brmultinom
Method for computing Wald confidence intervals for one or more regression parameters in a '"brnb"' objectconfint.brnb
Method for computing confidence intervals for one or more regression parameters in a '"mdyplFit"' objectconfint.mdyplFit
Histology grade and risk factors for 79 cases of endometrial cancerendometrial
Liver enzyme dataenzymes
Estimate the exponential of parameters of generalized linear models using various methodsbrglmFit_expo expo expo.brglmFit expo.glm
Post-transfusion hepatitis: impact of non-A, non-B hepatitis surrogate testshepatitis
Habitat preferences of lizardslizards
Auxiliary function for 'glm()' fitting using the 'brglmFit()' method.mdyplControl mdypl_control
Fitting function for 'glm()' for maximum Diaconis-Ylvisaker prior penalized likelihood estimation of logistic regression modelsmdyplFit mdypl_fit
A '"link-glm"' object for misclassified responses in binomial regression modelsmis
Multiple features dataMultipleFeatures
Ordinal superiority scores of Agresti and Kateri (2017)ordinal_superiority ordinal_superiority.bracl
Penalized likelihood ratio test for '"mdyplFit"' objectsplrtest plrtest.mdyplFit
Predict method for bracl fitspredict.bracl
Predict method for brmultinom fitspredict.brmultinom
Residuals for multinomial logistic regression and adjacent category logit modelsresiduals.bracl residuals.brmultinom
MDYPL state evolution functions with no interceptse0
Logistic ridge regression state evolution functions with no interceptse0_ridge
MDYPL state evolution functions with interceptse1
Method for simulating a data set from '"brmultinom"' and '"bracl"' objectssimulate.brmultinom
Simulate Responsessimulate.brnb
Estimate the corrupted signal strength in a model with (sub-)Gaussian covariatessloe
Solve the MDYPL state evolution equations with or without intercept, with signal strength or contaminated signal strengthsolve_se
Opinion on stem cell research and religious fundamentalismstemcell
'summary()' method for '"brglmFit"' objectsprint.summary.brglmFit summary.brglmFit
'summary()' method for '"brnb"' objectsprint.summary.brnb summary.brnb
Summary method for '"mdyplFit"' objectsprint.summary.mdyplFit summary.mdyplFit
Return the variance-covariance matrix for the regression parameters in a 'brglmFit()' objectvcov.brglmFit
Extract model variance-covariance matrix from '"brnb"' objectsvcov.brnb