data(ivDiag)
controls <- c('altitudine', 'escursione', 'costal', 'nearsea', 'population',
'pop2', 'gini_land', 'gini_income')
ltz_out <- ltz(data = gsz, Y = "totassoc_p", D = "libero_comune_allnord",
Z = "bishopcity", controls = controls, weights = "population",
prior = c(0.178, 0.137))
plot_ltz(ltz_out)
library(testthat)
test_that("Check local-to-zero adjustment", {
expect_equal(as.numeric(ltz_out$ltz[1]), 3.6088)
})
Run the code above in your browser using DataLab