me: Marginal effects at unique combinations of independent...

View source: R/me.R

meR Documentation

Marginal effects at unique combinations of independent variables

Description

me computes the marginal effects of variable x for the distinct values of variables x and over.

Usage

me(
  x,
  over = NULL,
  model = NULL,
  data = NULL,
  formula = NULL,
  link = NULL,
  coefficients = NULL,
  vcov = NULL,
  discrete = FALSE,
  discrete_step = 1,
  at = NULL,
  mc = FALSE,
  pct = c(lb = 2.5, ub = 97.5),
  iter = 1000,
  weights = NULL
)

Arguments

x

a character string representing the name of the main variable of interest. Marginal effects will be computed for this variable.

over

a character string representing the name of the conditionning variable. DAME will be computed for the bins long the range of this variable.

model

fitted model object. The package works best with GLM objects and will extract the formula, dataset, family, coefficients, and the QR components of the design matrix if arguments formula, data, link, coefficients, and/or vcov are not explicitly specified.

data

the dataset to be used to compute marginal effects (if not specified, it is extracted from the fitted model object).

formula

the formula used in estimation (if not specified, it is extracted from the fitted model object).

link

the name of the link function used in estimation (if not specified, it is extracted from the fitted model object).

coefficients

the named vector of coefficients produced during the estimation (if not specified, it is extracted from the fitted model object).

vcov

the variance-covariance matrix to be used for computing standard errors (if not specified, it is extracted from the fitted model object).

discrete

A logical variable. If TRUE, the function will compute the effect of a discrete change in x. If FALSE, the function will compute the partial derivative of x.

discrete_step

The size of a discrete change in x used in computations (used only if discrete=TRUE).

at

an optional named list of values of independent variables. These variables will be set to these value before computations. The remaining numeric variables (except x and over) will be set to their means. The remaining factor variables will be set to their modes.

mc

logical. If TRUE, the standard errors and confidence intervals will be computed using simulations. If FALSE (default), the delta method will be used.

pct

a named numeric vector with the sampling quantiles to be output with the DAME estimates (the names are used as the new variable names). Default = c(lb=2.5,ub=97.5).

iter

the number of interations used in Monte-Carlo simulations. Default = 1,000.

weights

an optional vector of sampling weights.

Value

me returns a data frame with the estimates of the marginal effects for each combination of the variables specified in x and over, along with standard errors, confidence intervals, and the used values of the independent variables. All quantitative variable not included in at, x and over are set to their means, and all qualitative variables (except those listed in at, x and over) are converted to factors and set to their modes.

Examples

##poisson regression with 2 variables and interaction between them
#fit the regression first
data <- data.frame(y = rpois(10000, 10), x2 = rpois(10000, 5), x1 = rpois(10000, 3))
y <- glm(y ~ x1 + x2 + x1*x2, data = data, family = "poisson")
me(model = y, x = "x1", over = "x2")

andreizhirnov/DAME documentation built on Nov. 21, 2022, 5:55 p.m.