data-raw/generate_example_data.R

#linear substantive model with quadratic covariate effect
set.seed(1234)
n <- 1000
x <- rnorm(n)
z <- rnorm(n)
#v is an auxiliary variable
v <- x+rnorm(n)
y <- 1+z+x+x^2+rnorm(n)

#make some x values missing
xobsxb <- (y-mean(y))/sd(y)
xobspr <- exp(xobsxb)/(1+exp(xobsxb))
x[runif(n)>xobspr] <- NA

ex_linquad <- data.frame(y,z,x,xsq=x^2,v)

usethis::use_data(ex_linquad, overwrite=TRUE)
jwb133/bootmi documentation built on Sept. 22, 2024, 10:14 p.m.