Avansert logistisk regresjonsanalyse
Eksempelet illustrerer hvordan man kan gå frem for å analysere sannsynligheten for å komme i jobb og tjene over 500000kr ett år etter at en befinner seg i en tilstand uten jobb. Aldersgruppen vi ser på er 16-60.
I analysen blir det kontrollert for diverse demografiske kjennetegn samt status på arbeidsmarkedet (arbeidsledig, ordinære arbeidsmarkedstiltak, yrkeshemmet, andre arbeidssøkertilstander, samt arbeidsuførhet).
Det lages først en del deskriptiv statistikk, og tilslutt kjøres en logit-analyse inkludert marginaleffekter (opsjonen mfx(dydx)
brukes til dette).
//Kobler til databank
require no.ssb.fdb:23 as db
//Oppretter populasjonen personer 16-60 år som er uten jobb i november 2018, og bosatt i Norge per 1. januar 2019
create-dataset demografidata
import db/BEFOLKNING_FOEDSELS_AAR_MND as faarmnd
import db/BEFOLKNING_STATUSKODE 2019-01-01 as regstat
import db/REGSYS_ARB_ARBMARK_STATUS 2018-11-16 as arbmstatus
generate alder = 2018 - int(faarmnd / 100)
generate jobb = 0
replace jobb = 1 if arbmstatus == '1' | arbmstatus == '2'
keep if alder >= 16 & alder <= 60 & regstat == '1' & jobb == 0
histogram alder, discrete
//Henter øvrige variabler (demografidata finnes det bare data på per 1/1 hvert år)
import db/BEFOLKNING_KJOENN as kjønn
import db/BEFOLKNING_INVKAT as innvkat
import db/SIVSTANDFDT_SIVSTAND 2018-11-16 as sivstand
import db/BEFOLKNING_BARN_I_REGSTAT_FAMNR 2019-01-01 as antbarn
import db/NUDB_BU 2018-11-16 as utd
import db/NUDB_SOSBAK as sosbakgrunn
import db/BEFOLKNING_KOMMNR_FAKTISK 2019-01-01 as bosted
import db/ARBSOEK2001FDT_HOVED 2018-11-16 as arbsøkerstatus
import db/UFOERP2011FDT_GRAD 2018-11-16 as uføregrad
import db/INNTEKT_BRUTTOFORM 2018-12-31 as formue
import db/INNTEKT_WYRKINNT 2019-12-31 as inntekt19
//Lager en avhengig variabel med to utfall (dummyvariabel): Høy inntekt vs. lav inntekt
histogram inntekt19, width(100000) freq
summarize inntekt19
generate høyinnt = 0
replace høyinnt = 1 if inntekt19 > 500000
piechart høyinnt
//Tilrettelegger de uavhengige variablene slik at de passer med den statistiske modellen (innebærer at de fleste variabler gjøres som til dummy’er)
generate mann = 0
replace mann = 1 if kjønn == '1'
piechart mann
destring sivstand
generate gift = 0
replace gift = 1 if sivstand == 2
replace gift = sivstand if sysmiss(sivstand)
piechart gift
generate innvandrer = 0
replace innvandrer = 1 if innvkat == 'B'
piechart innvandrer
tabulate antbarn, missing
generate ettbarn = 0
replace ettbarn = 1 if antbarn == 1
generate flerebarn = 0
replace flerebarn = 1 if antbarn > 1
destring utd
generate høyutd = 0
replace høyutd = 1 if utd >= 700000 & utd < 900000
replace høyutd = utd if sysmiss(utd)
piechart høyutd
generate høyutd_foreldre = 0
replace høyutd_foreldre = 1 if sosbakgrunn == '1'
piechart høyutd_foreldre
generate oslo = 0
replace oslo = 1 if bosted == '0301'
generate bergen = 0
replace bergen = 1 if bosted == '1201'
generate stavanger = 0
replace stavanger = 1 if bosted == '1103'
generate trondheim = 0
replace trondheim = 1 if bosted == '5001'
barchart(sum) oslo bergen stavanger trondheim
destring arbsøkerstatus
tabulate arbsøkerstatus, missing
generate ledig = 0
replace ledig = 1 if arbsøkerstatus == 1
generate tiltak = 0
replace tiltak = 1 if arbsøkerstatus == 3
generate yrkeshemmet = 0
replace yrkeshemmet = 1 if arbsøkerstatus == 5 | arbsøkerstatus >= 10
generate andrearbsøk = 0
replace andrearbsøk = 1 if arbsøkerstatus == 2 | arbsøkerstatus == 4 | arbsøkerstatus == 7
generate ufør = 1
replace ufør = 0 if sysmiss(uføregrad)
barchart(sum) ledig tiltak yrkeshemmet andrearbsøk ufør
histogram formue, width(100000) freq
summarize formue
generate formuehøy = 0
replace formuehøy = 1 if formue > 1000000
replace formuehøy = formue if sysmiss(formue)
piechart formuehøy
//Bruk sankey for å vise overganger
sankey arbsøkerstatus høyinnt
sankey høyutd høyinnt
//Kjører selve logit-analysen der den avhengige variabelen listes først (må være dummy)
logit høyinnt mann gift alder innvandrer ettbarn flerebarn høyutd høyutd_foreldre oslo bergen stavanger trondheim ledig tiltak yrkeshemmet andrearbsøk ufør formuehøy, mfx(dydx)