# Run this code to load the required packages
suppressMessages(suppressWarnings(suppressPackageStartupMessages({
library(mosaic)
library(supernova)
library(Lock5withR)
})))
# Adjust the plots to be a bit smaller
options(repr.plot.width = 6, repr.plot.height = 4)
CensusSchool <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vSVaWnM4odSxy0mlnhWvvGbeLtiKoZmsbqC6KLzXtBOjQfrF9EVKuX4RVh3XbP3iw/pub?gid=2100178416&single=true&output=csv", header = TRUE)
str(CensusSchool)
CensusSchool <- filter(CensusSchool, Doing_Things_With_Family_Hours <= 16*7)
CensusSchool <- filter(CensusSchool, Preferred_Status!= "NA")
gf_histogram(~Doing_Things_With_Family_Hours, data = CensusSchool)%>%
gf_facet_grid(Preferred_Status ~ .)
gf_jitter(Preferred_Status ~ Doing_Things_With_Family_Hours, data = CensusSchool)
lm(Doing_Things_With_Family_Hours ~ NULL, data = CensusSchool) -> Empty_model
predict(Empty_model) -> CensusSchool$Empty.pred
Hours.Model <- lm(Doing_Things_With_Family_Hours ~ Preferred_Status, data = CensusSchool)
predict(Hours.Model) -> CensusSchool$Empty.pred
Hours.Model
Empty_model
gf_jitter(Doing_Things_With_Family_Hours ~ Preferred_Status, data = CensusSchool) %>%
gf_jitter(predict(Hours.Model) ~ Preferred_Status, color = "dodgerblue", height = 0)
gf_jitter(Doing_Things_With_Family_Hours ~ Preferred_Status, data = CensusSchool) %>%
gf_jitter(predict(Empty_model) ~ Preferred_Status, color = "blue")%>%
gf_jitter(predict(Hours.Model) ~ Preferred_Status, color = "orange")
supernova(Hours.Model)
confint(Hours.Model)
samplePRE <- PRE(Doing_Things_With_Family_Hours ~ Preferred_Status, data = CensusSchool)
sdoPRE <- do(1000) * PRE(Doing_Things_With_Family_Hours ~ shuffle(Preferred_Status), data = CensusSchool)
gf_histogram(~PRE, data = sdoPRE)%>%
gf_point(0 ~ samplePRE, color = "darkgreen")
tally(~PRE > samplePRE, data = sdoPRE)