Sam Power's Avatar

Sam Power

@spmontecarlo.bsky.social

Lecturer in Maths & Stats at Bristol. Interested in probabilistic + numerical computation, statistical modelling + inference. (he / him). Homepage: https://sites.google.com/view/sp-monte-carlo Seminar: https://sites.google.com/view/monte-carlo-semina

2,413 Followers  |  1,971 Following  |  1,332 Posts  |  Joined: 22.11.2023  |  2.2277

Latest posts by spmontecarlo.bsky.social on Bluesky

On the other hand, I feel vindicated w.r.t. a past confusion, whereby an Adaptive MCMC paper was more-or-less celebrating replacing "adaptation by stochastic approximation" with "adaptation by reinforcement learning". Even then, my heart told me that no replacement had really occurred.

04.12.2025 14:06 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

It certainly makes "highly-offline" learning seem rather quaint by comparison.

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

Reassuring to read this (www.argmin.net/p/defining-r...) presentation, since I had confused myself a bit in the past w.r.t. i) just how general the attached is (in principle), and ii) observing that "what people mean when they say they're doing RL" has (historically anyways) been rather narrower.

04.12.2025 14:03 β€” πŸ‘ 6    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

(There are indeed some solutions which are o(N), though not necessarily super simple ones; I am yet to see anything which is O(1). In any case, the answer depends a bit on how well you know Q.)

04.12.2025 08:50 β€” πŸ‘ 0    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

A fun problem, which I came across recently:

Suppose that you can generate random variates X ~ P by rejection sampling from some Q.

For a given integer N, let X_(N) denote the sample maximum of N independent draws from P.

Can you simulate realisations of X_(N) at a cost of o(N)? Or even O(1)?

04.12.2025 08:49 β€” πŸ‘ 4    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0
Post image

Quite fun

01.12.2025 04:28 β€” πŸ‘ 26    πŸ” 4    πŸ’¬ 2    πŸ“Œ 0

It's a fair point! Will keep it in mind.

28.11.2025 06:43 β€” πŸ‘ 0    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

another day volunteering at the topology museum.

27.11.2025 15:24 β€” πŸ‘ 5    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

This project was initiated during our residence at a 2024 INI programme on "Stochastic Systems for Anomalous Diffusion" (www.newton.ac.uk/event/ssd/), wherein many presented works treated specific approaches to the 'robustness' issue. It seemed an opportune time to take stock of things.

27.11.2025 10:35 β€” πŸ‘ 5    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

With Giorgos Vasdekis, we have written a manuscript - arxiv.org/abs/2511.21563 - which surveys the state of affairs within this literature, outlining signals for anticipating non-robustness, principles for improving robustness, and examples of contemporary methods which confront these issues.

27.11.2025 10:35 β€” πŸ‘ 15    πŸ” 1    πŸ’¬ 1    πŸ“Œ 0

In response to this, there have been a range of proposed MCMC strategies which aim to

i) perform acceptably when these conditions hold, but
ii) degrade gracefully when these conditions start to break down,

collectively giving rise to a burgeoning literature on 'robust MCMC'.

27.11.2025 10:35 β€” πŸ‘ 6    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

Usual MCMC algorithms are typically guaranteed to work well when used to sample from target distributions for which

i) mass is reasonably well-concentrated in the centre of the state space, and
ii) the log-density is smooth and of moderate growth.

Outside of this setting, things can go poorly.

27.11.2025 10:35 β€” πŸ‘ 33    πŸ” 6    πŸ’¬ 1    πŸ“Œ 0

I'll ask around!

26.11.2025 18:08 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image

On behalf of some friends, let me quickly advertise an event taking place in London, January 12-13 2026 (sites.google.com/view/lpd-tnn), with an overall focus on 'Geometric methods in probability'. Registration is free but required, and closes on December 1 (i.e. next Monday) - exciting stuff!

26.11.2025 10:57 β€” πŸ‘ 9    πŸ” 1    πŸ’¬ 1    πŸ“Œ 0

humans in the 1960s / potentially slightly earlier

25.11.2025 10:09 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

(The pictured set is conjecturally optimal for the problem described)

24.11.2025 13:29 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image

I mean, lol

24.11.2025 13:29 β€” πŸ‘ 22    πŸ” 4    πŸ’¬ 3    πŸ“Œ 2

I'm also slightly reminded of the Beta function, but without any particular conclusions for now.

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

I like it! It hadn't occurred to me that it is somehow 'really' a product of two (rather than three) factors, but I think that I'm on board. A fan of falling / rising factorials, in any case.

21.11.2025 18:00 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

cf. bsky.app/profile/scie...

21.11.2025 16:47 β€” πŸ‘ 1    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

Tour coming to an end, as I settle in for 5 hours of train journey! Had lots of fun talking about Random Walk Metropolis, Gradient Flows, and Skill Rating in Sports (among other chats). Slides from all talks are saved at github.com/sampower88/t....

21.11.2025 16:46 β€” πŸ‘ 7    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0
Post image Post image

πŸ‰

21.11.2025 11:18 β€” πŸ‘ 5    πŸ” 0    πŸ’¬ 2    πŸ“Œ 0
Post image

splendid

20.11.2025 22:46 β€” πŸ‘ 27    πŸ” 3    πŸ’¬ 1    πŸ“Œ 1

Super gripping (and fun!) lectures here:
youtu.be/OHDYdmuLMW0?...
"Fourier Analysis & Beyond I" - Mini-course
- Stefan Steinerberger

19.11.2025 23:52 β€” πŸ‘ 11    πŸ” 3    πŸ’¬ 1    πŸ“Œ 0

great stuff, right up my alley:

arxiv.org/abs/2511.11497
'A Recursive Theory of Variational State Estimation: The Dynamic Programming Approach'
- Filip Tronarp

19.11.2025 23:04 β€” πŸ‘ 10    πŸ” 1    πŸ’¬ 0    πŸ“Œ 1
Preview
Inference in Stochastic Epidemic Models via Multinomial Approximations We introduce a new method for inference in stochastic epidemic models which uses recursive multinomial approximations to integrate over unobserved variables and thus circumvent likelihood intractabili...

I think this - arxiv.org/abs/2006.13700 - is a good starting point, IIRC.

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

Very nice! Some of the latter part seemed natural to think of in terms of Markov kernels / channels. A friend had some recent work on epidemiological models where this conjugacy between Poisson / Multinomial and 'colouring' channels was quite computationally useful.

19.11.2025 12:02 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

Prompted by glancing at arxiv.org/abs/2511.14200.

19.11.2025 10:18 β€” πŸ‘ 2    πŸ” 0    πŸ’¬ 0    πŸ“Œ 0

One application of this would be to compare the features of different biased random walks, since the convex ordering is preserved under independent summation, and allows for the control of { variance, large deviations, etc. }.

19.11.2025 10:16 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0
Post image

However, this changes if you allow for a dilation factor. In particular, there is always some Οƒ = Οƒ(p, q) such that Z(p) is dominated by Οƒβ‹…Z(q). After a bit of book-keeping and changing variables for convenience, one can deduce the attached formula.

19.11.2025 10:16 β€” πŸ‘ 3    πŸ” 0    πŸ’¬ 1    πŸ“Œ 0

@spmontecarlo is following 20 prominent accounts