import numpy as np
from scipy.stats import bootstrap

latency_ms = np.array([
    117, 121, 113, 125,
    119, 128, 115, 122,
    120, 118, 124, 116,
])
rng = np.random.default_rng(20260625)


def mean_latency(sample, axis):
    return np.mean(sample, axis=axis)


result = bootstrap(
    (latency_ms,),
    mean_latency,
    confidence_level=0.95,
    n_resamples=9999,
    method="BCa",
    rng=rng,
)

ci = result.confidence_interval
sample_mean = np.mean(latency_ms)

print(f"sample mean: {sample_mean:.2f} ms")
print(f"95% BCa confidence interval: {ci.low:.2f} to {ci.high:.2f} ms")
print(f"bootstrap standard error: {result.standard_error:.2f} ms")
print(f"resamples: {result.bootstrap_distribution.shape[-1]}")
print(f"sample mean inside interval: {ci.low < sample_mean < ci.high}")