Daily_reported_cases.jl

############################# Calculate number of daily reported cases ####################################


y_tmp = deepcopy(y);               # creates an independent copy, changes in y_tmp won't affect y
rows,cols = size(y_tmp)
dfrows = nrow(y_tmp);
name = names(y_tmp);

y_daily = similar(y_tmp,dfrows-1); # copy the structure to an empty dataframe with dfrows-1 rows

for j = 1:length(name)
    for i = 1:dfrows-1
        y_daily[!,name[j]][i] = y_tmp[!,name[j]][i+1] - y_tmp[!,name[j]][i] # calculate number of daily increase in reported cases by subtracting numbers of previous day
    end
end

x_daily = deepcopy(x);
popfirst!(x_daily);                # remove first entry, daily increase can only be reported from the next date

gr(size=(900,600))
display(@df y_daily bar(x_daily, cols(1:cols), 
    label = reshape(names(y),(1,length(names(y)))),
    xlabel = "Time",
    ylabel = "Daily number of reported cases",
    xticks = x[1:7:end],
    xrotation = 45,
    legend = :topleft,
    grid = false,
    framestyle = :semi,
    legendfontsize = 9,
    formatter = :plain))

savefig("Daily_cases.png")