% uiopen('owid-covid-data.csv',1)
countries(1)=owidcoviddata.location(1);
c=1;
% d(c)=1;
for i=2:size(owidcoviddata,1)
    if owidcoviddata.location(i) ~= countries(c)
        c=c+1;
        countries(c)=owidcoviddata.location(i);
    %     d(c)=1;
    % else 
    %     d(c)=d(c)+1;
    end
    % deaths(c,d(c))=owidcoviddata.new_deaths_smoothed_per_million(i);
    % dates(c,d(c))=owidcoviddata.date(d(c));
end


for i=1:size(countries,2)
    kk=owidcoviddata(find(owidcoviddata.location==countries(i)),:);
    d(i)=size(kk,1);
    deaths(i,1:size(kk,1))=kk.new_deaths_smoothed;
    dates(i,1:size(kk,1))=kk.date;
end


deaths(isnan(deaths))=0;

[m,c]=max(d);
dates_voc=dates(c,:);

vocabulary=string(dates_voc);


figure
plot(dates_voc,deaths')

mdl=fitlda(floor(deaths),5);

plot(mdl.TopicWordProbabilities)