reporte/plots/fig-gen.gp

# 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"
f(x) = a*exp(b*x)
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
# 6 Modelo lineal Michoacan
clear
file2="../datos/datos-covid-michoacan.csv"
set output "1-6.png"
f(x) = a*x+b
fit f(x) file2  using 2:3 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("f(x)=%+-f*x%+-f", a, b) at graph 0.01,0.85 font ',12'
set title 'Figura 6: Modelo lineal para los casos en Michoacán' font ',14'
show title
set key left top font ',10'  box 3
plot file2 using 2:3 title 'Casos' w p ls 1, \
     f(x) with lines ls 2 title 'Modelo'
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