Hello! This is my thesis project. I want to know the influence of THI (temperature and humidity index) on the EIPH (excercise induced pulmonary haemorrhage) in horses. We did an average of THI for each month. I want to do the same relation for each part (in the program) but separating it by year and by month. This analysis part 1: I did the relation with THI and Positive/negatives to bleeding and part 2: I divided the grades of eiph into slightly/serious bleeding. This may help you out to understand it. I tried many ways but it always came with an error in the model. I just took the class, so I only know the basics. If you can help me out here, I'll appreciate it!
FILENAME REFFILE '/folders/myfolders/sasuser.v94/DatosFinales_ITH.xlsx';
PROC IMPORT DATAFILE=REFFILE DBMS=XLSX OUT=WORK.IMPORT;
GETNAMES=YES;
RUN;
PROC CONTENTS DATA=WORK.IMPORT;
proc print;
RUN;
libname eiph '/folders/myfolders/sasuser.v94/';
data eiph.import;
set import;
run;
proc sgplot data=eiph.import pctlevel=group;
vbar thi /stat=percent group=eiph grouporder=data;
title "Proporciones por grado de EIPH según el ITH";
proc sgplot data=eiph.import;
vline mes/ response=thi stat=mean markers group=year;
title 'VALORES DE ITH POR MES Y POR AÑO';
Proc freq data=eiph.import;
Tables thi*eiph/nocol nopercent chisq;
Run;
data eiphpos;
set IMPORT;
if eiph=0 then
eiphpos=0;
if eiph>0 then
eiphpos=1;
proc glimmix data=eiphpos;
class thi;
model eiphpos(descending)=THI / dist=binary link=logit solution;
output out=grado pred(ilink)=eiphp;
lsmeans thi/ pdiff lines ilink;
proc sgplot data=eiphpos pctlevel=group;
vbar thi /stat=percent group=eiphpos grouporder=data;
Proc freq data=eiphpos;
Tables thi*eiphpos/nocol nopercent chisq;
Run;
proc sort data=grado;
by thi eiphp;
proc sgplot data=grado;
series x=THI y=eiphp;
scatter x=THI y=eiphpos / jitter;
run;
data eiph.grad;
set IMPORT;
if eiph>=4 then
eiphgrad=1;
if eiph<4 then
eiphgrad=0;
if eiph=0 then
eiphgrad=.;
proc glimmix data=eiph.grad;
class thi;
model eiphgrad(descending)=THI / dist=binary link=logit solution;
output out=gradog pred(ilink)=eiphg;
lsmeans thi / pdiff lines ilink;
proc sgplot data=eiph.grad pctlevel=group;
vbar thi /stat=percent group=eiphgrad grouporder=data;
Proc freq data=eiph.grad;
Tables thi*eiphgrad/nocol nopercent chisq;
Run;
proc sort data=gradog;
by thi eiphg;
proc sgplot data=gradog;
series x=THI y=eiphg;
scatter x=THI y=eiphgrad / jitter;
run;
... View more
