tests/testthat/test-parties.R

In pos5737/pos5737data: Data for POS 5737

test_that("parties dataset replicates Clark and Golder's (2006) results", {
  # load data
  data(parties)
  # regression model from their table 2, pooled analysis, whole sample
  fit <- lm(enep ~ log(average_magnitude)*eneg + upper_tier*eneg + en_pres*proximity, data = parties)
  # check the intercept
  expect_equivalent(round(coef(fit)["(Intercept)"], 2), 2.81)
  # check the number of observations
  expect_equivalent(nobs(fit), 555)
  # check R^2
  expect_equivalent(round(summary(fit)$r.squared, 2), 0.30)
})