Introduction

val sample = doubleArrayOf(2.0, 4.0, 4.0, 5.0, 7.0, 9.0)

val summary = sample.describe()
summary.mean              // 5.1667
summary.standardDeviation // 2.4833

val normality = shapiroWilkTest(sample)
normality.pValue          // 0.8933

val fitted = NormalDistribution(mu = summary.mean, sigma = summary.standardDeviation)
fitted.cdf(6.0)           // 0.6335