Closed nicholasjhorton closed 1 year ago
I have some speculations based on some preliminary analyses. I'm hoping that we can discuss today.
@nicholasjhorton will ask Eileen about this when we visit the Holyoke Public Library History Room.
Can this question be woven into the #24 work?
I'm still not 100% certain what's happening with the currency, but it looks like the units are 1/6 of a dollar? Note that $R^2$ is nearly one:
suppressPackageStartupMessages(library(tidyverse))
suppressPackageStartupMessages(library(mosaic))
register <- readr::read_csv(
"https://nhorton.people.amherst.edu/stat210/Parsons_page253.csv") |>
mutate(daily_pay = total / number_days)
#> Rows: 44 Columns: 12
#> ── Column specification ────────────────────────────────────────────────────────
#> Delimiter: ","
#> chr (4): lastname, firstname, signature_x, note
#> dbl (6): number_days, price1, price2, total, rent_or_board, balance_due
#> lgl (2): number_pounds, number_reams
#>
#> ℹ Use `spec()` to retrieve the full column specification for this data.
#> ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
glimpse(register)
#> Rows: 44
#> Columns: 13
#> $ lastname <chr> "Johnson", "O'Brien", "Connor", "Kennedy", "Casey", "Gal…
#> $ firstname <chr> "William", "John W.", "William", "Martin", "James", "Tho…
#> $ number_pounds <lgl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
#> $ number_reams <lgl> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, …
#> $ number_days <dbl> 24.50, 21.50, 25.00, 23.00, 25.50, 25.00, 25.00, 27.00, …
#> $ price1 <dbl> 16, 9, 16, 9, 12, 9, 12, 9, 12, 12, 12, 10, 7, 21, 15, 1…
#> $ price2 <dbl> 6, 9, 6, 9, 0, 0, 0, 0, 0, 0, 0, 6, 6, 0, 0, 6, 9, 6, 0,…
#> $ total <dbl> 67.37, 34.95, 68.75, 37.38, 51.00, 37.50, 50.00, 40.50, …
#> $ rent_or_board <dbl> NA, 4.67, 7.25, 0.00, 0.00, 0.00, 0.00, 0.00, 7.25, 0.00…
#> $ balance_due <dbl> 67.37, 30.28, 61.50, 37.38, 51.00, 37.50, 50.00, 40.50, …
#> $ signature_x <chr> "yes", "yes", "yes", "no", "no", "no", "no", "no", "no",…
#> $ note <chr> NA, "Brown and Flood?", NA, NA, NA, NA, NA, NA, NA, NA, …
#> $ daily_pay <dbl> 2.749796, 1.625581, 2.750000, 1.625217, 2.000000, 1.5000…
model <- lm(daily_pay ~ I(price1/6) + price2, data = register)
mosaic::msummary(model)
#> Estimate Std. Error t value Pr(>|t|)
#> (Intercept) 0.025724 0.016927 1.52 0.137
#> I(price1/6) 0.983465 0.007530 130.61 < 2e-16 ***
#> price2 0.013562 0.001197 11.33 6.71e-14 ***
#>
#> Residual standard error: 0.0289 on 39 degrees of freedom
#> (2 observations deleted due to missingness)
#> Multiple R-squared: 0.9979, Adjusted R-squared: 0.9978
#> F-statistic: 9282 on 2 and 39 DF, p-value: < 2.2e-16
Created on 2023-03-02 with reprex v2.0.2
These are of the form:
9/- or 6/9 or 10/6 or similar. What do they denote?