Udregninger til kritik af rapport fra Justitsministeriets forskningskontor

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