I created the following boxlopt:
boxplot(LCZ_fm[,'SVF_mean'] ~ LCZ_fm[,'LCZ'], main='Sky View Factor pro LCZ-Klasse', ylab='SVF', xlab='LCZ', na.action = NULL )
Now I would like to color specific ranges for each boxplot (like I drafted in green), showing the 'optimal range' of values for each group ('LCZ', on the x axis) according to literature, to illustrate how my data deviates from these theoretical values.
Any suggestions? Could I maybe use some sort of faceting technique from ggplot2?
Bonus question: My goal is to regroup my data, so that the values fit better into the theoretical ranges. Meaning, some data points are going to be assigned to a new group (LCZ). Does someone have an idea how I could compare the before and after data using statistics to quantify how much more accurate the grouping is, i.e. if the deviation from theoretical value ranges has decreased overall?