Displaying 2 results from an estimated 2 matches for "pmva".
Did you mean:
mva
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data?
...ctmi<-rep(0,sum(air$nctrl))
capm<-rep(pm10, mi)
ctpm<-rep(pm10, air$nctrl)
cast<-rep(air$stratn, mi)
ctst<-rep(air$stratn, air$nctrl)
cady<-rep(dates, mi)
ctdy<-rep(dates, air$nctrl)
cases<-c(cami, ctmi)
stranu<-c(cast, ctst)
days<-c(cady, ctdy)
pmva<-c(capm, ctpm)
air2<-data.frame(cases, days, pmva,stranu)
air.cl<-clogit(cases~pmva + strata(stranu), method="approximate",
data=air2)
air.glm<-glm(mi~pm10, family=poisson, data=air)
#air.gam<-gam(mi~s(pm10), family=poisson)
qres[q]<-as.numeric(coef(air.cl))
glm...
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data? Unformatted text-version, with an additional note
...<-rep(1,sum(mi))
ctmi<-rep(0,sum(air$nctrl))
capm<-rep(pm10, mi)
ctpm<-rep(pm10, air$nctrl)
cast<-rep(air$stratn, mi)
ctst<-rep(air$stratn, air$nctrl)
cady<-rep(dates, mi)
ctdy<-rep(dates, air$nctrl)
cases<-c(cami, ctmi)
stranu<-c(cast, ctst)
days<-c(cady, ctdy)
pmva<-c(capm, ctpm)
air2<-data.frame(cases, days, pmva,stranu)
air.cl<-clogit(cases~pmva + strata(stranu), method="approximate", data=air2)
air.glm<-glm(mi~pm10, family=poisson, data=air)
#air.gam<-gam(mi~s(pm10), family=poisson)
qres[q]<-as.numeric(coef(air.cl))
glmres[q]&l...