% import_starych_danych;
rows = owidcoviddata.iso_code == 'POL';
Polska = owidcoviddata(rows,:);

%
% Teraz skasujemy niepotrzebne zmienne rows i owidcoviddata

clear rows owidcoviddata
figure
plot(Polska.date, Polska.new_cases);

%
% Wykres jest „poszarpany”, więc skorzystamy z wartości wygładzonej

figure
plot(Polska.date, Polska.new_cases_smoothed);

%
% Wybierzemy teraz dane od  1 września do końca roku 2020

r2020 = Polska.date >= datetime(2020, 9, 1) & Polska.date <= datetime(2020,12,31);

% Teraz to samo dla roku 2021 i kolejnych lat

r2021 = Polska.date >= datetime(2021, 9, 1) & Polska.date <= datetime(2021,12,31);
r2022 = Polska.date >= datetime(2022, 9, 1) & Polska.date <= datetime(2022,12,31);
r2023 = Polska.date >= datetime(2023, 9, 1) & Polska.date <= datetime(2023,12,31);
r2024 = Polska.date >= datetime(2024, 9, 1) & Polska.date <= datetime(2024,12,31);

%
% Chciałbym porównać uśrednione wykresy z pięciu lat. Wybierzemy dane za 2020

figure
p2020 = Polska(r2020,:);
p2021 = Polska(r2021,:);
p2022 = Polska(r2022,:);
p2023 = Polska(r2023,:);
p2024 = Polska(r2024,:);
plot(p2020.date,p2020.new_cases_smoothed)
title('2020')
figure
plot(p2021.date,p2021.new_cases_smoothed)
title('2021')
figure
plot(p2022.date,p2022.new_cases_smoothed)
title('2022')
figure
plot(p2023.date,p2023.new_cases_smoothed)
title('2023')
figure
plot(p2024.date,p2024.new_cases_smoothed)
title('2024')

%
% za rok 2024 nie ma już żadnnych danych w tym okressie
% To co chciałbym osiągnąć to oba zestawy danych na jednym wykresie (te same jednostki na osiach). Pierwszy problem jest taki, że osie czasowe danych są różne…
% Kolejnym problemem jest to, że zestawy danych są różnej długości. Najprościej skrócić zestaw dłuższy

%a4=size(p2024)
a3=size(p2023);
a2=size(p2022);
a1=size(p2021);
a0=size(p2020);
N = min([a3(1,1), a2(1,1), a1(1,1), a0(1,1)]);

%
% Teraz możemy rysować wykres danych. Jako współrzędne czasowe możemy przyjąć daty z 2020 lub 2021 albo 2022 — zawsze będzie drobne zafałszowanie…

figure
plot(p2020.date(1:N), p2020.new_cases_smoothed(1:N), p2020.date(1:N), p2021.new_cases_smoothed(1:N), p2020.date(1:N), p2022.new_cases_smoothed(1:N), p2020.date(1:N), p2023.new_cases_smoothed(1:N))
legend('2020','2021','2022','2023','Location',"best")

%
% A teraz średnia z M dni

M = 30;
plot(p2020.date,movmean(p2020.new_cases,M))
title('2020')
plot(p2021.date,movmean(p2021.new_cases,M))
title('2021')
plot(p2022.date,movmean(p2022.new_cases,M))
title('2022')
plot(p2023.date,movmean(p2023.new_cases,M))
title('2023')
figure
plot(p2020.date(1:N), movmean(p2020.new_cases(1:N),M), p2020.date(1:N), movmean(p2021.new_cases(1:N),M),p2020.date(1:N), movmean(p2022.new_cases(1:N),M), p2020.date(1:N), movmean(p2023.new_cases(1:N),M))
title("Średnia liczona inaczej")
legend('2020','2021','2022','2023','Location',"best")

% Na koniec zapiszemy tablicę p2021 w pliku do innych analiz i prezentacji

save('p2021.mat',"p2021")

%
% A teraz liczba zgonów

figure
plot(p2020.date(1:N), p2020.new_deaths_smoothed(1:N), p2020.date(1:N), p2021.new_deaths_smoothed(1:N), p2020.date(1:N), p2022.new_deaths_smoothed(1:N),p2020.date(1:N), p2023.new_deaths_smoothed(1:N))
legend('2020', '2021','2022','2023','Location',"best")
title("Liczba zgonów")

%
%
%
%
%
%