Posts

It is with great pleasure that I can announce that Ryacas version 1.1.0 has now been accepted into Journal of Open Source Software and …
Last updated on 2 min read R

Back in 2017, David Balding and I published the paper “How convincing is a matching Y-chromosome profile?”. One of the key parameters …
Last updated on 5 min read Forensic genetics

In inexact line search (a numerical optimisation technique) the step length (or learning rate) must be decided. In connection to that …
Last updated on 6 min read Statistics

It is with great pleasure that I can announce that Ryacas version 1.0.0 is now released to CRAN …
Last updated on 4 min read R, R-bloggers

Update Oct 14, 2019: Michael Höhle caught a mistake and notified me on Twitter. Thanks! The problem is that I used \(\text{Unif}(-10, …
Last updated on 4 min read Statistics, R-bloggers

In a recent blog post I tried to get yacas to solve a system of polynomial equations. Unfortunately it could not do that, so I solved …
Last updated on 3 min read R, R-bloggers

Finding prediction intervals (for future observations) is something different than finding confidence intervals (for unknown population …
Last updated on 4 min read R, R-bloggers

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 - …
Last updated on 3 min read R

Let \(\rho_{XY}\) be the correlation between the stochastic variables \(X\) and \(Y\) and similarly for \(\rho_{XZ}\) and …
Last updated on 6 min read Statistics, R, R-bloggers

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 …
Last updated on 3 min read R, R-bloggers