Displaying 1 result from an estimated 1 matches for "sort_ev".
Did you mean:
sortrev
2005 Jul 28
1
conversion from SAS
...3.14*dayno/365);
sin1=sin(2*3.14*dayno/365);
cos2=cos(4*3.14*dayno/365);
sin2=sin(4*3.14*dayno/365);
cos3=cos(6*3.14*dayno/365);
sin3=sin(6*3.14*dayno/365);
cos4=cos(8*3.14*dayno/365);
sin4=sin(8*3.14*dayno/365);
bloom=0;
w_chla=1/chla/chla;
run;
ODS listing close;
%macro sort_event(cut_off,last=0);
/*proc glm data=sort_dataset;
class year;
model logchla=year cos1 sin1 cos2 sin2 cos3 sin3
cos4 sin4 /solution;
by station;
where bloom=0;
output out=chla_res predicted=pred student=studres
cookd=cookd rstudent=rstudent u95=u95;
lsmeans year / at (cos1 sin1 cos...