Mark Rieke's Avatar

Mark Rieke

@markjrieke.bsky.social

certified big nerd | big #rstats dweeb | accidental python fan | he/him thedatadiary.net

506 Followers  |  280 Following  |  773 Posts  |  Joined: 29.10.2024  |  1.9839

Latest posts by markjrieke.bsky.social on Bluesky

some enterprising nerd with more time than me could compare the CRPS from the forecast publications for a better comparison. Even better still would be to see if the outlets (rip 538) could share the actual draws & look at pr(obs | model)

17.02.2026 16:02 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image

I don't think this metric is necessarily any better than brier scores, but it is ironic that the SB model would've lost if used to inform bets on polymarket, given Silver's affiliation
dl.acm.org/doi/10.1145/...

17.02.2026 15:59 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

as a workout, Fran has the highest density of suffering / time spent moving

17.02.2026 14:42 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image

sometimes they put the yee haw right into the programming language

15.02.2026 22:20 β€” πŸ‘ 6    πŸ” 1    πŸ’¬ 1    πŸ“Œ 0

olympic bobsled should add podracing noises send tweet

15.02.2026 17:21 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

Also, Yglesias lecturing people on uncertainty when it comes to interpreting polling data is laughable

15.02.2026 17:11 β€” πŸ‘ 233    πŸ” 21    πŸ’¬ 7    πŸ“Œ 0
Someone touches grass, realized they don't like computer science anymore

Someone touches grass, realized they don't like computer science anymore

Or, as written up in SIGBOIVK last year:

15.02.2026 15:42 β€” πŸ‘ 8    πŸ” 5    πŸ’¬ 1    πŸ“Œ 0
The stupidity, the criminal vandalism, the wanton destruction of information involved in dichotomisation | Statistical Modeling, Causal Inference, and Social Science

The stupidity, the criminal vandalism, the wanton destruction of information involved in dichotomisation
statmodeling.stat.columbia.edu/2026/02/14/t...

14.02.2026 18:03 β€” πŸ‘ 11    πŸ” 4    πŸ’¬ 0    πŸ“Œ 1

I have concluded that the gamma distribution has too many parameterizations and simply needs to chill, thank you

13.02.2026 17:32 β€” πŸ‘ 4    πŸ” 1    πŸ’¬ 0    πŸ“Œ 0

after my sufficient gamma crashout last night, I've done some simulations & realized:
- the pdf below works for a sufficient gamma
- the original pdf I shared ALSO works for a sufficient statistic under a different parameterization
- weve already done this at work under a DIFFERENT parameterization

13.02.2026 17:28 β€” πŸ‘ 5    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

we have a slightly different formulation at work that gets around this --- I'll ping you on slack when I get online !

13.02.2026 12:57 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

ooh! no, wasn't aware --- aside from the super simple sufficient forms (binomial & counts of poisson obs), I haven't done too much w/sufficiency (I've implemented a sufficient normal like one time remember it being a struggle lol)

13.02.2026 04:32 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

oh yeah that's why my gut says that sum(x) ~ Gamma(n * alpha, theta) isn't a sufficient stat.

Link here shows what I get if I just repeatedly multiply the gamma pdf given new elements of x. Need to double check but it feels like this might actually give sufficiency?

bsky.app/profile/mark...

13.02.2026 04:18 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 2    πŸ“Œ 0

if I just sit down and bake out what a repeated product of the gamma pdf is I end up with this --- needs to be checked with simulation, but certainly passes the gut check

13.02.2026 04:14 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 2

I haven't yet simulated, but IIRC my thinking was that by summarizing (x1, x2, ... xn) to sum(x), you lose some information in that there are many combinations of alpha/theta that can reasonably be fit to a single observation of sum(x)

13.02.2026 04:07 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

someone also has pointed out that x ~ Gamma(alpha, theta) is equivalent to sum(x) ~ Gamma(n * alpha, theta). At one point I had convinced myself that there was a loss of individual information here, but now I can't quite see why I came to that conclusion --- gotta simulate some stuff to see !

12.02.2026 23:36 β€” πŸ‘ 4    πŸ” 0    πŸ’¬ 2    πŸ“Œ 0

oops xprod and xsum should index n here (xprod[n] & xsum[n])

12.02.2026 22:29 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

spitballing a function and maybe this works? would need to test it out against a non-sufficient gamma

12.02.2026 22:28 β€” πŸ‘ 0    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

moreso scary in the sense that I've got ptsd from thinking the sufficient normal was straightforward then banging my head against a wall until I got a working solution lol

(the sufficient normal was also my crash course in "sometimes the centered parameterization samples better")

12.02.2026 22:03 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image

pinging the #bayesian / #rstats / #pydata hivemind --- has anyone implemented a sufficient formulation of a gamma distribution/willing to share code (ideally in stan or pymc)? supposedly this equation is a density function for the sufficient gamma, but ngl it's scary

12.02.2026 21:42 β€” πŸ‘ 5    πŸ” 4    πŸ’¬ 1    πŸ“Œ 4

point estimate probability says 100% chance of bayesian takeover 😈

12.02.2026 16:02 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

it is a plot, we will get all of you, eventually, to embrace bayesian inference, resistance is futile, etc. etc.

12.02.2026 13:12 β€” πŸ‘ 4    πŸ” 0    πŸ’¬ 2    πŸ“Œ 0
butterfly meme
ME
ANYTHING
IS THIS A GAUSSIAN PROCESS

butterfly meme ME ANYTHING IS THIS A GAUSSIAN PROCESS

New gaussian process slides going well

11.02.2026 16:48 β€” πŸ‘ 95    πŸ” 8    πŸ’¬ 7    πŸ“Œ 2

I have reached the "writing yaml specifications" stage of my career

10.02.2026 17:14 β€” πŸ‘ 0    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

why is jazz so fuckin' good

08.02.2026 17:39 β€” πŸ‘ 4    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

"Is my sample large enough for Bayesian statistics?" is a weird question coming from someone who only ever uses frequentist stats that are only approximately correct for large samples.

08.02.2026 15:26 β€” πŸ‘ 12    πŸ” 1    πŸ’¬ 3    πŸ“Œ 1

"production notebook" is an oxymoron and should be shunned from pleasant society

06.02.2026 23:28 β€” πŸ‘ 1    πŸ” 1    πŸ’¬ 0    πŸ“Œ 0

this series is mostly horny nonsense (positive review), but this scene is legitimately incredible and moving

06.02.2026 17:25 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

nothing so captivating as watching someone exert complete mastery over a skill with ease

06.02.2026 13:00 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

coulda left it at "please don't email me"

05.02.2026 19:52 β€” πŸ‘ 0    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

@markjrieke is following 20 prominent accounts