For at finde et eksempel på den marginale forskel i sandsynligheden for tilbagefald til kriminalitet for en person, der får samfundstjeneste sammenlignet med en person, der får fængselsstraf, har vi taget udgangspunkt i en person, der har gennemsnitsværdier på metriske og dikotome variable og er i referencekategorien for kategoriske variable med flere end to værdier. Gennemsnittene er fundet via forskellige tabeller i rapporten. Vi har hevet koefficienter ud fra den første model i rapportens bilagsmateriale.
Nedenfor har vi kopieret R-koden ind. Koden starter her:
##Aflæs gennemsnitsværdier fra tabeller i rapporten
##Kategoriske variable med mere end to kategorier sættes til referencekategorien
##Kategoriske variable med to kategorier sættes til deres gennemsnit
alder <- 30
##Alder: Sat til 2006 (midt i observationsperioden) - 1976 (gennemsnitligt fødeår). Ens for grupperne
tidkrim <- (1021*11 + 589*16.1)/(1021+589)
##Tidligere kriminelle forhold vægtet efter gruppernes størrelse
domlængde <- (1021*11.9 + 589*22.9)/(1021+589) ##Domslængde
alkm <- (1021*7 + 589*10)/(100*(1021+589)) ##Alkoholmisbrug
narko <- 0.08 ##Stofmisbrug (ens i grupperne)
mænd <- (1021*87 + 589*95)/(100*(1021+589)) ##Mand el kvinde
besk <- 1 ##I beskæftigelse
grund <- 1 ##Uddannelse = Grundskole
enligU <-1 ##Enlig uden børn
gift <- 1 ##Gift
kbh <- 1 ##Afdeling for vurdering = København
emi <- (1021*13 + 589*17)/(100*(1021+589)) ##Emigreret efter 15. år
## Saml kontrolvariable og tilføj 1 for konstantled
vars <- c(alder,tidkrim,domlængde,alkm,narko,mænd,besk,grund,enligU,gift,kbh,emi,1)
##Hiv koefficienter ind fra tabel 5.1.
coefs <- c(0.722,0.932,1.014,1.002,1.263,1.086,
1.878,1,1,1,1,1,0.928,1.757)
##Tag logaritme til koefficienterne for at kunne beregne som lineær funktion
lncoefs <- log(coefs)
##Find forudsagt værdi for samfundstjeneste
logitSamf <- lncoefs %*% c(1,vars)
##Find forudsagt sandsynlighed
pSamf <- 1/(1+exp(-logitSamf))
##Gen tag for fængselsgruppen
logitKrim <- lncoefs %*% c(0,vars)
pKrim <- 1/(1+exp(-logitKrim))
#Find marginal sandsynlighed
pSamf - pKrim ## i pct.point
(pSamf - pKrim)/pKrim
