This article is now at http://stla.github.io/stlapblog/posts/KantorovichWithR.html.

In this article you will firstly see how to get rational numbers arbitrary close to \( \pi \) by performing the binary splitting algorithm with the gmp package. The binary splitting algorithm fastly calculates the partial sums of a rational hypergeometric series by manipulating only integer numbers. But these integer numbers are generally gigantic

GLPK library In the article 'Using R to compute the Kantorovich distance' I have shown how to use the cddlibb C library through R with the help of the rccd R package to compute the Kantorovich distance between two probability measures (on a finite set). In the present article I show how to do so using three different ways with Julia: GLPK: Similarl

Sample size determination for a mean Sample size determination for a mean This article explains the methodology implemented in the Shiny application availbale at http://glimmer.rstudio.com/stla/samplesize_mean/ Statement of the problem Consider a preliminary experiment \( {\cal E}_0 \), whose issue is a sample \( y_0 \) of size \( n_0=5 \) generate

A reactive sliced 3D surface response A reactive sliced 3D surface response In my previous article I showed an interactive 3D surface response fitted from a model with two continous predictors. But when there is more than two continuous predictors, since we can use only two predictors at time in the image, we can only show a surface plot depending

This article is now at http://stla.github.io/stlapblog/posts/rgl_knitr.html. - CanvasMatrix.js - testgl1snapshot.png - testgl2snapshot.png - CanvasMatrix - Copie.js

In order to make a presentation, I was wondering how to display the variance of a distribution, or the variance, of a sample on a graphic. Finally, I've found this solution: What is this “ellipse” with an arrow ? This is a picture commonly used in classical mechanics to represent the moment of inertia of a body spinning around an axis of rotati