divorce_df <- bacondecomp::divorce
stevenson_wolfson_replicates <- c(
1.6,
-1.5,
-1.5,
-3.0,
-8.0,
-10.0,
-11.9,
-12.8,
-13.3,
-16.4,
-18.7
)
model_fit <- lm(data = divorce_df[divorce_df$sex == 2, ],
suiciderate_elast_jag ~ factor(chyrspos) + factor(st) + factor(year))
model_coefs <- (model_fit$coefficients)[2:12]
names(model_coefs) <- NULL
model_coefs <- round(model_coefs, 3)
model_coefs <- model_coefs * 100
test_that("Table 1 col 1 Stevenson and Wolfers Replicates",
{
expect_equal(
model_coefs,
stevenson_wolfson_replicates
)
})
## TODO We don't actually test anything related to bacon here...
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.