Statistical routines in Julia