Statistics, computation

…and all that jazz

Recent Posts

In a recent blog post I used something like this (for use in a call to optim()): obj_fun_expr <- expression((x^4 + 4 * x^2 * y^2 - …

Let \(\rho_{XY}\) be the correlation between the stochastic variables \(X\) and \(Y\) and similarly for \(\rho_{XZ}\) and …

It is important to test software. One approach is unit-testing, and for R packages this can e.g. be done using testthat. It is also …

Update Aug 10, 2019: I wrote a new blog post about the same as below but using a simulation approach. Let \(\rho_{XY}\) be the …

As you may know, I am maintaining the Ryacas package (with online documentation) for doing symbolic mathematics (and other stuff) in R …