Displaying 2 results from an estimated 2 matches for "cos3".
Did you mean:
cos
2005 Jul 28
1
conversion from SAS
...sort data=chla_italian; by station; run;
/* Check bloom for seasonal cycle outliers */
data sort_dataset;
set chla_italian;
chla=chl_a;
dayno=date-mdy(1,1,year)+1;
cos1=cos(2*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 co...
2007 Dec 19
0
leaps
...DF, p-value: < 2.2e-16*
newmat <- cosmat[,-c(5,8,10,12,13,15,16,18,20:44)]
newmat <- cbind(newmat,sinmat[,-c(7:11,13:44)])
Regmod <- regsubsets(newmat, yy)
rs <- summary(Regmod)
which.max(rs$adjr)
*[1] 8
*rs$which[which.max(rs$adjr), ]
*(Intercept) cos1 cos2 cos3 cos4 cos6
TRUE TRUE TRUE FALSE TRUE FALSE
cos7 cos9 cos11 cos14 cos17 cos19
TRUE FALSE FALSE FALSE FALSE FALSE
sin1 sin2 sin3 sin4 sin5...