Displaying 4 results from an estimated 4 matches for "air2".
Did you mean:
air
2012 Nov 25
1
Error : Error in if (antipodal(p1, p2))
...ot;#303030", fill=TRUE, bg="black", lwd=0.05, xlim=xlim,
ylim=ylim)
fsub <- flights[flights$airline == "AA",]
fsub <- fsub[order(fsub$cnt),]
maxcnt <- max(fsub$cnt)
for (j in 1:length(fsub$airline))
{
air1 <- airports[airports$iata == fsub[j,]$airport1,]
air2 <- airports[airports$iata == fsub[j,]$airport2,]
inter <- gcIntermediate(c(air1[1,]$long, air1[1,]$lat), c(air2[1,]$long,
air2[1,]$lat), n=100, addStartEnd=TRUE)
colindex <- round( (fsub[j,]$cnt / maxcnt) * length(colors) )
lines(inter, col=colors[colindex], lwd=0.8)
}
--
View thi...
2013 Feb 09
1
R maps Help
...ot;#00FF00", "#FF0000"))
colors <- pal(100)
fsub <- flights[flights$airline == "FSR",]
fsub <- flights[order(flights$cnt),]
maxcnt <- max(fsub$cnt)
for (j in 1:length(fsub$airline)) {
air1 <- airports[airports$iata == fsub[j,]$airport1,]
air2 <- airports[airports$iata == fsub[j,]$airport2,]
p1 <- c(air1[1,]$long, air1[1,]$lat)
p2 <- c(air2[1,]$long, air2[1,]$lat)
inter <- clean.Inter(p1,p2,n=100, addStartEnd=TRUE)
colindex <- round( (fsub[j,]$cnt / maxcnt) * length(colors))
lines(inter, col=colors[co...
2008 Mar 07
0
How to do a time-stratified case-crossover analysis for air pollution data?
...trl))
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]<-as.numeric(coef(...
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(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]<-as.numeric(coef(air.g...