% Unanswered
t1 = time1527(~isnan(data1527));
d1 = data1527(~isnan(data1527));

% Answered/Not Accepted
t2 = time1528(~isnan(data1528));
d2 = data1528(~isnan(data1528));

% Answered/Accepted
t3 = time1529(~isnan(data1529));
d3 = data1529(~isnan(data1529));

tn = unique([t1; t2; t3]);
tn(tn < max([t1(1) t2(1) t3(1)])) = [];
tn(tn > min([t1(end) t2(end) t3(end)])) = [];
d1n = interp1(t1,d1,tn);
d2n = interp1(t2,d2,tn);
d3n = interp1(t3,d3,tn);

area(tn,[d3n d2n d1n],'EdgeColor',[215 215 215]/255);
datetick
legend({'Answered/Accepted','Answers/Not Accepted','Unanswered'},'Location','NorthOutside')
colormap([178 171 115; 235 200 121; 201 118 92]/255)
box off
% xlim([now-90 now])
dall = d1(end)+d2(end)+d3(end);
text(tn(end),d3(end), ...
  sprintf(' %.1f%%\n answered &\n accepted',100*d3(end)/dall))
text(tn(end),d2(end)+d3(end), ...
  sprintf(' %.1f%%\n answered',100*(d2(end)+d3(end))/dall))