search for: dayno

...L2000 REPLACE; GETNAMES=YES; RUN; data chla_italian; set chla_italian; year=year(datepart(date)); month=month(datepart(date)); run; proc 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/...

...I am using package mgcv to describe presence of a migratory bird species as a function of several variables, including year, day number (i.e. day-of-the-year), duration of survey, latitude and longitude. Thus, the "global model" is: global_model<-gam(present ~ as.factor(year) + s(dayno, k=5) + s(duration, k=5) + s(x, k=5) + s(y, k=5), family = binomial(link = logit), data = presence, gamma =1.4) The model works fine, suggesting that all the variables have strong, non-linear, influences on the probability of detecting the species. My main interest is in the effect "day...

hi, R-users I have a data.frame for example test$newdataday24 and test$newdataday48 I can plot them by plot(test$newdataday24) but now i want to plot different data by define a variable to describe them dayno<-c(24,48) newnam<-paste("test$newdataday",dayno,sep="") plot(newnam[1]) but i failed,the error message said that something wrong with plot.window what can i do to fix my script ? thanks - TANG Jie [[alternative HTML version deleted]]

Hello there fellow R-users, I have received some data which comes in the following format: example1<-"200301" The first 4 digits correspond to the year and the remaining 2 digits correspond to the week of the year. I have tried to convert this to a date by using strptime as follows: strptime(example1,format="%Y%U") where U (looking up strptime) is the week of the