data(AustralianElectionPolling)
xyplot(ALP ~ startDate | org,
data=AustralianElectionPolling,
layout=c(1,5),
type="b",
xlab="Start Date",
ylab="ALP")
## test for house effects
library(gam)
y <- AustralianElectionPolling$ALP/100
v <- y*(1-y)/AustralianElectionPolling$sampleSize
w <- 1/v
m1 <- gam(y ~ lo(startDate,span=1/10),
weight=w,
data=AustralianElectionPolling)
m2 <- update(m1, ~ . + org)
anova(m1,m2)Run the code above in your browser using DataLab