calc_stats_2x3 function

Calculate p-values for 2x3-way analysis