test_that('Test output object', {
  expect_s3_class(
    mars(formula = effect ~ 1 + (1|district/study),
         data = school, studyID = 'district',
         variance = 'var', varcov_type = 'multilevel',
         structure = 'multilevel',
         estimation_method = "MLE"),
    "mars")
})