Udregninger til kritik af rapport fra Justitsministeriets forskningskontor (http://politiskkompas.blogspot.com/2014/11/udregninger-til-kritik-af-rapport-fra.html)

udregninger.r

##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