Via @johnhenry/math's Statistics module (mean, median, standard deviation, five-number summary), plus Distributions.normal/binomial/poisson/studentT/chiSquare for interval probability — every distribution factory exposes the same cdf(x) shape regardless of continuous/discrete, so P(lower ≤ X ≤ upper) is computed identically across all five. An interactive draggable-marker axis (GeoGebra's Probability Calculator UX) is a later extension.
Uses the data entered in Descriptive statistics above, in entry order.
Blue dots = raw data, colored line = smoothed. The first/last 2 point(s) at each edge are trimmed (a "same"-mode convolution boundary sample there averages against zero-padding, not real neighboring data).
= 0.682689 (CDF(1) = 0.841345, CDF(-1) = 0.158655)
Uses the data entered in Descriptive statistics above as the (first) sample.
statistic = 6.6144, df = 7.00, p = 0.000300
Reject H₀ at α = 0.05 (p = 0.0003 < 0.05)