Pseudo-Random Variate Engine
Generates 50,000 independent identically distributed uniform random draws using the 64-bit Mersenne Twister algorithm.
Compute empirical sample moments, central moments, sample skewness, and excess kurtosis across controlled data distributions, and quantify the moment recurrence relation. Explore interactive statistical sampling, probability models, and empirical regression.
Ready to sample with Order k = 1.0.
Adjust Order k, then execute statistical trial.
Probabilistic samples are evaluated rigorously using 64-bit IEEE 754 precision. Residuals reflect ordinary least-squares regression across successive observations. Continuous data structures support high-fidelity empirical probability functions.
Higher-order moments (k ≥ 4) have large sampling variance and are volatile in small datasets. The empirical data validates the theoretical moment sequence E[X^k] = 1/(k+1) for uniform distributions, yielding an exact linear relationship y = 1.000x + 1.000 with R² = 1.000000.
Generates 50,000 independent identically distributed uniform random draws using the 64-bit Mersenne Twister algorithm.
Computes sample raw moments m_k' = (1/N) ∑ X_i^k with compensated Kahan summation to prevent numerical cancellation.
Applies binomial expansion transformations to compute central moments μ_k = E[(X - μ)^k] from raw sample moments.
Evaluates standard sample skewness g_1 = m_3 / (m_2)^(3/2) to detect directional distribution asymmetry.
Computes standardized fourth central moment g_2 = [m_4 / (m_2)^2] - 3 to measure tail heaviness relative to Gaussian mesokurtosis.
Performs linear regression of inverse expectation versus moment order to verify theoretical moment recurrence.
Compute empirical sample moments, central moments, sample skewness, and excess kurtosis across controlled data distributions, and quantify the moment recurrence relation. Standardized parameter settings: 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 order.
Set the pseudo-random generator seed to 42 and allocate a continuous floating-point memory buffer for 50,000 independent samples drawn from the standard uniform distribution U[0, 1].
Set moment order k = 1. Accumulate the arithmetic sample mean m_1' = (1/N) ∑ X_i and record the expectation value on the digital measurement dashboard.
Center all sample observations by subtracting the arithmetic mean (X_i - x̄) and evaluate the second central moment m_2, third central moment m_3, and fourth central moment m_4.
Compute Fisher-Pearson standardized sample skewness g_1 and excess kurtosis g_2, logging the resulting empirical coefficients alongside their asymptotic standard errors.
Increment the power moment order k across standardized integer steps k ∈ {1, 2, 3, 4, 5, 6}, observing that for U[0, 1], theoretical expectation E[X^k] = 1/(k+1), such that 1/E[X^k] = k + 1.
Log all six empirical trials into the digital laboratory notebook. Fit the linear regression model y = mx + b on transformed coordinates (x = k, y = 1/E[X^k]) to confirm unit slope m = 1.000.
Points are your recorded observations. The solid regression line validates Sample Moments & Kurtosis mathematical scaling.
| # | Raw Moment Order k (order) | Reciprocal Normalized Expectation (1/E[X^k]) | Moment Order k | Inverse Uniform Moment 1 / E[X^k] | Remove |
|---|
No readings yet. Execute a statistical trial and record it.