s1Logical = month(giulia_TT.Time) == 12 | month(giulia_TT.Time) == 1;
s2Logical = month(giulia_TT.Time) >= 3 & month(giulia_TT.Time) <= 10;
s3Logical = month(giulia_TT.Time) == 11;
s4Logical = month(giulia_TT.Time) == 2;
sAll = s1Logical*1 + s2Logical*2 + s3Logical*3 + s4Logical*4;
T = timetable2table(giulia_TT);
data = T(:,{'Var5','Season'});
statarray = grpstats(data, 'Season', {'mean' 'std'})
yMean = statarray.mean_Var5;
ySTD = statarray.std_Var5;
b = bar(x, yMean, 'facecolor', [.8 .8 .9]);
set(gca, 'YLim', [0 120]);
eb = errorbar(x, yMean, ySTD, ySTD);
s = sprintf('%.2f,', statarray.mean_Var5);
text(b.XEndPoints, b.YData * 0.4, s, 'b', 'fontsize', 12, 'horizontalalignment', 'center');