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
card.svg |card.png
brglm2/json (API)
NEWS

# 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.42 score 34 stars 14 packages 222 scripts 14k downloads 15 mentions 21 exports 8 dependencies

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

TargetResultTimeFilesSyslog
linux-devel-arm64OK156
linux-devel-x86_64OK155
source / vignettesOK211
linux-release-arm64OK151
linux-release-x86_64OK128
macos-release-arm64OK129
macos-release-x86_64OK226
macos-oldrel-arm64OK106
macos-oldrel-x86_64OK201
windows-develOK170
windows-releaseOK114
windows-oldrelOK116
wasm-releaseOK111

Exports:braclbrglm_controlbrglm_fitbrglmControlbrglmFitbrmultinombrnbcheck_infinite_estimatesdetect_separationexpomdypl_controlmdypl_fitmdyplControlmdyplFitmisordinal_superiorityplrtestse0se1sloesolve_se

Dependencies:enrichwithlatticeMASSMatrixnleqslvnnetnumDerivstatmod

Adjacent category logit models using brglm2

Rendered fromadjacent.Rmdusingknitr::rmarkdownon Jun 05 2026.

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

Bias reduction in generalized linear models

Rendered fromiteration.Rmdusingknitr::rmarkdownon Jun 05 2026.

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

Estimating the exponential of regression parameters using brglm2

Rendered fromexpo.Rmdusingknitr::rmarkdownon Jun 05 2026.

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

Multinomial logistic regression using brglm2

Rendered frommultinomial.Rmdusingknitr::rmarkdownon Jun 05 2026.

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

Negative binomial regression using brglm2

Rendered fromnegativeBinomial.Rmdusingknitr::rmarkdownon Jun 05 2026.

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

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