Presuming you have the Statistics Toolbox which defines the two statistics of interest m-files--
G = groupsummary(T,groupvars,method,datavars)
where
method={'skewness','kurtosis'};
You'll define the other inputs to match your wishes on grouping and data variables.
If you don't have the requisite TB, then write your own functions and refer to them in their place.
