Modelado prediccion enfermedades cultivos - edUTecNe ...

Modelado para la predicción de enfermedades en cultivos de alto valor comercial
data Conc1;
SET c1;by anio dj;
tmed=(tmax+tmin)/2;
pr1=lag1(prec);pr2=lag2(prec);
pr3=lag3(prec);pr4=lag4(prec);pr5=lag5(prec);pr6=lag6(prec);pr7=lag7(p
rec);pr8=lag8(prec);pr9=lag9(prec);pr10=lag10(prec);pr11=lag11(prec);p
r12=lag12(prec);pr13=lag13(prec);pr14=lag14(prec);
hrd1=lag1(hrd);hrd2=lag2(hrd);hrd3=lag3(hrd);hrd4=lag4(hrd);hrd5=lag5(
hrd);hrd6=lag6(hrd);hrd7=lag7(hrd);hrd8=lag8(hrd);hrd9=lag9(hrd);hrd10
=lag10(hrd);hrd11=lag11(hrd);hrd12=lag12(hrd);hrd13=lag13(hrd);hrd14=l
ag14(hrd);
tmin1=lag1(tmin);tmin2=lag2(tmin);tmin3=lag3(tmin);tmin4=lag4(tmin);tm
in5=lag5(tmin);tmin6=lag6(tmin);tmin7=lag7(tmin);tmin8=lag8(tmin);tmin
9=lag9(tmin);tmin10=lag10(tmin);tmin11=lag11(tmin);tmin12=lag12(tmin);
tmin13=lag13(tmin);tmin14=lag14(tmin);
tmax1=lag1(tmax);tmax2=lag2(tmax);tmax3=lag3(tmax);tmax4=lag4(tmax);tm
ax5=lag5(tmax);tmax6=lag6(tmax);tmax7=lag7(tmax);tmax8=lag8(tmax);tmax
9=lag9(tmax);tmax10=lag10(tmax);tmax11=lag11(tmax);tmax12=lag12(tmax);
tmax13=lag13(tmax);tmax14=lag14(tmax);
run;
data Fc2;set conc1;
DPrHRT =0;DTxnN=0;DTnN=0;DPrHR=0;FPr=0; DTxnP=0;
DHR=0; retain DPrHRT DTxnP DTxnN DTnN DPrHR FPr DHR;
TnM=(tmin1+tmin2+tmin3+tmin4+tmin5+tmin6+tmin7+tmin8+tmin9+tmin10+tmin
11+tmin12+tmin13+tmin14)/14;
TxM=(tmax1+tmax2+tmax3+tmax4+tmax5+tmax6+tmax7+tmax8+tmax9+tmax10+tmax
11+tmax12+tmax13+tmin14)/14;
if hrd1>65 and pr1>0.2 then do;if tmax116 then
DPrHRT+1;end;if hrd2>65 and pr2>0.2 then do;if tmax216
then DPrHRT +1;end;if hrd3>65 and pr3>0.2 then do;if tmax316 then DPrHRT +1;end;if hrd4>65 and pr4>0.2 then do;if tmax416 then DPrHRT +1;end;if hrd5>65 and pr5>0.2 then do;if
tmax516 then DPrHRT +1;end;if hrd6>65 and pr6>0.2 then
do;if tmax616 then DPrHRT +1;end;if hrd7>65 and pr7>0.2
then do;if tmax716 then DPrHRT +1;end;if hrd8>65 and
pr8>0.2 then do;if tmax816 then DPrHRT +1;end;
if hrd9>65 and pr9>0.2 then do;if tmax916 then DPrHRT
+1;end;if hrd10>65 and pr10>0.2 then do;if tmax1016
then DPrHRT +1;end;if hrd11>65 and pr11>0.2 then do;if tmax1116 then DPrHRT +1;end;if hrd12>65 and pr12>0.2 then do;if
tmax1216 then DPrHRT +1;end;if hrd13>65 and pr13>0.2
then do;if tmax1316 then DPrHRT +1;end;if hrd14>65 and
pr14>0.2 then do;if tmax1416 then DPrHRT +1;end;
if tmax116 then DTxnP+1;if tmax216 then
DTxnP+1;if tmax316 then DTxnP+1;if tmax416
then DTxnP+1;if tmax516 then DTxnP+1;if tmax616 then DTxnP+1;if tmax716 then DTxnP+1;if
tmax816 then DTxnP+1;if tmax916 then
DTxnP+1;if tmax1016 then DTxnP+1;if tmax1116 then DTxnP+1;if tmax1216 then DTxnP+1;if
tmax1316 then DTxnP+1;if tmax1416 then
DTxnP+1;
if hrd1>65 and pr1>0.2 then do;DPrHR+1;end;if hrd2>65 and pr2>0.2 then
do; DPrHR+1;end;if hrd3>65 and pr3>0.2 then do;DPrHR+1;end;
if hrd4>65 and pr4>0.2 then do;DPrHR+1;end;if hrd5>65 and pr5>0.2 then
do; DPrHR+1;end;if hrd6>65 and pr6>0.2 then do;DPrHR+1;end;
if hrd7>65 and pr7>0.2 then do;DPrHR+1;end;if hrd8>65 and pr8>0.2 then
do; DPrHR+1;end;if hrd9>65 and pr9>0.2 then do;DPrHR+1;end;
173