packages <- c(“tidyverse”, “infer”, “fst”, “modelsummary”, “effects”, “survey”, “performance”, “flextable”, “broom”, “scales”, “ggeffects”, “marginaleffects”) # add any you need here
new_packages <- packages[!(packages %in% installed.packages()[,“Package”])] if(length(new_packages)) install.packages(new_packages)
lapply(packages, library, character.only = TRUE) install.packages(“fst”)
library(fst)
library(modelsummary)
## `modelsummary` 2.0.0 now uses `tinytable` as its default table-drawing
## backend. Learn more at: https://vincentarelbundock.github.io/tinytable/
##
## Revert to `kableExtra` for one session:
##
## options(modelsummary_factory_default = 'kableExtra')
##
## Change the default backend persistently:
##
## config_modelsummary(factory_default = 'gt')
##
## Silence this message forever:
##
## config_modelsummary(startup_message = FALSE)
library(infer)
library(fst)
file_path <- "/Users/Jessica/Downloads/norway_data.fst"
Norway_data <- read_fst("norway_data.fst")
# Null hypothesis (H0): There is no relationship between income inequality and social trust in Norway.
N0 <- "There is no significant relationship between income inequality and social trust in Norway."
# Alternative hypothesis (H1): Higher income inequality is associated with lower levels of social trust in Norway.
H1 <- "Higher income inequality is associated with lower levels of social trust in Norway."
# Interaction hypothesis (H2): The relationship between income inequality and social trust is moderated by social cohesion.
H2 <- "The relationship between income inequality and social trust in Norway is moderated by social cohesion"
library(modelsummary)
summary_table <- datasummary_skim(Norway_data)
## Warning: These variables were omitted because they include more than 50 levels:
## cntbrthb, cntbrthd, fbrncnt, fbrncnta, fbrncntc, lnghom1, mbrncnt, mbrncnta,
## mbrncntc.
## Warning: These variables were omitted because they are entirely missing:
## prtvde1, prtvade1, prtvbde1, prtvcde1, prtvdde1, prtvede1, prtvfde1, prtvde2,
## prtvade2, prtvbde2, prtvcde2, prtvdde2, prtvede2, prtvfde2, prtvlt1, prtvalt1,
## prtvblt1, prtvclt1, prtvlt2, prtvalt2, prtvblt2, prtvclt2, prtvlt3, prtvalt3,
## prtvblt3, prtvclt3, prtvtal, prtvtat, prtvtaat, prtvtbat, prtvtcat, prtvtbe,
## prtvtabe, prtvtbbe, prtvtcbe, prtvtdbe, prtvtebe, prtvtbg, prtvtabg, prtvtbbg,
## prtvtcbg, prtvtdbg, prtvtebg, prtvtch, prtvtach, prtvtbch, prtvtcch, prtvtdch,
## prtvtech, prtvtfch, prtvtgch, prtvthch, prtvtcy, prtvtacy, prtvtbcy, prtvtccy,
## prtvtcz, prtvtacz, prtvtbcz, prtvtccz, prtvtdcz, prtvtecz, prtvtdk, prtvtadk,
## prtvtbdk, prtvtcdk, prtvtddk, prtvtee, prtvtaee, prtvtbee, prtvtcee, prtvtdee,
## prtvteee, prtvtfee, prtvtgee, prtvthee, prtvtes, prtvtaes, prtvtbes, prtvtces,
## prtvtdes, prtvtees, prtvtfes, prtvtfi, prtvtafi, prtvtbfi, prtvtcfi, prtvtdfi,
## prtvtefi, prtvtfr, prtvtafr, prtvtbfr, prtvtcfr, prtvtdfr, prtvtefr, prtvtgb,
## prtvtagb, prtvtbgb, prtvtcgb, prtvtdgb, prtvtgr, prtvtagr, prtvtbgr, prtvtcgr,
## prtvtdgr, prtvthr, prtvtahr, prtvtbhr, prtvthu, prtvtahu, prtvtbhu, prtvtchu,
## prtvtdhu, prtvtehu, prtvtfhu, prtvtghu, prtvtie, prtvtaie, prtvtbie, prtvtcie,
## prtvtdie, prtvtil, prtvtail, prtvtbil, prtvtcil, prtvtdil, prtvtis, prtvtais,
## prtvtbis, prtvtcis, prtvtdis, prtvtit, prtvtait, prtvtbit, prtvtcit, prtvtdit,
## prtvtlu, prtvtlv, prtvtalv, prtvtnl, prtvtanl, prtvtbnl, prtvtcnl, prtvtdnl,
## prtvtenl, prtvtfnl, prtvtgnl, prtvthnl, prtvtpl, prtvtapl, prtvtbpl, prtvtcpl,
## prtvtdpl, prtvtepl, prtvtpt, prtvtapt, prtvtbpt, prtvtcpt, prtvtdpt, prtvtro,
## prtvtaro, prtvtru, prtvtaru, prtvtbru, prtvtcru, prtvtdru, prtvtse, prtvtase,
## prtvtbse, prtvtcse, prtvtdse, prtvtsi, prtvtasi, prtvtbsi, prtvtcsi, prtvtdsi,
## prtvtesi, prtvtfsi, prtvtsk, prtvtask, prtvtbsk, prtvtcsk, prtvtdsk, prtvtesk,
## prtvttr, prtvtatr, prtvtua, prtvtaua, prtvtbua, prtvtcua, prtvtxk, prtvtrs,
## prtvtars, prtvtme, prtvtame, prtclal, prtclat, prtclaat, prtclcat, prtcldat,
## prtclbe, prtclabe, prtclbbe, prtclcbe, prtcldbe, prtclebe, prtclbg, prtclabg,
## prtclbbg, prtclcbg, prtcldbg, prtclebg, prtclch, prtclach, prtclbch, prtclcch,
## prtcldch, prtclech, prtclfch, prtclgch, prtclhch, prtclsk, prtclask, prtclbsk,
## prtclcsk, prtcldsk, prtclesk, prtclcy, prtclacy, prtclbcy, prtclccy, prtclcz,
## prtclacz, prtclbcz, prtclccz, prtcldcz, prtclecz, prtclde, prtclade, prtclbde,
## prtclcde, prtcldde, prtclede, prtclfde, prtcldk, prtcladk, prtclbdk, prtclcdk,
## prtclddk, prtclee, prtclaee, prtclbee, prtclcee, prtcldee, prtcleee, prtclfee,
## prtclgee, prtclhee, prtcles, prtclaes, prtclbes, prtclces, prtcldes, prtclees,
## prtclfes, prtclges, prtclfi, prtclafi, prtclbfi, prtclcfi, prtcldfi, prtclefi,
## prtclffi, prtclfr, prtclafr, prtclbfr, prtclcfr, prtcldfr, prtclefr, prtclffr,
## prtclgb, prtclagb, prtclbgb, prtclcgb, prtcldgb, prtclgr, prtclagr, prtclbgr,
## prtclcgr, prtcldgr, prtclhr, prtclahr, prtclbhr, prtclhu, prtclahu, prtclbhu,
## prtclchu, prtcldhu, prtclehu, prtclfhu, prtclghu, prtclhhu, prtclie, prtclaie,
## prtclbie, prtclcie, prtcldie, prtcleie, prtclfie, prtclil, prtclail, prtclbil,
## prtclcil, prtcldil, prtcleil, prtclis, prtclais, prtclbis, prtclcis, prtcldis,
## prtclit, prtclait, prtclbit, prtclcit, prtcldit, prtcleit, prtcllt, prtclalt,
## prtclblt, prtclclt, prtcllu, prtcllv, prtclalv, prtclnl, prtclanl, prtclbnl,
## prtclcnl, prtcldnl, prtclenl, prtclfnl, prtclgnl, prtclpl, prtclapl, prtclbpl,
## prtclcpl, prtcldpl, prtclepl, prtclfpl, prtclgpl, prtclhpl, prtclipl, prtclpt,
## prtclapt, prtclbpt, prtclcpt, prtcldpt, prtclept, prtclfpt, prtclro, prtclaro,
## prtclru, prtclaru, prtclbru, prtclcru, prtcldru, prtclse, prtclase, prtclbse,
## prtclcse, prtcldse, prtclsi, prtclasi, prtclbsi, prtclcsi, prtcldsi, prtclesi,
## prtclfsi, prtcltr, prtclatr, prtclua, prtclaua, prtclbua, prtclcua, prtcldua,
## prtclxk, prtclrs, prtclars, prtclme, prtclame, prtmbat, prtmbaat, prtmbbe,
## prtmbabe, prtmbbbe, prtmbcbe, prtmbbg, prtmbabg, prtmbbbg, prtmbch, prtmbach,
## prtmbbch, prtmbcch, prtmbcy, prtmbcz, prtmbacz, prtmbbcz, prtmbde, prtmbade,
## prtmbbde, prtmbcde, prtmbdk, prtmbadk, prtmbbdk, prtmbee, prtmbaee, prtmbbee,
## prtmbcee, prtmbes, prtmbaes, prtmbbes, prtmbfi, prtmbafi, prtmbbfi, prtmbfr,
## prtmbafr, prtmbbfr, prtmbcfr, prtmbgb, prtmbagb, prtmbgr, prtmbagr, prtmbbgr,
## prtmbcgr, prtmbhr, prtmbhu, prtmbahu, prtmbbhu, prtmbchu, prtmbie, prtmbaie,
## prtmbbie, prtmbil, prtmbail, prtmbbil, prtmbis, prtmbit, prtmbait, prtmblt,
## prtmblu, prtmblv, prtmbnl, prtmbanl, prtmbbnl, prtmbcnl, prtmbpl, prtmbapl,
## prtmbbpl, prtmbcpl, prtmbdpl, prtmbpt, prtmbapt, prtmbbpt, prtmbro, prtmbaro,
## prtmbru, prtmbaru, prtmbbru, prtmbse, prtmbase, prtmbsi, prtmbasi, prtmbbsi,
## prtmbcsi, prtmbsk, prtmbask, prtmbbsk, prtmbtr, prtmbatr, prtmbua, prtmbaua,
## prtmbbua, prtmbcua, vteurmmb, vteumbgb, rlgdnal, rlgdnat, rlgdnbat, rlgdnbe,
## rlgdnch, rlgdnach, rlgdncy, rlgdnacy, rlgdnde, rlgdnade, rlgdndk, rlgdnfi,
## rlgdnafi, rlgdngb, rlgdngr, rlgdnagr, rlgdnhu, rlgdnie, rlgdnil, rlgdnis,
## rlgdnais, rlgdnlt, rlgdnlv, rlgdnnl, rlgdnanl, rlgdnpl, rlgdnapl, rlgdnbpl,
## rlgdnpt, rlgdnro, rlgdnru, rlgdnaru, rlgdnse, rlgdnase, rlgdnsi, rlgdnsk,
## rlgdnask, rlgdnbsk, rlgdnua, rlgdnrs, rlgdme, rlgdeal, rlgdeat, rlgdebat,
## rlgdebe, rlgdech, rlgdeach, rlgdecy, rlgdeacy, rlgdede, rlgdeade, rlgdedk,
## rlgdefi, rlgdeafi, rlgdegb, rlgdegr, rlgdeagr, rlgdehu, rlgdeie, rlgdeil,
## rlgdeis, rlgdeais, rlgdelt, rlgdelv, rlgdenl, rlgdeanl, rlgdepl, rlgdeapl,
## rlgdebpl, rlgdept, rlgdero, rlgderu, rlgdearu, rlgdese, rlgdease, rlgdesi,
## rlgdesk, rlgdeask, rlgdebsk, rlgdeua, rlgders, rlgdeme, gndr16, gndr17, gndr18,
## gndr19, gndr20, gndr21, gndr22, gndr23, gndr24, rshipa10, rshipa11, rshipa12,
## rshipa13, rshipa14, rshipa15, rshipa16, rshipa17, rshipa18, rshipa19, rshipa20,
## rshipa21, rshipa22, rshipa23, rshipa24, yrbrn16, yrbrn17, yrbrn18, yrbrn19,
## yrbrn20, yrbrn21, yrbrn22, yrbrn23, yrbrn24, occf14ie, occm14ie, wrkctrhu,
## edagegb, edlvat, edlveat, eduat1, eduat2, edlvbe, edlvabe, edlvdbe, edlvebe,
## edlvbg, edlvdbg, edlvebg, edlvch, edlvach, edlvbch, edlvcch, edlvdch, edlvcy,
## edlvacy, edlvdcy, edlvecy, edlvgcy, edlvcz, edlvdcz, edlvdal, edlvde, edlvade,
## edlvdfi, edlvdis, edlvdk, edlvadk, edlvddk, edlvdxk, edlvee, edlvaee, edlvbee,
## edlvdee, edlves, edlvaes, edlvdes, edlvees, edlvges, edlvhes, edlvfr, edlvafr,
## edlvbfr, edlvdfr, edlvgb, edlvagb, edlvgr, edlvagr, edlvdgr, edlvegr, edlvhr,
## edlvdhr, edlvehr, edlvhu, edlvahu, edlvbhu, edlvdhu, edlvdahu, edlvie, edlvaie,
## edlvbie, edlvdie, edlvil, edlvail, eduil1, eduail1, edubil1, edlvit, edlvait,
## edlvdit, edlveit, edlvlt, edlvdlt, edlvlu, edlvlv, edlvdlv, edlvelv, edlvnl,
## edlvdnl, edlvenl, edlvpl, edlvapl, edlvbpl, edlvdpl, edlvepl, edlvgpl, edlvhpl,
## edlvpt, edlvapt, edlvbpt, edlvdpt, edlvro, edlvru, edlvdru, edlvse, edlvase,
## edlvdse, edlvsi, edlvasi, edlvdsi, edlvesi, edlvsk, edlvask, edlvdsk, edlvtr,
## edlvua, edlvaua, edlvdua, edude1, eduade1, edubde1, educde1, edude2, eduade2,
## edubde2, edude3, eduade3, edugb2, eduagb2, edubgb2, eduil2, eduail2, edupl2,
## eduyrpl, edugb1, eduagb1, edubgb1, educgb1, edlvdrs, edlvdme, edagepgb,
## edlvpat, edlvpeat, edlvpfat, edlvpbe, edlvpdbe, edlvpebe, edlvpch, edlvpdch,
## edlvpcy, edlvpdcy, edlvpecy, edlvpgcy, edlvpcz, edlvpdcz, edlvpdbg, edlvpebg,
## edlvpdfi, edlvpdis, edlvpdit, edlvpeit, edlvpdk, edlvpddk, edlvpdxk, edlvpee,
## edlvpdee, edlvpes, edlvpdes, edlvpees, edlvpfes, edlvpges, edlvpfr, edlvpdfr,
## edlvpgb, edlvpgr, edlvpdgr, edlvpegr, edlvphr, edlvpdhr, edlvpehr, edlvphu,
## edlvpdhu, edlvpdahu, edlvpie, edlvpdie, edlvpil, edupil1, edupail1, edupbil1,
## edlvplt, edlvpdlt, edlvplv, edlvpdlv, edlvpelv, edlvpnl, edlvpdnl, edlvpenl,
## edlvppl, edlvpdpl, edlvpepl, edlvpfpl, edlvpgpl, edlvppt, edlvpdpt, edlvpro,
## edlvpru, edlvpdru, edlvpse, edlvpdse, edlvpsi, edlvpdsi, edlvpesi, edlvpsk,
## edlvpdsk, edlvptr, edlvpua, edlvpdua, edupde1, edupade1, edupbde1, edupcde1,
## edupde2, edupade2, edupbde2, edupde3, edupade3, edupgb2, edupagb2, edupbgb2,
## edupil2, edupail2, eduppl2, eduyrppl, edlvpdal, edupgb1, edupagb1, edupbgb1,
## edupcgb1, edlvpdrs, edlvpdme, edlvfdal, edlvfat, edlvfeat, edlvfbe, edlvfdbe,
## edlvfebe, edlvfdbg, edlvfebg, edlvfhr, edlvfdhr, edlvfehr, edlvfcy, edlvfdcy,
## edlvfecy, edlvfgcy, edlvfcz, edlvfdcz, edlvfdk, e
print(summary_table)
##
## +----------+------------+-------+-------+
## | | | N | % |
## +==========+============+=======+=======+
## | name | ESS10e03_1 | 1411 | 8.8 |
## +----------+------------+-------+-------+
## | | ESS1e06_6 | 2036 | 12.7 |
## +----------+------------+-------+-------+
## | | ESS2e03_6 | 1760 | 11.0 |
## +----------+------------+-------+-------+
## | | ESS3e03_7 | 1750 | 10.9 |
## +----------+------------+-------+-------+
## | | ESS4e04_5 | 1549 | 9.6 |
## +----------+------------+-------+-------+
## | | ESS5e03_4 | 1548 | 9.6 |
## +----------+------------+-------+-------+
## | | ESS6e02_5 | 1624 | 10.1 |
## +----------+------------+-------+-------+
## | | ESS7e02_2 | 1436 | 8.9 |
## +----------+------------+-------+-------+
## | | ESS8e02_2 | 1545 | 9.6 |
## +----------+------------+-------+-------+
## | | ESS9e03_1 | 1406 | 8.8 |
## +----------+------------+-------+-------+
## | proddate | 01.12.2018 | 10079 | 62.7 |
## +----------+------------+-------+-------+
## | | 10.12.2020 | 1545 | 9.6 |
## +----------+------------+-------+-------+
## | | 17.02.2021 | 1406 | 8.8 |
## +----------+------------+-------+-------+
## | | 19.01.2023 | 1624 | 10.1 |
## +----------+------------+-------+-------+
## | | 30.06.2023 | 1411 | 8.8 |
## +----------+------------+-------+-------+
## | | NO | 16065 | 100.0 |
## +----------+------------+-------+-------+
## | cntbrth | 66 | 1903 | 11.8 |
## +----------+------------+-------+-------+
## | | AT | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | BE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | CZ | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | DK | 17 | 0.1 |
## +----------+------------+-------+-------+
## | | ES | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 13 | 0.1 |
## +----------+------------+-------+-------+
## | | GH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | ID | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | IR | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | KE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | KR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | KZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LB | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | MA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | RO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | RW | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 14 | 0.1 |
## +----------+------------+-------+-------+
## | | SG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | VN | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | YU | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 14029 | 87.3 |
## +----------+------------+-------+-------+
## | cntbrtha | 66 | 3258 | 20.3 |
## +----------+------------+-------+-------+
## | | 99 | 58 | 0.4 |
## +----------+------------+-------+-------+
## | | AF | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | AL | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | BE | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | CG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | CN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CS | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | DK | 17 | 0.1 |
## +----------+------------+-------+-------+
## | | ER | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ES | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ET | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | HR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HU | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | ID | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | IR | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | KZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | MT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | MX | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | PE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 21 | 0.1 |
## +----------+------------+-------+-------+
## | | SO | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | UA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 18 | 0.1 |
## +----------+------------+-------+-------+
## | | UY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | VN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 12555 | 78.2 |
## +----------+------------+-------+-------+
## | cntbrthc | 66 | 4063 | 25.3 |
## +----------+------------+-------+-------+
## | | 88 | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | 99 | 283 | 1.8 |
## +----------+------------+-------+-------+
## | | BA | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | BY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 23 | 0.1 |
## +----------+------------+-------+-------+
## | | DK | 19 | 0.1 |
## +----------+------------+-------+-------+
## | | ER | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | ET | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 9 | 0.1 |
## +----------+------------+-------+-------+
## | | GR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | IR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IT | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | KR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | LT | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | LV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 45 | 0.3 |
## +----------+------------+-------+-------+
## | | RO | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 49 | 0.3 |
## +----------+------------+-------+-------+
## | | SG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SO | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | VN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | YE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ZM | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 11460 | 71.3 |
## +----------+------------+-------+-------+
## | ctzship | 66 | 1975 | 12.3 |
## +----------+------------+-------+-------+
## | | AT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DK | 13 | 0.1 |
## +----------+------------+-------+-------+
## | | ES | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | GH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | MA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PS | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RW | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | SG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 14029 | 87.3 |
## +----------+------------+-------+-------+
## | ctzshipa | 66 | 3373 | 21.0 |
## +----------+------------+-------+-------+
## | | 99 | 35 | 0.2 |
## +----------+------------+-------+-------+
## | | AL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | BE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CH | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | DK | 13 | 0.1 |
## +----------+------------+-------+-------+
## | | ET | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 14 | 0.1 |
## +----------+------------+-------+-------+
## | | SO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | UA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | NA | 12555 | 78.2 |
## +----------+------------+-------+-------+
## | ctzshipb | 66 | 2919 | 18.2 |
## +----------+------------+-------+-------+
## | | 88 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | 99 | 34 | 0.2 |
## +----------+------------+-------+-------+
## | | BA | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | BE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CD | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | DK | 17 | 0.1 |
## +----------+------------+-------+-------+
## | | ER | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | HU | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | IT | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LV | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | MA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 21 | 0.1 |
## +----------+------------+-------+-------+
## | | RS | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 25 | 0.2 |
## +----------+------------+-------+-------+
## | | SO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 12968 | 80.7 |
## +----------+------------+-------+-------+
## | ctzshipc | 66 | 4282 | 26.7 |
## +----------+------------+-------+-------+
## | | 99 | 176 | 1.1 |
## +----------+------------+-------+-------+
## | | DE | 18 | 0.1 |
## +----------+------------+-------+-------+
## | | DK | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | ER | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | IR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | IT | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | LT | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | LV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 40 | 0.2 |
## +----------+------------+-------+-------+
## | | RO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 32 | 0.2 |
## +----------+------------+-------+-------+
## | | US | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 11460 | 71.3 |
## +----------+------------+-------+-------+
## | ctzshipd | 6666 | 1337 | 8.3 |
## +----------+------------+-------+-------+
## | | 8888 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BG | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | DK | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | ER | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | ET | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | JP | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | NO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NP | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PL | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | PT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SD | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 9 | 0.1 |
## +----------+------------+-------+-------+
## | | SO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SY | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 14659 | 91.2 |
## +----------+------------+-------+-------+
## | fbrncntb | 66 | 3983 | 24.8 |
## +----------+------------+-------+-------+
## | | 88 | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | 99 | 302 | 1.9 |
## +----------+------------+-------+-------+
## | | BA | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | BY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | CU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 30 | 0.2 |
## +----------+------------+-------+-------+
## | | DK | 32 | 0.2 |
## +----------+------------+-------+-------+
## | | DZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ER | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 17 | 0.1 |
## +----------+------------+-------+-------+
## | | GR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | HU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 9 | 0.1 |
## +----------+------------+-------+-------+
## | | IR | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | IT | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | LT | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | MK | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | MY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | NO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | PL | 45 | 0.3 |
## +----------+------------+-------+-------+
## | | PT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RO | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | RS | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 54 | 0.3 |
## +----------+------------+-------+-------+
## | | SO | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | VN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | ZA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 11460 | 71.3 |
## +----------+------------+-------+-------+
## | lnghoma | 777 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | 888 | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | 999 | 29 | 0.2 |
## +----------+------------+-------+-------+
## | | ALB | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | AMH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ARA | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | BAL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BOS | 13 | 0.1 |
## +----------+------------+-------+-------+
## | | CHE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CHI | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DAN | 14 | 0.1 |
## +----------+------------+-------+-------+
## | | DUT | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | ENG | 46 | 0.3 |
## +----------+------------+-------+-------+
## | | FIN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | FRE | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | GER | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | HIN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HUN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ICE | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | IND | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IRA | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | KUR | 15 | 0.1 |
## +----------+------------+-------+-------+
## | | LAV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NNO | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | NOB | 32 | 0.2 |
## +----------+------------+-------+-------+
## | | NOR | 6754 | 42.0 |
## +----------+------------+-------+-------+
## | | PAN | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | PER | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | PHI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | POL | 14 | 0.1 |
## +----------+------------+-------+-------+
## | | POR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | RUS | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SCC | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | SCR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SIN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SME | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | SMI | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | SOM | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | SPA | 17 | 0.1 |
## +----------+------------+-------+-------+
## | | SWA | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SWE | 22 | 0.1 |
## +----------+------------+-------+-------+
## | | TAI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TAM | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | TGL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | THA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TIR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TUR | 9 | 0.1 |
## +----------+------------+-------+-------+
## | | URD | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | VIE | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 8970 | 55.8 |
## +----------+------------+-------+-------+
## | lnghomb | 000 | 6846 | 42.6 |
## +----------+------------+-------+-------+
## | | 888 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | 999 | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | ALB | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ARA | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | BOS | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | DAN | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | DUT | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | ENG | 75 | 0.5 |
## +----------+------------+-------+-------+
## | | FIN | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | FRE | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | GER | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | GRE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HIN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HUN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ICE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IRA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ITA | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | NNO | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NOR | 79 | 0.5 |
## +----------+------------+-------+-------+
## | | PAN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | PER | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | PHI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | POL | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | RUS | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | SGN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SME | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SMI | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SOM | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SPA | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | SWA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SWE | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | TAM | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | TGL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TUR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | UKR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | URD | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | VIE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 8970 | 55.8 |
## +----------+------------+-------+-------+
## | lnghom2 | 000 | 5522 | 34.4 |
## +----------+------------+-------+-------+
## | | 777 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | 888 | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | 999 | 2942 | 18.3 |
## +----------+------------+-------+-------+
## | | ALB | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ARA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BEM | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BEN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BOS | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | BUL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CHI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DAN | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | DUT | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | ENG | 187 | 1.2 |
## +----------+------------+-------+-------+
## | | FIL | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | FIN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | FRE | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | GER | 21 | 0.1 |
## +----------+------------+-------+-------+
## | | HEB | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | HIN | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | HRV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ICE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | ITA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | KUR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | LIT | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | LTZ | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NNO | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NOB | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | NOR | 169 | 1.1 |
## +----------+------------+-------+-------+
## | | PAN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | PER | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | PHI | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | POL | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | POR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | RUP | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | RUS | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | SMI | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | SOM | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SPA | 20 | 0.1 |
## +----------+------------+-------+-------+
## | | SRP | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | SWE | 12 | 0.1 |
## +----------+------------+-------+-------+
## | | TAM | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | TGL | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | THA | 5 | 0.0 |
## +----------+------------+-------+-------+
## | | TIG | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | TIR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TUR | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | URD | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | VIE | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 7095 | 44.2 |
## +----------+------------+-------+-------+
## | mbrncntb | 66 | 3995 | 24.9 |
## +----------+------------+-------+-------+
## | | 88 | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | 99 | 309 | 1.9 |
## +----------+------------+-------+-------+
## | | AU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BA | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | BR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | BY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CH | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CL | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | CN | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | CV | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | DE | 29 | 0.2 |
## +----------+------------+-------+-------+
## | | DK | 34 | 0.2 |
## +----------+------------+-------+-------+
## | | ER | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | ES | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | FI | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | FR | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | GB | 11 | 0.1 |
## +----------+------------+-------+-------+
## | | HU | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | IN | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IQ | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | IR | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | IS | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | IT | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | JP | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | LK | 4 | 0.0 |
## +----------+------------+-------+-------+
## | | LT | 10 | 0.1 |
## +----------+------------+-------+-------+
## | | MD | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | MY | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NL | 8 | 0.0 |
## +----------+------------+-------+-------+
## | | PH | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | PK | 9 | 0.1 |
## +----------+------------+-------+-------+
## | | PL | 48 | 0.3 |
## +----------+------------+-------+-------+
## | | RO | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | RU | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | SE | 49 | 0.3 |
## +----------+------------+-------+-------+
## | | SO | 7 | 0.0 |
## +----------+------------+-------+-------+
## | | TH | 2 | 0.0 |
## +----------+------------+-------+-------+
## | | TR | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | UA | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | US | 6 | 0.0 |
## +----------+------------+-------+-------+
## | | VN | 3 | 0.0 |
## +----------+------------+-------+-------+
## | | NA | 11460 | 71.3 |
## +----------+------------+-------+-------+
## | region | 99999 | 1 | 0.0 |
## +----------+------------+-------+-------+
## | | NO01 | 1716 | 10.7 |
## +----------+------------+-------+-------+
## | | NO02 | 651 | 4.1 |
## +----------+------------+-------+-------+
## | | NO03 | 1442 | 9.0 |
## +----------+------------+-------+-------+
## | | NO04 | 1115 | 6.9 |
## +----------+------------+-------+-------+
## | | NO05 | 1329 | 8.3 |
## +----------+------------+-------+-------+
## | | NO06 | 856 | 5.3 |
## +----------+------------+-------+-------+
## | | NO07 | 805 | 5.0 |
## +----------+------------+-------+-------+
## | | NO08 | 536 | 3.3 |
## +----------+------------+-------+-------+
## | | NO09 | 186 | 1.2 |
## +----------+------------+-------+-------+
## | | NO0A | 333 | 2.1 |
## +----------+------------+-------+-------+
## | | NA | 7095 | 44.2 |
## +----------+------------+-------+-------+
data_numeric <- Norway_data
summary_table <- datasummary_skim(data_numeric[c("ppltrst", "hinctnt", "hincfel")])
print(summary_table, title = "Table 1: Predictor Variables and Outcome variable")
##
## +---------+--------+--------------+------+------+-----+--------+------+--------------------------------------------------+
## | | Unique | Missing Pct. | Mean | SD | Min | Median | Max | Histogram |
## +=========+========+==============+======+======+=====+========+======+==================================================+
## | ppltrst | 14 | 0 | 6.8 | 3.4 | 0.0 | 7.0 | 99.0 |  |
## +---------+--------+--------------+------+------+-----+--------+------+--------------------------------------------------+
## | hinctnt | 15 | 65 | 10.9 | 14.0 | 1.0 | 9.0 | 88.0 |  |
## +---------+--------+--------------+------+------+-----+--------+------+--------------------------------------------------+
## | hincfel | 7 | 0 | 1.5 | 0.8 | 1.0 | 1.0 | 9.0 |  |
## +---------+--------+--------------+------+------+-----+--------+------+--------------------------------------------------+
library(infer)
library(ggplot2)
# Reminder the N0 is N0 <- "There is no significant relationship between income inequality and social trust in Norway."
# Choose a test statistic: mean difference between the two groups as our test statistic.
null_distribution <- replicate(1000, {
shuffled_data <- sample(c(Norway_data$ppltrst, Norway_data$hincfel))
mean_diff <- mean(shuffled_data[1:length(Norway_data$ppltrst)]) - mean(shuffled_data[(length(Norway_data$ppltrst) + 1):length(shuffled_data)])
return(mean_diff)
})
#Calculate and Visualize the null distribution
observed_statistic <- mean(Norway_data$ppltrst) - mean(Norway_data$hincfel)
# Histogram of the null distribution
null_plot <- ggplot() +
geom_histogram(aes(x = null_distribution), binwidth = 0.1, fill = "lightblue", color = "green") +
geom_vline(xintercept = observed_statistic, color = "black", linetype = "dashed") +
labs(title = "Distribution of Mean difference for H0") +
xlab("Mean") +
ylab("Frequency of Occurance")
print(null_plot)
# Regression Modelling
# Single Predictor Model
model_single <- lm(ppltrst ~ hincfel, data = Norway_data)
# Multivariate Regression Model
model_multivariate <- lm(ppltrst ~ hincfel + hinctnt, data = Norway_data)
# Interaction Model
model_interaction <- lm(ppltrst ~ hincfel * hinctnt, data = Norway_data)
# Regression Model Table: Single predictor model
single_model <- lm(ppltrst ~ hincfel, data = Norway_data)
single_table1 <- modelsummary(single_model, stars = TRUE)
print(single_table1)
##
## +-------------+------------+
## | | (1) |
## +=============+============+
## | (Intercept) | 7.068*** |
## +-------------+------------+
## | | (0.061) |
## +-------------+------------+
## | hincfel | -0.174*** |
## +-------------+------------+
## | | (0.036) |
## +-------------+------------+
## | Num.Obs. | 16065 |
## +-------------+------------+
## | R2 | 0.001 |
## +-------------+------------+
## | R2 Adj. | 0.001 |
## +-------------+------------+
## | AIC | 85276.0 |
## +-------------+------------+
## | BIC | 85299.1 |
## +-------------+------------+
## | Log.Lik. | -42635.012 |
## +-------------+------------+
## | RMSE | 3.44 |
## +=============+============+
## | + p < 0.1, * p < 0.05, |
## | ** p < 0.01, *** p < |
## | 0.001 |
## +=============+============+
# Regression Model Table: Multivariate Regression Model
multivariate_model <- lm(ppltrst ~ hincfel + hinctnt, data = Norway_data)
multivariate_table2 <- modelsummary(multivariate_model, stars = TRUE)
print(multivariate_table2)
##
## +-------------+------------+
## | | (1) |
## +=============+============+
## | (Intercept) | 6.867*** |
## +-------------+------------+
## | | (0.088) |
## +-------------+------------+
## | hincfel | -0.130** |
## +-------------+------------+
## | | (0.048) |
## +-------------+------------+
## | hinctnt | 0.005* |
## +-------------+------------+
## | | (0.003) |
## +-------------+------------+
## | Num.Obs. | 5546 |
## +-------------+------------+
## | R2 | 0.002 |
## +-------------+------------+
## | R2 Adj. | 0.002 |
## +-------------+------------+
## | AIC | 26474.7 |
## +-------------+------------+
## | BIC | 26501.2 |
## +-------------+------------+
## | Log.Lik. | -13233.344 |
## +-------------+------------+
## | RMSE | 2.63 |
## +=============+============+
## | + p < 0.1, * p < 0.05, |
## | ** p < 0.01, *** p < |
## | 0.001 |
## +=============+============+
#Regression Model Table: Interaction Model
interaction_model <- lm(ppltrst ~ hincfel * hinctnt, data = Norway_data)
interaction_table3 <- modelsummary(interaction_model, stars = TRUE)
print(interaction_table3)
##
## +-------------------+------------+
## | | (1) |
## +===================+============+
## | (Intercept) | 7.569*** |
## +-------------------+------------+
## | | (0.102) |
## +-------------------+------------+
## | hincfel | -0.534*** |
## +-------------------+------------+
## | | (0.056) |
## +-------------------+------------+
## | hinctnt | -0.041*** |
## +-------------------+------------+
## | | (0.004) |
## +-------------------+------------+
## | hincfel × hinctnt | 0.026*** |
## +-------------------+------------+
## | | (0.002) |
## +-------------------+------------+
## | Num.Obs. | 5546 |
## +-------------------+------------+
## | R2 | 0.032 |
## +-------------------+------------+
## | R2 Adj. | 0.031 |
## +-------------------+------------+
## | AIC | 26310.9 |
## +-------------------+------------+
## | BIC | 26344.0 |
## +-------------------+------------+
## | Log.Lik. | -13150.433 |
## +-------------------+------------+
## | RMSE | 2.59 |
## +===================+============+
## | + p < 0.1, * p < 0.05, ** p < |
## | 0.01, *** p < 0.001 |
## +===================+============+
# Another visual: Interaction Effect Visual
library(ggplot2)
interaction_plot <- ggplot(Norway_data, aes(x = hincfel, y = ppltrst, color = hinctnt)) +
geom_point() +
stat_smooth(method = "lm", se = FALSE, aes(group = hinctnt)) +
labs(x = "Living Comfortably", y = "Social Trust", color = "Mean Income") +
theme_minimal()
print(interaction_plot)
## `geom_smooth()` using formula = 'y ~ x'