Displaying 2 results from an estimated 2 matches for "glmres".
Did you mean:
glares
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data?
...s)
moeff<-.5*sin(2*pi*moeff/365 + 127)
#new version; to have strictly stationary levels of air-pollution
#within strata
pm10<-2 + .5*moeff[tymofa] + 0.1*ar1
#intercept and proportionality to pm10 coefficients of log-rate
al<- -2.5
be<- 1.4
qres<-numeric(samp)
glmres<-qres
gamres<-qres
for(q in 1:samp)
{
rate<-exp(al + be*pm10)
mi<-rpois(length(rate),lambda=rate)
date1<-as.Date(startdate)
dates<-date1 + (ti-1)
tyda<-weekdays(dates)
tymo<-months(dates)
tyyear<-as.character(dates)
for (i in 1:n)
{
ask<-tyyear[...
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data? Unformatted text-version, with an additional note
...atch(tymo,monthnames)
moeff<-.5*sin(2*pi*moeff/365 + 127)
#new version; to have strictly stationary levels of air-pollution
#within strata
pm10<-2 + .5*moeff[tymofa] + 0.1*ar1
#intercept and proportionality to pm10 coefficients of log-rate
al<- -2.5
be<- 1.4
qres<-numeric(samp)
glmres<-qres
gamres<-qres
for(q in 1:samp)
{
rate<-exp(al + be*pm10)
mi<-rpois(length(rate),lambda=rate)
date1<-as.Date(startdate)
dates<-date1 + (ti-1)
tyda<-weekdays(dates)
tymo<-months(dates)
tyyear<-as.character(dates)
for (i in 1:n)
{
? ask<-tyyear[i]
? tyyear[i]<-...