Consider the two following data frames:
df <- data.frame(REGION = c("REG01","REG02","REG03","REGSUM"),