mean
Deciles
intersect
FOO <- function(x){
SAP <- x$mean[x$Population == "SAP"]
names(SAP) <- x$Deciles[x$Population == "SAP"]
SH <- x$mean[x$Population == "SH"]
names(SH) <- x$Deciles[x$Population == "SH"]
TD6 <- x$mean[x$Population == "TD6"]
names(TD6) <- x$Deciles[x$Population == "TD6"]
SAPvsSH <- sum(SH)*100 / sum(SAP[intersect(names(SAP), names(SH))]) - 100
SAPvsTD6 <- sum(TD6)*100 / sum(SAP[intersect(names(SAP), names(TD6))]) - 100
SHvsTD6 <- sum(TD6[intersect(names(SH), names(TD6))])*100 / sum(SH[intersect(names(SH), names(TD6))]) - 100
}