New Cases in Italy

p = plot(new_national$nuovi_positivi ~ new_national$data, type="l", lwd=6, pch=16, cex=1.6, col=c("#3B3176"))
#text(new_national$nuovi_positivi ~ new_national$data, labels=new_national$nuovi_positivi, pos=2, cex=1.3)
abline(v = as.Date("2020-10-25"), col="#FF0000", lwd=3)

grid(p, col = "black", lty="dotted")
rect(as.Date("2020-03-09"), -200, as.Date("2020-05-19"), 7000, border = "#AA0055", col = "#AA0055", density=2)
abline(h = new_national[27,]$nuovi_positivi, col="#FF0000", lwd=3)

<<EOS
<div id="heatmap-total"></div>
<script>
calendar_heatmap.create({
    data: [
    #{data.map { |x| "{ day: \"%s\", count: %d }" % [x[0], x[1].to_i] }.join(",") }
    ],
    date_var: "day",
    fill_var: "count",
    target: "#heatmap-total"
});
</script>
EOS

require 'date'

def find_new_cases_for date, table
  table.each do |row|
    return row[1] if date == Date.parse(row[0])
  end
  return ""
end

year_grid = Array.new
legenda = (0..31).map { |x| x }
legenda[0] = ""
year_grid << legenda

(1..12).each do |month|
  month_row = Array.new(32)
  month_row[0] = Date.new(2020, month, 1).strftime("%b")
  (1..31).each do |day|
    if Date.valid_date?(2020, month, day) 
      date = Date.new(2020, month, day)
      cases = find_new_cases_for(date, table)
      month_row[day] = cases
    else
      month_row[day] = ""
    end
  end
  year_grid << month_row
end
year_grid

current_month=$(date +"%b")
current_month_no=$(date +"%m")

egrep 2020-${current_month_no} ./data/dpc-covid19-ita-regioni.csv | sed -e "s/T17:00:00//g" | sed -e "s/2020-${current_month_no}-/${current_month} /g" | datamash -t, crosstab 4,1 sum 13

# how many rows and columns?
par(mfrow=c(11, 2))
regions <- c("Valle d'Aosta", "Piemonte", "Liguria", "Lombardia", "Veneto",
             "P.A. Trento", "P.A. Bolzano", "Friuli Venezia Giulia",
             "Emilia-Romagna", "Toscana", "Marche", "Umbria",
             "Lazio", "Abruzzo", "Molise", "Campania",
             "Puglia", "Basilicata", "Calabria", "Sicilia",
             "Sardegna")
for (region in regions) {
  my_plot(region,
          data, 
          max=max(data$nuovi_positivi),
          filter=TRUE, 
          textlabels=FALSE,
          variables=c("nuovi_positivi"),
          graphtypes=c("l"),
          colors=c("black"))
}