In UCLATALL/coursekata-r: Packages and Functions for CourseKata Courses

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.path = "man/figures/README-",
  out.width = "100%"
)

coursekata

Overview

CourseKata Statistics and Data Science, is an innovative interactive online textbook for teaching introductory statistics and data science in colleges, universities, and high schools. Part of CourseKata’s Better Book Project, we are leveraging research and student data to guide continuous improvement of online learning resources. The coursekata package is designed to make it easy to install and load the packages, functions, and data used in the book and supplementary materials.

Learn more about CourseKata and its free services and materials at CourseKata.org.

This package makes it easy to install and load all packages and functions used in CourseKata courses. It additionally provides a handful of helper functions and augments some generic functions to provide cohesion between the network of packages. This package was inspired by the tidyverse meta-package.

Installation

# Install the development version from GitHub
# install.packages("remotes")
remotes::install_github("coursekata/coursekata-r")

Note that installing the package will install all of the functions that are used in the course, but by default a couple of packages will not be installed: fivethirtyeight and fivethirtyeightdata. These packages only contain data, so the R package building process complains when functions are not imported from them. The first time you call library(coursekata) you will be prompted to install the packages if they are not already installed.

Loading Packages Used in CourseKata Courses

library(coursekata) will load these core packages in addition to the functions and theme included in the coursekata package:

library(coursekata)

supernova, for
creating ANOVA tables.
tools for extracting information from fitted models (b0(), b1(), PRE(), fVal())
an augmented print.lm() which prints the fitted equation as well
… and more!
mosaic, for a unified interface to most statistical tools.
ggformula, for a formula interface to ggplot2.
dplyr, for data manipulation.
Metrics, for model evaluation.

In addition to useful functions, a great deal of data sets are used by instructors who teach the course. This package installs these:

Functions and Theme

This package also comes with a variety of functions useful for teaching statistics and data science, and it has an automatically set ggplot2 theme complete with colorblind-friendly palettes and other improvements to aid perception and clarity of plots.

Estimate Extraction (and Bootstrapping)

Extracting an estimate is as easy as passing a fitted linear model to one of the extraction functions:

fit <- lm(mpg ~ hp, data = mtcars)

# the estimate for β₀, the intercept
b0(fit)

# the estimate for β₁, the slope
b1(fit)

# the proportional reduction in error
PRE(fit)

# Fisher's F value
fVal(fit)

# the sum of the squared errors from the model
SSE(fit)

# the sum of squares model/regression
SSM(fit)

The estimate extraction functions help to simplify the ability to create bootstrapped sampling distributions of those estimates. Here is an example of bootstrapping the slope:

# use mosaic package to repetitively resample to bootstrap a distribution
samp_dist_of_b1 <- do(1000) * b1(lm(mpg ~ hp, data = resample(mtcars)))

# plot the bootstrapped estimates
gf_histogram(~ samp_dist_of_b1$b1)

Other estimates and terms can be bootstrapped using the same technique, but you will need to calculate the values yourself. Here's an example of doing that for a term that doesn't have a dedicated extraction function:

samp_dist_of_hp <- do(1000) * {
  # create a new model from the resampled data
  model <- lm(mpg ~ disp * hp, data = resample(mtcars))

  # extract the desired estimate, here the coefficient for hp
  coef(model)[["hp"]]
}

# plot the bootstrapped estimates
gf_histogram(~ samp_dist_of_hp$result)

Sectioning a Distribution

When teaching about hypothesis testing, F, and p-value, it is useful to mark different portions of a distribution as inside or outside the critical zone. middle(), upper(), and lower() each take a distribution of values and return whether the value was in, e.g. the middle 95% of the distribution. Use this with a plotting function to shade in those areas:

# shade in the middle 80% of the Thumb distribution
gf_histogram(~Thumb, data = Fingers, fill = ~ middle(Thumb, .80))

Toggling the Theme

The ggplot2 theme is loaded by default but can be toggled on and off via load_coursekata_themes() and restore_default_themes(). The actual plot theme and scale components are also provided for advanced users as theme_coursekata() and scale_discrete_coursekata() (viridis is used for continuous color scales).

Contributing

If you see an issue, problem, or improvement that you think we should know about, or you think would fit with this package, please let us know on our issues page. Alternatively, if you are up for a little coding of your own, submit a pull request:

Fork it!
Create your feature branch: git checkout -b my-new-feature
Commit your changes: git commit -am 'Add some feature'
Push to the branch: git push origin my-new-feature
Submit a pull request :D