View source: R/standardization.R

standardization | R Documentation |

'standardization' uses a standard `glm`

linear model to perform parametric standardization
by adjusting bias through including all confounders as covariates. The model will calculate during training both the risk difference
and the risk ratio. Both can be accessed from the model as well as estimates of the counterfactuals of treatment.

standardization( data, f = NA, family = gaussian(), simple = pkg.env$simple, n.boot = 50, ... )

`data` |
a data frame containing the variables in the model.
This should be the same data used in |

`f` |
(optional) an object of class "formula" that overrides the default parameter |

`family` |
the family to be used in the general linear model.
By default, this is set to |

`simple` |
a boolean indicator to build default formula with interactions. If true, interactions will be excluded. If false, interactions will be included. By default, simple is set to false. NOTE: if this is changed, the coefficient for treatment may not accurately represent the average causal effect. |

`n.boot` |
an integer value that indicates number of bootstrap iterations to calculate standard error. |

`...` |
additional arguments that may be passed to the underlying |

`standardization`

returns an object of `class "standardization"`

.

The functions `print`

, `summary`

, and `predict`

can be used to interact with
the underlying `glm`

model.

An object of class `"standardization"`

is a list containing the following:

`call` |
the matched call. |

`formula` |
the formula used in the model. |

`model` |
the underlying glm model. |

`ATE` |
a data frame containing estimates of the treatment effect of the observed, counterfactuals, and risk metrics. |

`ATE.summary` |
a data frame containing the ATE, SE, and 95% CI of the ATE. |

library(causaldata) data(nhefs) nhefs.nmv <- nhefs[which(!is.na(nhefs$wt82)), ] nhefs.nmv$qsmk <- as.factor(nhefs.nmv$qsmk) confounders <- c( "sex", "race", "age", "education", "smokeintensity", "smokeyrs", "exercise", "active", "wt71" ) init_params(wt82_71, qsmk, covariates = confounders, data = nhefs.nmv ) # model using all defaults model <- standardization(data = nhefs.nmv) print(model) summary(model) print(model$ATE.summary) print(model$ATE.summary$Estimate[[2]] - model$ATE.summary$Estimate[[3]]) # manually calculate risk difference

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.