Я столкнулся с этим фрагментом кода R, который имитирует сильную инструментальную переменную Z в соотношении между X и Y.
U<- MASS::mvrnorm(1000, c(20, 15), matrix(c(1, 0.5, 0.5, 1), 2, 2))
Xprime <- U[, 1]
C <- U[, 2]
Z <- rnorm(2000)
X <- Xprime + Z
Y <- 1 + X + C + rnorm(1000, 0, 1)
Если сейчас я хотел бы изменить код, чтобы показать, как ошибочен инструмент влияет на отношение, я знаю, что я должен изменить Z так, чтобы он больше не был экзогенным. У меня вопрос, как я могу изменить код, чтобы сделать Z эндогенным?