In ProjectMOSAIC/mosaic: Project MOSAIC Statistics and Mathematics Teaching Utilities

require(knitr)
require(mosaic)
opts_chunk$set(dev=c("pdf", "png"), fig.path = "Figures/")
options(digits=2)
trellis.par.set(theme = col.mosaic())
trellis.par.set(plot.symbol = list(cex=0.4, alpha=0.6))
trellis.par.set(
  superpose.symbol = 
    list(cex=0.4, alpha=0.6, 
         pch=c(16, 21), 
         col=c("navy", "skyblue"), 
         fill=c("transparent", "transparent"))
)
trellis.par.set(par.strip.text = list(col="transparnt"))
trellis.par.set(fontsize.text = 8)
trellis.par.set(layout.heights=list(strip=0.5))
trellis.par.set(box.dot = list(cex=0.2))
trellis.par.set(box.rectangle = list(lwd=0.6))
trellis.par.set(box.umbrella = list(lty=1, lwd=0.6))
trellis.par.set( 
  strip = strip.custom(
    style=1, 
    par.strip.text=list(col="transparent"))
)
noscales = list(draw = FALSE)

xyplot( riders ~ highT, data = Riders, 
        scales=noscales, xlab="", ylab="")

xyplot( riders ~ hi | weekday, data = Riders, 
        scales=noscales, xlab="", ylab="", 
        strip = 
          strip.custom(style=1, 
                       par.strip.text=list(col="transparent")) )

xyplot( riders ~ highT, groups = weekday, 
        data = Riders, 
        scales=noscales, xlab="", ylab="") 

bwplot( riders ~ weekday, data = Riders, 
        scales=noscales, xlab="", ylab="" 
       )

bwplot( riders ~ weekday | ntiles(clouds), data = Riders, 
        scales=noscales, xlab="", ylab="", coef = 3,
        layout = c(3,1),
        strip = 
          strip.custom(style=1, 
                       par.strip.text=list(col="transparent")) 
       )

do(3) * mean(ct ~ shuffle(wday), data = Riders)

diffmean(ct ~ wday, data = Riders)
Riders.Null <- 
  do(1000) * diffmean(ct ~ shuffle(wday), data = Riders)
histogram( ~ diffmean, data = Riders.Null, v=-80)

Riders.Boot <- 
  do(1000) * diffmean(ct ~ wday, data = resample(Riders, group=wday))
histogram( ~ diffmean, data = Riders.Boot)
confint(Riders.Boot)
confint(Riders.Boot, method = "quantile")

model <- lm(ct ~ hi + clouds, data = Riders)
lm.Boot <- do(1000) * relm(model)
head(lm.Boot, 3)
histogram(~hi, data = lm.Boot)
confint(model)
confint(lm.Boot)

ProjectMOSAIC/mosaic documentation built on Feb. 21, 2024, 2:11 a.m.