inst/doc/compare-smm-fits.R

In OneSampleMR: One Sample Mendelian Randomization and Instrumental Variable Analyses

## ----include=FALSE------------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(OneSampleMR)

## -----------------------------------------------------------------------------
set.seed(12345)
n    <- 5000
psi0 <- 0.5
psi1 <- 0.2
Z    <- rbinom(n, 1, 0.5)
X    <- rbinom(n, 1, 0.7*Z + 0.2*(1 - Z)) 
m0   <- plogis(1 + 0.8*X - 0.39*Z)
Y    <- rbinom(n, 1, plogis(psi0*X + log(m0/(1 - m0)))) 
dat  <- data.frame(Z, X, Y)

## -----------------------------------------------------------------------------
fit02 <- msmm(Y ~ X | Z, data = dat)
summary(fit02)

## -----------------------------------------------------------------------------
fit03 <- msmm(Y ~ X | Z, data = dat, estmethod = "gmmalt")
summary(fit03)

## -----------------------------------------------------------------------------
fit04 <- msmm(Y ~ X | Z, data = dat, estmethod = "tsls")
summary(fit04)

## -----------------------------------------------------------------------------
fit05 <- msmm(Y ~ X | Z, data = dat, estmethod = "tslsalt")
summary(fit05)