AER: Interfaces for AER package for data science pipelines.

Description Usage Arguments Details Value Author(s) Examples

Description

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

Usage

1
2

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
## Not run: 
library(intubate)
library(magrittr)
library(AER)


## ntbt_ivreg: Instrumental-Variable Regression
data("CigarettesSW", package = "AER")
CigarettesSW$rprice <- with(CigarettesSW, price/cpi)
CigarettesSW$rincome <- with(CigarettesSW, income/population/cpi)
CigarettesSW$tdiff <- with(CigarettesSW, (taxs - tax)/cpi)

## Original function to interface
ivreg(log(packs) ~ log(rprice) + log(rincome) | log(rincome) + tdiff + I(tax/cpi),
      data = CigarettesSW, subset = year == "1995")

## The interface puts data as first parameter
ntbt_ivreg(CigarettesSW,
           log(packs) ~ log(rprice) + log(rincome) | log(rincome) + tdiff + I(tax/cpi),
           subset = year == "1995")

## so it can be used easily in a pipeline.
CigarettesSW %>%
  ntbt_ivreg(log(packs) ~ log(rprice) + log(rincome) | log(rincome) + tdiff + I(tax/cpi),
             subset = year == "1995")


## ntbt_tobit: Tobit Regression
data("Affairs")

## Original function to interface
tobit(affairs ~ age + yearsmarried + religiousness + occupation + rating,
      data = Affairs)

## The interface puts data as first parameter
ntbt_tobit(Affairs,
           affairs ~ age + yearsmarried + religiousness + occupation + rating)

## so it can be used easily in a pipeline.
Affairs %>%
  ntbt_tobit(affairs ~ age + yearsmarried + religiousness + occupation + rating)

## End(Not run)

Example output

Loading required package: car
Loading required package: carData
Loading required package: lmtest
Loading required package: zoo

Attaching package:zooThe following objects are masked frompackage:base:

    as.Date, as.Date.numeric

Loading required package: sandwich
Loading required package: survival

Call:
ivreg(formula = log(packs) ~ log(rprice) + log(rincome) | log(rincome) +     tdiff + I(tax/cpi), data = CigarettesSW, subset = year ==     "1995")

Coefficients:
 (Intercept)   log(rprice)  log(rincome)  
      9.8950       -1.2774        0.2804  


Call:
ivreg(formula = log(packs) ~ log(rprice) + log(rincome) | log(rincome) +     tdiff + I(tax/cpi), data = CigarettesSW, subset = year ==     "1995")

Coefficients:
 (Intercept)   log(rprice)  log(rincome)  
      9.8950       -1.2774        0.2804  


Call:
ivreg(formula = log(packs) ~ log(rprice) + log(rincome) | log(rincome) +     tdiff + I(tax/cpi), data = ., subset = year == "1995")

Coefficients:
 (Intercept)   log(rprice)  log(rincome)  
      9.8950       -1.2774        0.2804  


Call:
tobit(formula = affairs ~ age + yearsmarried + religiousness + 
    occupation + rating, data = Affairs)

Coefficients:
  (Intercept)            age   yearsmarried  religiousness     occupation  
       8.1742        -0.1793         0.5541        -1.6862         0.3261  
       rating  
      -2.2850  

Scale: 8.247 


Call:
tobit(formula = affairs ~ age + yearsmarried + religiousness + 
    occupation + rating, data = Affairs)

Coefficients:
  (Intercept)            age   yearsmarried  religiousness     occupation  
       8.1742        -0.1793         0.5541        -1.6862         0.3261  
       rating  
      -2.2850  

Scale: 8.247 


Call:
tobit(formula = affairs ~ age + yearsmarried + religiousness + 
    occupation + rating, data = .)

Coefficients:
  (Intercept)            age   yearsmarried  religiousness     occupation  
       8.1742        -0.1793         0.5541        -1.6862         0.3261  
       rating  
      -2.2850  

Scale: 8.247 

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