# Run this code to load the required packages
suppressMessages(suppressWarnings(suppressPackageStartupMessages({
library(mosaic)
library(supernova)
library(Lock5withR)
})))
CensusSchool <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vSVaWnM4odSxy0mlnhWvvGbeLtiKoZmsbqC6KLzXtBOjQfrF9EVKuX4RVh3XbP3iw/pub?gid=2100178416&single=true&output=csv", header = TRUE)
str(CensusSchool)
gf_bar(~Schoolwork_Pressure, data=CensusSchool, color="black", fill="light blue")
filter(CensusSchool, Home_Occupants <30 )-> Home_Occupants2
gf_histogram(~Home_Occupants, data=Home_Occupants2)%>%gf_facet_grid(Schoolwork_Pressure~.)
tally(~Home_Occupants, data= Home_Occupants2)
Home_Occupants2$DataYear<2020->Home_Occupants2$PreCovid
gf_dhistogram(~Home_Occupants, data=Home_Occupants2)%>%gf_facet_grid(Schoolwork_Pressure~PreCovid)
recode(CensusSchool$Schoolwork_Pressure, "None"= 1, "Very little"= 2, "Some"= 3, "A lot"=4) -> SWP
recode(CensusSchool$Home_Occupants, "None"= 1, "Very little"= 2, "Some"= 3, "A lot"=4)-> HO
lm(HO~SWP, data=CensusSchool)
lm(HO~SWP, data=CensusSchool)-> Census.model1
confint(Census.model1)
sdob1<-do(1000)*b1(HO~shuffle(SWP), data=CensusSchool)
sampleb1<-b1(HO~SWP, data=CensusSchool)
gf_histogram(~b1, data=sdob1)%>%
gf_point(0~sampleb1)