m <- lava::lvm(y ~ x+z, r ~ x) |>
lava::distribution(~ r, value = lava::binomial.lvm()) |>
transform(y0~r+y, value = \(x) { x[x[,1]==0,2] <- NA; x[,2] })
d <- lava::sim(m,1e3,seed=1)
aipw(y0 ~ x, data=d)
Run the code above in your browser using DataLab