ml_logistic_loglinear: Maximum Likelihood with Logistic Regression Disease Model and...

In vandomed/meuc: Measurement Error and Unmeasured Confounding

Description

Calculates maximum likelihood estimates for measurement error/unmeasured confounding scenario where true disease model is logistic regression and measurement error model is log-transformed linear regression.

Usage

ml_logistic_loglinear(all_data = NULL, main = NULL, internal = NULL,
  external = NULL, y_var, z_var, d_vars = NULL, c_vars = NULL,
  b_vars = NULL, tdm_covariates = NULL, mem_covariates = NULL,
  approx_integral = TRUE, integrate_tol = 1e-08,
  integrate_tol_hessian = integrate_tol, estimate_var = TRUE, ...)

Arguments

all_data	Data frame with data for main study and validation study.
main	Data frame with data for the main study.
internal	Data frame with data for internal validation study.
external	Data frame with data for the external validation study.
y_var	Character string specifying name of Y variable.
z_var	Character string specifying name of Z variable.
d_vars	Character string specifying name of D variables.
c_vars	Character vector specifying names of C variables.
b_vars	Character vector specifying names of variables in true disease model but not in measurement error model.
tdm_covariates	Character vector specifying variables in true disease model. The Z variable is automatically included whether you include it in tdm_covariates or not.
mem_covariates	Character vector specifying variables in measurement error model.
approx_integral	Logical value for whether to use the probit approximation for the logistic-normal integral, to avoid numerically integrating Z's out of the likelihood function.
integrate_tol	Numeric value specifying tol input to hcubature for numerical integration.
integrate_tol_hessian	Same as integrate_tol, but for use when estimating the Hessian matrix only. Sometimes using a smaller value than for likelihood maximization helps prevent cases where the inverse Hessian is not positive definite.
estimate_var	Logical value for whether to return variance-covariance matrix for parameter estimates.
...	Additional arguments to pass to nlminb.

Details

The true disease model is:

logit[P(Y = 1)] = beta_0 + beta_z Z + beta_c^T C + beta_b^T B

The measurement error model is:

ln(Z) = alpha_0 + alpha_d^T D + alpha_c^T C + d, d ~ N(0, sigsq_d)

There should be main study data with (Y, D, C, B) as well as internal validation data with (Y, Z, D, C, B) and/or external validation data with (Z, D, C).

Value

List containing:

1. Numeric vector of parameter estimates.

2. Variance-covariance matrix (if estimate_var = TRUE).

3. Returned nlminb object from maximizing the log-likelihood function.

4. Akaike information criterion (AIC).

vandomed/meuc documentation built on May 12, 2019, 6:17 p.m.