Standing on one foot!

Load libraries

library(mosaic)
library(dplyr)
library(googledrive)
library(googlesheets4)

Read in data

  sheetnm <- "https://docs.google.com/spreadsheets/d/1gRTJeV0oAvaxxQOtesHszOUrHzh852c3mO0yWJnI8w4/edit?usp=sharing"
  tippynew <- read_sheet(sheetnm)

✔ Reading from "StandingTall2026 (Responses)".

✔ Range 'Form Responses 1'.

  names(tippynew)<-c("Timestamp", "Height", "Time")
  tippynew<- mutate(tippynew, year = as.character(substr(Timestamp, 1, 4)))
  tippynew <- filter(tippynew, year=="2026")
  tippynew$Height2<-"Taller"
  tippynew$Height2[tippynew$Height=="I am shorter than 5 foot 9"]<-"Smaller"
  tippynew$Height2[tippynew$Height=="I have no idea how tall I am"]<-"NotSure!"
  tippynew <- tippynew %>% select(Height, Time, Height2, year)

Read in data from past years and combine

  tippyold<-read.csv( file="data/TippyAll2025.csv") 
  tippy<-rbind(tippynew, tippyold)

Plot

  gf_boxplot(Time~Height2|year, data=tippy)

  gf_histogram(~Time|year, data=tippy)

`stat_bin()` using `bins = 30`. Pick better value `binwidth`.

Plot

  gf_density(~Time|year, fill=~Height, data=tippy)

Warning: Groups with fewer than two data points have been dropped.
Groups with fewer than two data points have been dropped.

Warning: Removed 2 rows containing missing values or values outside the scale range
(`geom_area()`).

  favstats(~Time|year, data=tippy)  

  year  min      Q1 median     Q3   max     mean       sd  n missing
1 2018 2.95  7.3900 12.945 24.140 329.0 27.66531 57.65088 32       0
2 2019 2.40  5.8500 10.300 27.865  33.0 15.80593 11.13617 27       0
3 2020 3.39 11.0000 16.000 39.000 380.0 37.39600 63.54834 85       0
4 2023 1.00  7.5125 14.655 23.250  98.4 22.14250 22.56505 44       0
5 2024 2.00  5.3500 16.660 31.200  85.0 21.79649 20.72853 37       0
6 2025 1.00  5.0000 18.000 41.000  77.0 23.46054 21.75645 37       0
7 2026 3.00 15.3750 25.525 37.100 157.0 32.49294 30.36879 34       0

Write out the data

   write.csv(tippy, file = "data/TippyAll2026.csv", row.names= FALSE)