# Run this code to load the required packages
suppressMessages(suppressWarnings(suppressPackageStartupMessages({
library(mosaic)
library(supernova)
library(Lock5withR)
})))
# To make slightly smaller plots
options(repr.plot.width = 5, repr.plot.height = 3)
CensusSchool <- read.csv("https://docs.google.com/spreadsheets/d/e/2PACX-1vSVaWnM4odSxy0mlnhWvvGbeLtiKoZmsbqC6KLzXtBOjQfrF9EVKuX4RVh3XbP3iw/pub?gid=2100178416&single=true&output=csv", header = TRUE)
str(CensusSchool)
'data.frame': 10113 obs. of 60 variables:
$ Country : chr "USA" "USA" "USA" "USA" ...
$ Region : chr "SD" "PA" "CO" "MN" ...
$ DataYear : int 2011 2020 2015 2018 2014 2011 2011 2011 2021 2011 ...
$ ClassGrade : int 9 9 9 9 9 9 9 9 9 9 ...
$ Gender : chr "Female" "Female" "Male" "Female" ...
$ Ageyears : num 16 14 15 14 16 14 15 14 14 14 ...
$ Handed : chr "Right-Handed" "Left-Handed" "Right-Handed" "Right-Handed" ...
$ Height_cm : num 61 170 68 16256 385 ...
$ Footlength_cm : num 19 30 24 8 22 800 22 25 25 23 ...
$ Armspan_cm : num 157 67 164 NA 147 168 150 166 172 148 ...
$ Languages_spoken : num 0 1 1 1 1 1 1 1 1 2 ...
$ Travel_to_School : chr "Bus" "Bus" "Car" "Bus" ...
$ Travel_time_to_School : num 45 45 NA 20 12 20 15 10 6 20 ...
$ Reaction_time : num 0.422 0.31 0.34 0.491 0.433 0.297 0.375 61 0.721 0.469 ...
$ Score_in_memory_game : num 44 47 41 46 48 35 37 60 44 41 ...
$ Favourite_physical_activity : chr "Swimming" "Basketball" "Baseball/Softball" "Dancing" ...
$ Importance_reducing_pollution : num 5 1000 1000 1000 900 350 25 0 5 NA ...
$ Importance_recycling_rubbish : num 400 900 1000 1000 0 350 15 0 21 1050 ...
$ Importance_conserving_water : num 600 900 1000 0 800 250 50 0 102 1000 ...
$ Importance_saving_energy : num NA 900 500 1000 1000 NA NA NA 1 NA ...
$ Importance_owning_computer : num 0 0 0 0 1000 0 0 0 5 NA ...
$ Importance_Internet_access : num 0 0 0 0 0 0 0 0 0 0 ...
$ Left_Footlength_cm : num 18 30 24 9 21 800 21 26 24 23 ...
$ Longer_foot : chr "Right foot" "Same length" "Same length" "Same length" ...
$ Index_Fingerlength_mm : num 10 8 90 3 60 3000 6 8 7 8 ...
$ Ring_Fingerlength_mm : num 6 7 95 5 60 3500 6 10 6.5 7 ...
$ Longer_Finger_Lefthand : chr "Index finger" "Index finger" "Index finger" "Ring finger" ...
$ Birth_month : chr "October" "June" "August" "April" ...
$ Favorite_Season : chr "Spring" "Summer" "Summer" "Fall" ...
$ Allergies : chr "Yes" "Yes" "No" "Yes" ...
$ Vegetarian : chr "No" "No" "No" "No" ...
$ Favorite_Food : chr "Pizza/Pasta" "Desserts" "Meat" "No favorite" ...
$ Beverage : chr "Water" "Water" "Water" "Water" ...
$ Favorite_School_Subject : chr "History" "Other" "Mathematics and statistics" "History" ...
$ Sleep_Hours_Schoolnight : num 8 8 8 8 9 7 8 9 8 5 ...
$ Sleep_Hours_Non_Schoolnight : num 8 9 8 7 8 9 11 7 NA 13 ...
$ Home_Occupants : num 5 4 5 6 4 5 2 4 8 4 ...
$ Home_Internet_Access : chr "No internet connection" "Yes - broadband connection" "Yes - other" "Yes - other" ...
$ Communication_With_Friends : chr "In person" "Text messaging" "Text messaging" "Internet chat or instant messaging" ...
$ Text_Messages_Sent_Yesterday : num 0 10 3 8 200 500 12 250 1 249 ...
$ Text_Messages_Received_Yesterday: num 0 10 4 3 500 500 10 300 0 305 ...
$ Hanging_Out_With_Friends_Hours : num 2 NA 3 5 0 10 14 20 12 3 ...
$ Talking_On_Phone_Hours : num 0 1 2 7 10 0 0 15 3 1 ...
$ Doing_Homework_Hours : num 3 10 3 3 100 2 28 10 6 1 ...
$ Doing_Things_With_Family_Hours : num 1 4 15 7 10 9 28 15 12 0 ...
$ Outdoor_Activities_Hours : num 4 14 15 5 0 5 0 50 7 0 ...
$ Video_Games_Hours : num 0 0 3 0 6 0 25 0 0 2 ...
$ Social_Websites_Hours : num 1 21 3 7 10 0 7 7 0 3 ...
$ Texting_Messaging_Hours : num 1 5 2 7 56 10 24 10 0 2 ...
$ Computer_Use_Hours : num 1 25 2 5 6 1 24 5 0 12 ...
$ Watching_TV_Hours : num 2 18 5 3 1 1 28 0 0 2 ...
$ Paid_Work_Hours : num 5 2 0 0 0 0 0 0 0 3 ...
$ Work_At_Home_Hours : num 2 2 10 3 10 6 24 10 7 4 ...
$ Schoolwork_Pressure : chr "Some" "A lot" "Some" "Very little" ...
$ Planned_Education_Level : chr "Some college" "Undergraduate degree" "Some college" "Some college" ...
$ Favorite_Music : chr "Punk rock" "Pop" "Rap/Hip hop" "Rap/Hip hop" ...
$ Superpower : chr "Invisibility" "Fly" "Invisibility" "Telepathy" ...
$ Preferred_Status : chr "Healthy" "Happy" "Famous" "Healthy" ...
$ Role_Model_Type : chr "Friend" "Relative" "Sports person" "Actor" ...
$ Charity_Donation : chr "Health" "Health" "Health" "International aid" ...
select(CensusSchool, Sleep_Hours_Schoolnight, Texting_Messaging_Hours)
CensusSchool$Texting_Messaging_Hours
gf_jitter(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = CensusSchool)
Warning message:
“Removed 197 rows containing missing values (geom_point).”
newHours <- filter(CensusSchool, Texting_Messaging_Hours < 24*7)
NewSchoolHours <- filter(newHours, Sleep_Hours_Schoolnight < 12)
NewSchoolHours1 <- filter(NewSchoolHours, Sleep_Hours_Schoolnight >= 1)
Texting_Sleep <- filter(NewSchoolHours1, Sleep_Hours_Schoolnight != 0)
gf_bar(~ Texting_Messaging_Hours, data = Texting_Sleep, width = 20)
gf_bar(~ Sleep_Hours_Schoolnight, data = Texting_Sleep, width = 6)
Warning message:
“position_stack requires non-overlapping x intervals”
Warning message:
“position_stack requires non-overlapping x intervals”
gf_jitter(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep, alpha = 0.1)
gf_jitter(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep, alpha = 0.2, color = "goldenrod")
gf_jitter(Sleep_Hours_Schoolnight ~ shuffle(Texting_Messaging_Hours), data = Texting_Sleep, alpha = 0.2, color = "darkorchid1")
sleep_empty <- lm(Sleep_Hours_Schoolnight ~ NULL, data = Texting_Sleep)
sleep_empty
sleep_model <- lm(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep)
sleep_model
gf_jitter(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep, color = "magenta", alpha = 1) %>%
gf_hline(yintercept = 6.80, color = "orange")
Warning message:
“geom_hline(): Ignoring `mapping` because `yintercept` was provided.”
Sleep_pred <- predict(sleep_model)
gf_jitter(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep, color = "magenta", alpha = 1) %>%
gf_jitter(Sleep_pred ~ Texting_Messaging_Hours, color = "blue") %>%
gf_hline(yintercept = 6.80, color = "orange")
Warning message:
“geom_hline(): Ignoring `mapping` because `yintercept` was provided.”
supernova(sleep_model)
sampleF <- fVal(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep)
sdof <- do(1000) * fVal(Sleep_Hours_Schoolnight ~ shuffle(Texting_Messaging_Hours), data = Texting_Sleep)
gf_histogram(~ fVal, data = sdof)
gf_histogram(~ fVal, data = sdof,fill = ~fVal>sampleF) %>%
gf_point(0 ~ fVal(Sleep_Hours_Schoolnight ~ Texting_Messaging_Hours, data = Texting_Sleep))
confint(sleep_model)