R code to calculate dual-grouped median and arrange by descending medians

I have a data frame consisting of 3 columns: Site, Program, Result. Here a minimal repro dataset:

R code to calculate grouped median, arrange by descending medians, then re-arrange

I have a data frame consisting of 3 columns: Site, Program, Result. Here a minimal repro dataset: