set.seed(12345)
n <- 30
clean_x <- rnorm(n, mean = 10, sd = sqrt(7))
delta_x <- rnorm(n, mean = 0, sd = sqrt(3))
e <- rnorm(n, mean = 0, sd = sqrt(5))
y <- 20 + 10 * clean_x + e
dirty_x <- clean_x + delta_x
mydata <- data.frame(y = y, dirtyx = dirty_x)
result <- eive.cga.formula(
formula = y ~ dirtyx,
dirtyx.varname = "dirtyx",
data = mydata,
numdummies = 10
)
Run the code above in your browser using DataLab