My dataset looks like this
df = data.frame(site=c(rep(\'A\',95),rep(\'B\',110),rep(\'C\',250)), nps_score=c(floor(runif(455, min=0, max=10))),