# Preparing data
clear
reset
set datafile separator ","
set datafile missing "NaN"
set bmargin at screen 0.15
set lmargin at screen 0.15
file="../datos/datos-covid-mexico.csv"
set linestyle 1 lc rgb 'black' pt 6 ps 2  # circle
set linestyle 2 lt 0 lw 3.3
FIT_LIMIT = 1e-6
# 1 Expo
set terminal png
# set output "1-1.png"
# f(x) = a*exp(b*x)
# fit f(x) file using 2:3 via a,b
# set xlabel 'Días de Contagio' font ',12'
# set xtics font ", 11"
# show xlabel
# set ytics font ", 11"
# set ylabel 'Casos Confirmados' font ',14'
# set label sprintf("f(x)=%f*e(%f*x)", a, b) at graph 0.01,0.85 font ',12'
# set title 'Figura 1: Modelo Exponencial.' font ',16'
# show title
# set key left top font ',10'  box 3
# plot [1:][]file using 2:3 title 'Casos' w p ls 1, \
#      f(x) with lines ls 2 title 'Exponencial'
# set output 
# unset output
# unset label
# # 2 Polinomico
# clear
# set output "1-2.png"
# f(x) = a*x*x+b*x+c
# fit f(x) file using 2:3 via a,b,c
# set xtics font ", 13"
# show xlabel
# set ytics font ", 11"
# set ylabel 'Casos confirmados' font ',14'
# set label sprintf("f(x)=%f*x^2-%f*x+%f", a, b, c) at graph 0.01,0.85 font ',12'
# set title 'Figura 2: Modelo polinomial' font ',16'
# show title
# set key left top font ',10'  box 3
# plot file using 2:3 title 'Casos' w p ls 1, \
#      f(x) with lines ls 2 title 'Polinomial'
# set output 
# unset output
# unset label
# 3 Logaritmico 
clear
set output "1-3.png"
f(x) = a*x+b
fit f(x) file using 2:9 via a, b
set xlabel 'Días de contagio' font ',14'
set xtics font ", 13"
show xlabel
set ytics font ", 11"
set ylabel 'Casos confirmados' font ',14'
set label sprintf("log(f(x))=%f*x+%f", a, b) at graph 0.01,0.85 font ',12'
set title 'Figura 3: Modelo logarítmico' font ',16'
show title
set key left top font ',10'  box 3
plot file using 2:(log($3)) title 'Casos' w p ls 1, \
     f(x) with lines ls 2 title 'Log'
set output 
unset output
unset label
# 4 Factor de crecimiento mexico
clear
set output "1-4.png"
set linestyle 1 lc rgb 'black' pt 6 ps 2  # circle
set linestyle 2 lt 0 lw 3.3
set xlabel 'Días de contagio' font ',14'
set xtics font ", 13"
show xlabel
set ytics font ", 11"
set ylabel 'Factor de crecimiento' font ',14'
set title 'Figura 4: Factor de crecimeinto díario de COVID19 en México' font ',14'
show title
set key off
plot [1:][] file using 2:4 w lp ls 1
set output
unset output
unset label
# 5 Incremento diario
clear
set output "1-5.png"
g(x) = c*exp(d*x)
fit g(x) file using 2:6 via c, d
set xlabel 'Días de contagio' font ',14'
set xtics font ", 13"
show xlabel
set ytics font ", 11"
set ylabel 'Incremento de casos diario' font ',14'
set label sprintf("f(x)=%+-f*exp(%+-f*x)", c, d) at graph 0.01,0.85 font ',12'
set title 'Figura 5: Modelo lineal para incrmento de los casos diarios en México' font ',14'
show title
set key left top font ',10'  box 3
plot file using 2:6 title 'Incremento' w lp ls 1, \
     g(x) with lines ls 2 title 'Regresor'
set output 
unset output
unset label

# Factor de crecimiento Michoacan
clear
set output "1-7.png"
set xlabel 'Días de contagio' font ',14'
set xtics font ", 13"
show xlabel
set ytics font ", 11"
set ylabel 'Factor de crecimiento' font ',14'
set title 'Figura 7: Factor de crecimeinto díario de COVID19 en Michoacán' font ',14'
show title
set key off
plot [1:][] file2 using 2:6 w lp ls 1
set output 
unset output
unset label