Calculating Median and Mode of numerical variables for different subgroups in R
assuming your data is stored in an object named df you can do:

tapply(df$S_Calls, df$Emp_Stat, median)


As for the mode, oddly enough R does not have a built in function for that. You could define one yourself using: 

mode_stat &lt;- function(x) {
  ux &lt;- na.omit(unique(x))
  ux[which.max(tabulate(match(x, ux)))]
}


and then do in a similar fashion:

tapply(df$S_Calls, df$Emp_Stat, mode_stat)