Grant McDermott

@gmcd.bsky.social

5187 Followers 175 Following

Economics and data science | Former academic now working in the tech sector | Views my own | 🇿🇦 in 🇺🇸 | https://grantmcdermott.com/

NBC rugby pundits, Alex Corbisiero and Dan Kyle, inexplicitly in the same room rather than Skyping from their kitchens.

https://github.com/grantmcdermott/tinyplot/blob/main/vignettes/gallery_figs/polygon-valentine.R

tinyplot(iris, by = "Species", frames = TRUE)

https://github.com/grantmcdermott/tinyplot/blob/main/vignettes/gallery_figs/direct-labels-aq.R

https://github.com/grantmcdermott/tinyplot/blob/main/vignettes/gallery_figs/likert.R

https://github.com/grantmcdermott/tinyplot/blob/main/vignettes/gallery_figs/add-text-mtcars.R

https://github.com/grantmcdermott/tinyplot/blob/main/vignettes/gallery_figs/bubble-quakes.R

Grouped line (tiny)plot with direct labels.

Merged PR

tinyplot boxplot with jittered points

A nicer print.data.frame method showing column types, as well as a subset of rows. Inspired by data.table's print method.

dbreg example showing how to estimate a model against a large dataset, comprising parquet files on disk. The model takes < 1 to run, including the time needed to read the data from disk.

dbreg(
tip_amount ~ fare_amount + passenger_count | month + vendor_name,
path = "read_parquet('nyc-taxi/**/*.parquet')",
vcov = ~month, # clustered SEs
strategy = "compress" # skip auto strategy overhead
)

A simple dbbinsreg example from the dbreg package.

# Accounting for Diet "fixed effects" helps to resolve the situation
# (We'll also add a line and change the theme for a nicer plot)
dbbinsreg(weight ~ Time | Diet, data = ChickWeight, nbins = 10, cb = TRUE,
line = c(1, 1), theme = "clean")

WIP potting shed / structure.

A post in hole

It’s super impressive that agents were able to bootstrap a browser in a week! But holy crap! I hope nobody ever uses that thing or would try to build an actual browser out of it, at least with this generation of agents, it’s still pure slop with little oversight. It’s an impressive research and tech demo, not an approach to building software people should use. At least not yet

Brandolini's law: The amount of energy needed to refute bullshit is an order of magnitude bigger than that needed to produce it.

Bandit. The dad we all want to be when we grow up.

A houseplant casts a shadow of heart-shaped climbers in the background

Example of dodged pointrange (tiny)plot, with aligned ribbons layered on top.

A bloodied Damian McKenzie lining up a matching winning penalty against Scotland

Showing the minimum no. of voters that could have flipped the results from recent US presidential elections.

elec_marg(2016)
#> Winner: Republican
#> Min votes to flip Democratic: 77,744
#> 0.06% of total votes
#> 0.6% of votes from 3 states (Michigan, Pennsylvania, Wisconsin)

elec_marg(2020)
#> Winner: Democratic
#> Min votes to flip Republican: 123,473
#> 0.08% of total votes
#> 0.7% of votes from 4 states (Georgia, Arizona, Wisconsin, Pennsylvania)

elec_marg(2024)
#> Winner: Republican
#> Min votes to flip Democratic: 344,866
#> 0.2% of total votes
#> 2% of votes from 4 states (Wisconsin, Michigan, Pennsylvania, Georgia)
Comment

tinyplot showing US GDP since 1980 with recessions highlighted in the background

tinyplot bubble chart example with dual legend, using penguins data. Taken from https://grantmcdermott.com/tinyplot/vignettes/gallery.html

Smorgasbord, of sorts. (A bunch of veg.)

"Whoever designed the human visual system has a lot to answer for. Maybe not as much as the knee team, or the lower back group. But a lot." Quote from @kjhealy.co

Running dbreg::dbreg() on a full year of NYC taxi data... and it takes less than 2 seconds.

dbreg(
tip_amount ~ fare_amount + passenger_count | month + vendor_name,
path = "read_parquet('nyc-taxi/**/*.parquet')", ## path to hive-partitioned dataset
vcov = "hc1"
)
#> [dbreg] Estimating compression ratio...
#> [dbreg] Data has 178,544,324 rows and 24 unique FE groups.
#> [dbreg] Using strategy: compress
#> [dbreg] Executing compress strategy SQL
#>
#> Compressed OLS estimation, Dep. Var.: tip_amount
#> Observations.: 178,544,324 (original) | 70,782 (compressed)
#> Standard Errors: Heteroskedasticity-robust
#> Estimate Std. Error t value Pr(>|t|)
#> fare_amount 0.106744 0.000068 1564.742 < 2.2e-16 ***
#> passenger_count -0.029086 0.000106 -273.866 < 2.2e-16 ***
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Screenshot of a `tinyplot` figure displayed directly in an R console (launched from the terminal), thanks to the new `terminalgraphics` package.

Comparing the spacing of the same tinyplot without and without themes.

plt(wageinc ~ age | gender, data=svcensus, type='lm')
plt_add(type='p', alpha = 0.1)

tinyplot lm type without SE bands.

Raised garden beds

From the fixest NEWS file: the new default VCOV is iid for all estimations. To change the default to the way it was, place setFixest_vcov(all = "cluster", no_FE = "iid") in your .Rprofile.

Grouped line plot with direct labels

Screenshot of a Quarto reveal.js presentation using the "clean" theme.

library(broom)
mod = tidy(lm(Temp ~ 0 + factor(Month), data = airquality), conf.int = TRUE)

plt(
estimate ~ term, ymin = conf.low, ymax = conf.high,
data = mod,
type = "pointrange", flip = TRUE
)

plt(
Income ~ Population, data.frame(state.x77),
log = "x",
xaxl = "l", xaxb = 10^(2:5),
yaxl = "$"
)

plt(count ~ spray, data = InsectSprays, type = "violin")

library("tinyplot")
tinytheme("clean2")

plt(
weight ~ Time | Diet, ChickWeight,
type = "barplot", beside = TRUE
)

tinyplot example of flipped density plot with groups and faceting

tinyplot example showing use of plt_add() to layer OLS fit onto an existing plot

tinyplot example of a ridgline plot (type = "ridge") of monthly temperatures

tinyplot scatterplot showing earthquake locations off Fiji, using tinytheme("default")

Abstract art example (Peale's "Portrait of Rosalba") using parttree. https://grantmcdermott.com/parttree/articles/parttree-art.html

Example of a classification tree partition.

An air duct in the shape and coloring of an Arrakis sandworm.

Title sequence for Ronin (1998).

Screenshot demonstrating etwfe syntax for estimating an event study and then plotting the results.

A woman (Margaret Hamilton of NASA) stands next to a pile of code manuals as tall as herself.