tj mahr 🎃
@tjmahr.com
data scientist studying how kids learn to speak, dad, jump roper, bayesian, tjmahr.com
yes that’s what happened today www.axios.com/2025/11/12/n...
November 13, 2025 at 12:14 AM
yes that’s what happened today www.axios.com/2025/11/12/n...
library(ggplot2)
library(mgcv)
ggplot(mtcars) +
aes(x = hp, y = disp) +
geom_point() +
stat_smooth(
aes(color = "vanilla gam"),
method = "gam"
) +
stat_smooth(
aes(color = "monotonic constraint"),
method = mgcv::scasm,
formula = y ~ s(x, bs = "sc", xt = c("m+"))
)
library(mgcv)
ggplot(mtcars) +
aes(x = hp, y = disp) +
geom_point() +
stat_smooth(
aes(color = "vanilla gam"),
method = "gam"
) +
stat_smooth(
aes(color = "monotonic constraint"),
method = mgcv::scasm,
formula = y ~ s(x, bs = "sc", xt = c("m+"))
)
November 12, 2025 at 8:43 PM
library(ggplot2)
library(mgcv)
ggplot(mtcars) +
aes(x = hp, y = disp) +
geom_point() +
stat_smooth(
aes(color = "vanilla gam"),
method = "gam"
) +
stat_smooth(
aes(color = "monotonic constraint"),
method = mgcv::scasm,
formula = y ~ s(x, bs = "sc", xt = c("m+"))
)
library(mgcv)
ggplot(mtcars) +
aes(x = hp, y = disp) +
geom_point() +
stat_smooth(
aes(color = "vanilla gam"),
method = "gam"
) +
stat_smooth(
aes(color = "monotonic constraint"),
method = mgcv::scasm,
formula = y ~ s(x, bs = "sc", xt = c("m+"))
)
oh come on i wasnt even doing anything weird this time
November 3, 2025 at 7:36 PM
oh come on i wasnt even doing anything weird this time
i can never remember how pROC::roc() does things so my wrapper includes other options and wordy errors
![validate_roc_direction <- function(direction, auto_allowed) {
case_low <- c(">", "case-low", "control-high")
case_high <- c("<", "case-high", "control-low")
auto <- if (auto_allowed) c("auto") else character(0)
valid_directions <- c(case_low, case_high, auto)
bad_direction <- !is.element(direction, valid_directions)
if (is.null(direction) || isTRUE(bad_direction)) {
style_or <- list("vec-last" = ", or ", "vec-sep2" = " or ")
ch <- cli::cli_vec(case_high, style = style_or)
cl <- cli::cli_vec(case_low, style = style_or)
au <- cli::cli_vec(auto, style = style_or)
bullet3 <- if (auto_allowed) {
c("*" = "To guess the direction, use {.val {auto}}")
} else {
character(0)
}
cli::cli_abort(c(
"Provide a valid {.arg direction}",
"*" = "For control < case, use {.val {ch}}",
"*" = "For control > case, use {.val {cl}}",
bullet3
))
}
if (direction %in% case_low) {
direction <- case_low[1]
} else if (direction %in% case_high) {
direction <- case_high[1]
}
direction
}](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:tpgnsvgqsu5b752uud3wmtja/bafkreidbnz66xnzb7qbss3zvuqijbyxkq3t265uuoffnjifwut5tu5czo4@jpeg)
October 30, 2025 at 4:40 PM
i can never remember how pROC::roc() does things so my wrapper includes other options and wordy errors
maybe this will make it stop giving me unsolicited rewordings
October 27, 2025 at 4:59 PM
maybe this will make it stop giving me unsolicited rewordings
the upgraded hat was crocheted by @amandawalrus.bsky.social
October 27, 2025 at 3:48 AM
the upgraded hat was crocheted by @amandawalrus.bsky.social
i don't want to mean and i lapse into being didactic in peer reviews, so here's some free stats advice
October 24, 2025 at 3:18 PM
i don't want to mean and i lapse into being didactic in peer reviews, so here's some free stats advice
figuring out what numbers praat can/not parse, and "+.0" is a no
![dialog box:
Fatal error
Praat will crash. Notify the author ([email]) with the following information: "Assertion failed in the file "TextGridEditor.cpp" at line 76: isdefined(t)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:tpgnsvgqsu5b752uud3wmtja/bafkreib56edqkkxirq3gq3dtx5oq3ta6mtiysge4iw7k4vzfflebbziqoi@jpeg)
October 14, 2025 at 2:11 PM
figuring out what numbers praat can/not parse, and "+.0" is a no
and there it is. pointsize scales with font size
patchwork::wrap_plots(
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(),
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(19)
)
patchwork::wrap_plots(
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(),
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(19)
)
October 9, 2025 at 8:29 PM
and there it is. pointsize scales with font size
patchwork::wrap_plots(
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(),
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(19)
)
patchwork::wrap_plots(
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(),
ggplot(mtcars) +
geom_point(aes(wt, mpg)) +
theme_grey(19)
)
Which of the three layers of buttons should I use?
October 9, 2025 at 3:25 AM
Which of the three layers of buttons should I use?
i wrote this windows .bat file like 15 years ago and i still pull it up, change the source and destination, and save a new copy when i need to mirror data from network drives onto my local machine
October 2, 2025 at 8:08 PM
i wrote this windows .bat file like 15 years ago and i still pull it up, change the source and destination, and save a new copy when i need to mirror data from network drives onto my local machine
lol what the heck culvers
October 1, 2025 at 7:21 PM
lol what the heck culvers
anatomy of a grift, from the video
September 30, 2025 at 11:52 PM
anatomy of a grift, from the video
They should have stuck with the Hello Kitty x Ninja Turtles figures
September 24, 2025 at 12:55 AM
They should have stuck with the Hello Kitty x Ninja Turtles figures
man, those old-school programmers...
![a 46 phoneme x 46 diacritic compatibility table expressed as a set of 46 strings where each character marks the compatibility for that phoneme and diacritic. there are some elaborate code comments to mark the rows and columns of the table, include vertically write column headers.
first few rows
// table of diacritic compatibility with the various phonemes
// . = allowed; X = banned; s = syllabic; f = not initial
// x = banned--automatically delete when reading pepper transcript file
// c = banned--automatically change when reading pepper transcript file
const char* DiacUse[DIAC_COUNT+1] = {
"...............................................", // 0. Unsure (circled)
// td
// n szst d z o @ ae @ c 34
// mbjnwdphtgkgfvhhhhszhlrh?ei@uoUI4I^eIEUUIc^^a3D // ----------NASALITY----------
"x..x.....x.....................................", // 1. Nasalization
"...............................................", // 2. Nasal Emission
".XX.XXXXX.XXXXXXXXXXXXXXX......................", // 3. Denasalization
// td
// n szst d z o @ ae @ c 34
// mbjnwdphtgkgfvhhhhszhlrh?ei@uoUI4I^eIEUUIc^^a3D // -----------LIP--------------
"XXXXXXXXXXXXXXXXXXXXXXXXX...ccc.......cccccc...", // 4. Rounded Vowel
"XXXXXXXXXXXXXXXXXXXXXXXXXXXX...XXXXXXX......XXX", // 5. Unrounded Vowel
"..X......................XXXXXXXXXXXXXXXXXXXXXX", // 6. Labialized Consonant
"XXXX.XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX", // 7. Nonlabialized Consonant
"...............................................", // 8. Inverted](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:tpgnsvgqsu5b752uud3wmtja/bafkreie4im562tixqozl5yrzhr4q5jarlurn4dxm3kv4hzec27o37rztge@jpeg)
September 22, 2025 at 7:38 PM
man, those old-school programmers...
Costco has toy warehouses for cats
September 19, 2025 at 10:50 PM
Costco has toy warehouses for cats
here’s my YouTube history from a couple days ago lol
September 19, 2025 at 2:41 AM
here’s my YouTube history from a couple days ago lol
overwrite duckdb's show() method when it's loaded
![my_custom_print <- function(object) {
version <- DBI::dbGetInfo(object)[["db.version"]]
tables <- DBI::dbListTables(object)
cli::cli_bullets(
c(
"{.cls {class(object)}}",
"*" = "{.emph driver}: {duckdb:::drv_to_string(object@driver)}",
"*" = "{.emph db_version}: {.value {version}}",
"*" = "{.emph tables}: {.field {tables}}"
)
)
invisible(NULL)
}
# overwrite duckdb's show() method when it's loaded
setHook(
packageEvent("duckdb", "onLoad"),
function(...) setMethod("show", "duckdb_connection", my_custom_print),
action = "append"
)](https://cdn.bsky.app/img/feed_thumbnail/plain/did:plc:tpgnsvgqsu5b752uud3wmtja/bafkreidabxbxthfyqydd224wslqabldzacw2ym4aavkzwc5p3uuom5x2h4@jpeg)
September 17, 2025 at 3:58 PM
overwrite duckdb's show() method when it's loaded
weirdest lil guy in Silksong so far
September 15, 2025 at 12:41 AM
weirdest lil guy in Silksong so far
