Software, applied AI, scientific computing and statistics by Lincoln Colling.https://blog.colling.net.nz/en-NZhttps://blog.colling.net.nz/posts/bayes-factors-for-comparing-variances/https://blog.colling.net.nz/posts/bayes-factors-for-comparing-variances/Comparing a sample variance with a reference variance using Bayesplay, checking directed Bayes factors, and exploring why the result should not depend on the units of measurement.Sun, 23 Aug 2026 00:00:00 GMTstatisticsbayesrhttps://blog.colling.net.nz/posts/crusty-r/https://blog.colling.net.nz/posts/crusty-r/Why would you only want to write one language when you can write all the languages. In this short post I’ll briefly explain how, with the help of C, you can write Rust functions that you can call from R.Mon, 22 Aug 2022 00:00:00 GMTrustrchttps://blog.colling.net.nz/posts/speeding-up-bfda-part2/https://blog.colling.net.nz/posts/speeding-up-bfda-part2/In Part 1 of this series, I showed how Functional-style code allows you to separate the what and the how of your code, and how doing so opened up new ways of solving programming problems. Part 1 covered the reduce / accumulate family of higher-order functions, and in the current post I cover the apply / map family of higher-order functions. Finally, I show how changing from loops to the apply / map family of higher-order functions parallel programming becomes a trivial matter of slightly modifying a function name.Thu, 07 Jul 2022 00:00:00 GMTfunctional programmingrparallel programminghttps://blog.colling.net.nz/posts/simulating_pvalues/https://blog.colling.net.nz/posts/simulating_pvalues/The properties of p values can be difficult to understand. Therefore, one way to develop a good conceptual understanding of p values is through simulation. The document below allows you to simulate experiments and to examine how the distribution of data and p-values change.Mon, 06 Jun 2022 00:00:00 GMTstatisticsfrequentistteachinghttps://blog.colling.net.nz/posts/speeding-up-bfda-part1/https://blog.colling.net.nz/posts/speeding-up-bfda-part1/There’s recently been a bit of talk on twitter about the relative benefits of loops versus more functional-style code. Most of these discussions have centred around the relative performance benefits. However, I think this slightly misses the point. Functional-style code allows you to separate the what and the how of your code in a way that opens up new ways of thinking about problems. The end result is often faster code, but this benefit comes from writing cleaner code the encourages clearer thinking. In this post, I show how to speed up some of the simulation code in the {BFDA} R package using functional design patterns. Specifically, reductions and accumulations (aka scans).Tue, 24 May 2022 00:00:00 GMTfunctional programmingrhttps://blog.colling.net.nz/posts/adventures-in-wasm/https://blog.colling.net.nz/posts/adventures-in-wasm/Currently, robustness analysis for informed Bayesian t-tests is not implemented in JASP. So I thought I’d just build it! This is also a great excuse to play around with Rust WASM and Observable blocks in Quarto.Thu, 14 Apr 2022 00:00:00 GMTprogrammingstatisticsbayes