x = [120, 110.5, 123, 140, 142, 151, 121, 119, 110, 115, 130]';
y = [118, 110, 120, 135, 134, 140, 120, 116, 110, 113 ,127]';
g = [0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 1]';
m = [x y g];
t = array2table(m,'VariableNames',{'x','y','gender'});
t.gender = categorical(t.gender);
timeDesign = table([1 2]','VariableNames',{'time'});
timeDesign.time = categorical(timeDesign.time);
rm = fitrm(t,'x-y~gender','WithinDesign',timeDesign);
ranova(rm,'WithinModel','time')