% time1470 and data1470 available for use
% for Amazon best-seller's rank for #51-ranked book on 2/22/12
% 
% time1469 and data1469 available for use
% for Amazon best-seller's rank for #41-ranked book on 2/22/12
% 
% time1468 and data1468 available for use
% for Amazon best-seller's rank for #31-ranked book on 2/22/12
% 
% time1467 and data1467 available for use
% for Amazon best-seller's rank for #21-ranked book on 2/22/12
% 
% time1466 and data1466 available for use
% for Amazon best-seller's rank for #11-ranked book on 2/22/12
% 
% Amazon best-seller's rank for #1-ranked book on 2/22/12
%   time vector is: time1465
%   data vector is: data1465

for k = 65:70
  n = ['14',num2str(k)];
  eval(['idx = data',n,'==0;']);
  eval(['data',n,'(idx) = [];']);
  eval(['time',n,'(idx) = [];']);
end
% data1465(data1465==0) = NaN;
% data1466(data1466==0) = NaN;
% data1467(data1467==0) = NaN;
% data1468(data1468==0) = NaN;
% data1469(data1469==0) = NaN;
% data1470(data1470==0) = NaN;

subplot(2,1,1)
plot(time1465,data1465,time1466,data1466,time1467,data1467,time1468,data1468,time1469,data1469,time1470,data1470);
datetick
xlabel('Date'), ylabel('Rank')
set(gca,'YDir','reverse')
subplot(2,1,2)
plot(time1465,data1465,time1466,data1466,time1467,data1467,time1468,data1468,time1469,data1469,time1470,data1470);
datetick
ymx = median([max(data1465),max(data1466),max(data1467),max(data1468),max(data1469),max(data1470)]);
ylim([0,100*ceil(ymx/100)])
xlabel('Date'), ylabel('Rank')
set(gca,'YDir','reverse')