MASS::glm.nb(totalContacts~Country*Setting + hhnum + age)

Call:
MASS::glm.nb(formula = totalContacts ~ Country * Setting + q8_hhnum +
    q1_age, data = mergedDat, init.theta = 3.153674563, link = log)

Deviance Residuals:
    Min       1Q   Median       3Q      Max
-3.0063  -0.7961  -0.2385   0.4181   4.7016

Coefficients:
                                Estimate Std. Error z value           Pr(>|z|)
(Intercept)                    2.8096597  0.0260088 108.027            < 2e-16
CountryIndonesia               0.0297085  0.0237604   1.250           0.211177
CountryTaiwan                 -0.1656966  0.0441520  -3.753           0.000175
CountryThailand               -0.4521072  0.0281449 -16.064            < 2e-16
CountryVietnam                -0.6127563  0.0232558 -26.348            < 2e-16
SettingUrban                  -0.3945613  0.0279131 -14.135            < 2e-16
q8_hhnum                       0.0339148  0.0028337  11.969            < 2e-16
q1_age                        -0.0046803  0.0002523 -18.554            < 2e-16
CountryIndonesia:SettingUrban  0.2571554  0.0342107   7.517 0.0000000000000561
CountryTaiwan:SettingUrban     0.2725435  0.0508279   5.362 0.0000000822679310
CountryThailand:SettingUrban   0.3306912  0.0400571   8.255            < 2e-16
CountryVietnam:SettingUrban    0.5347755  0.0347029  15.410            < 2e-16

(Intercept)                   ***
CountryIndonesia
CountryTaiwan                 ***
CountryThailand               ***
CountryVietnam                ***
SettingUrban                  ***
q8_hhnum                      ***
q1_age                        ***
CountryIndonesia:SettingUrban ***
CountryTaiwan:SettingUrban    ***
CountryThailand:SettingUrban  ***
CountryVietnam:SettingUrban   ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for Negative Binomial(3.1537) family taken to be 1)

    Null deviance: 17204  on 14008  degrees of freedom
Residual deviance: 14650  on 13997  degrees of freedom
  (6 observations deleted due to missingness)
AIC: 93991

Number of Fisher Scoring iterations: 1


              Theta:  3.1537
          Std. Err.:  0.0471

 2 x log-likelihood:  -93964.6420

# A tibble: 10 x 3
# Groups:   Setting [?]
   Setting Country   `mean(totalContacts)`
   <fct>   <fct>                     <dbl>
 1 Rural   Cambodia                  17.3
 2 Rural   Indonesia                 17.8
 3 Rural   Taiwan                    13.7
 4 Rural   Thailand                  10.7
 5 Rural   Vietnam                    9.25
 6 Urban   Cambodia                  12.0
 7 Urban   Indonesia                 15.4
 8 Urban   Taiwan                    12.3
 9 Urban   Thailand                  10.0
10 Urban   Vietnam                   10.6