Displaying 2 results from an estimated 2 matches for "tyda".
Did you mean:
tmda
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data?
...al)
n<-2*365
samp<-100
ti<-1:n
ar1<-rnorm(1)
for (i in 2:n){
ar1[i]<-0.2*ar1[i-1]+rnorm(1)
}
#old version of air pollution
#pm10<-2 + .5*sin(2*pi*ti/365 + 127) + 0.1*ar1
startdate<-"1992-07-01"
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]<-substr(ask,1,4)
}
moeff<-cumsum(c(1,31,28,31,30,31,30,31,31,30,31,30))
a<-as.Date("1994-12-31")
monthnames<-months(a+moeff)
tymofa<-m...
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data? Unformatted text-version, with an additional note
...t:
library(survival)
n<-2*365
samp<-100
ti<-1:n
ar1<-rnorm(1)
for (i in 2:n){
? ar1[i]<-0.2*ar1[i-1]+rnorm(1)
}
#old version of air pollution
#pm10<-2 + .5*sin(2*pi*ti/365 + 127) + 0.1*ar1
startdate<-"1992-07-01"
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]<-substr(ask,1,4)
}
moeff<-cumsum(c(1,31,28,31,30,31,30,31,31,30,31,30))
a<-as.Date("1994-12-31")
monthnames<-months(a+moeff)
tymofa<-match(tymo,month...