clrsmp
Clearance-rate decomposition and Simpson’s paradox
clrsmp(clra, reca, clrb, recb)crimlab.clrsmp(clra, reca, clrb, recb)crimlab clrsmp --clra VALUE --reca VALUE --clrb VALUE --recb VALUESplits the change in an aggregate clearance rate between two periods,
areas or services (A and B) into a rate effect, from each offence type
clearing more or less, and a mix effect, from the workload moving between
types that clear at different rates (Kitagawa’s decomposition). The two
always add up to the change. When every type clears better in B and the
aggregate still falls, the fall is all mix: Simpson’s paradox, and the
reason an aggregate clearance rate cannot be read as performance without
the split by type. See proof("clrsmp").
Parameters
Section titled “Parameters”clraOffences cleared in A, by type: non-negative numbers, as a numeric vector, a 1-d table, a list of single numbers or a one-row data frame.
recaOffences recorded in A, by type: positive numbers, in the forms of
clra; named inputs are matched toclraby name, unnamed ones by position (only names protect against a mismatch).clrbOffences cleared in B, by type: non-negative numbers, in the forms of
clra, matched by name.recbOffences recorded in B, by type: positive numbers, in the forms of
clra, matched by name.
Returns
Section titled “Returns”A list. parts, a data frame with one row per type, in the order
of clra: part,
ratea, rateb (clearance rates), sharea, shareb (shares of the
recorded workload), rateeff and mixeff. total, a named numeric
vector: ratea, rateb (aggregate rates), change, rateeff,
mixeff (with change = rateeff + mixeff), up and down (types
whose rate rose or fell) and simpson: 2 when every type moved one way
and the aggregate the other, 1 when the aggregate change and the rate
effect differ in sign, 0 otherwise.
Each part’s mix term is measured from the mean of the two aggregate rates, so a part’s term is negative when work moved toward a part that clears below average; the centre changes how the mix effect is shared among the parts, never its total.
References
Section titled “References”Kitagawa, E. M. (1955). Components of a difference between two rates. Journal of the American Statistical Association, 50(272), 1168-1194. doi:10.1080/01621459.1955.10501299
Examples
Section titled “Examples”# both types clear better in B, yet the aggregate fallsclrsmp(clra = c(minor = 90, serious = 20), reca = c(100, 50), clrb = c(19, 60), recb = c(20, 140))Generated from the R help page: clrsmp in the R reference.