brglm: Interfaces for brglm package for data science pipelines.

In intubate: Interface to Popular R Functions for Data Science Pipelines

Description

Interfaces to brglm functions that can be used in a pipeline implemented by magrittr.

Usage

ntbt_brglm(data, ...)

Arguments

data	data frame, tibble, list, ...
...	Other arguments passed to the corresponding interfaced function.

Details

Interfaces call their corresponding interfaced function.

Value

Object returned by interfaced function.

Author(s)

Roberto Bertolusso

Examples

## Not run: 
library(intubate)
library(magrittr)
library(brglm)

## ntbt_brglm: Bias reduction in Binomial-response GLMs
data(lizards)

## Original function to interface
brglm(cbind(grahami, opalinus) ~ height + diameter +
        light + time, family = binomial(logit), data = lizards,
      method = "brglm.fit")
## The interface puts data as first parameter
ntbt_brglm(lizards, cbind(grahami, opalinus) ~ height + diameter +
             light + time, family = binomial(logit),
           method = "brglm.fit")

## so it can be used easily in a pipeline.
lizards %>%
  ntbt_brglm(cbind(grahami, opalinus) ~ height + diameter +
               light + time, family = binomial(logit),
             method = "brglm.fit")

## End(Not run)

Example output

Loading required package: profileModel
'brglm' will gradually be superseded by 'brglm2' (https://cran.r-project.org/package=brglm2), which provides utilities for mean and median bias reduction for all GLMs and methods for the detection of infinite estimates in binomial-response models.

Call:  brglm(formula = cbind(grahami, opalinus) ~ height + diameter +      light + time, family = binomial(logit), data = lizards, method = "brglm.fit") 

Coefficients:
 (Intercept)   height>=5ft  diameter>2in    lightshady    timemidday  
      1.9018        1.1064       -0.7536       -0.8177        0.2280  
    timelate  
     -0.7273  

Degrees of Freedom: 22 Total (i.e. Null);  17 Residual
Deviance:	    14.2462 
Penalized Deviance: -4.0065 	AIC: 83.0704 

Call:  brglm(formula = cbind(grahami, opalinus) ~ height + diameter +      light + time, family = binomial(logit), data = lizards, method = "brglm.fit") 

Coefficients:
 (Intercept)   height>=5ft  diameter>2in    lightshady    timemidday  
      1.9018        1.1064       -0.7536       -0.8177        0.2280  
    timelate  
     -0.7273  

Degrees of Freedom: 22 Total (i.e. Null);  17 Residual
Deviance:	    14.2462 
Penalized Deviance: -4.0065 	AIC: 83.0704 

Call:  brglm(formula = cbind(grahami, opalinus) ~ height + diameter +      light + time, family = binomial(logit), data = ., method = "brglm.fit") 

Coefficients:
 (Intercept)   height>=5ft  diameter>2in    lightshady    timemidday  
      1.9018        1.1064       -0.7536       -0.8177        0.2280  
    timelate  
     -0.7273  

Degrees of Freedom: 22 Total (i.e. Null);  17 Residual
Deviance:	    14.2462 
Penalized Deviance: -4.0065 	AIC: 83.0704 

intubate documentation built on May 2, 2019, 2:46 p.m.