Central Limit Theorem demonstrator

Generate sample means from a base distribution.

Usage

mrm_clt_demo(
  base_distribution = "unif",
  n_samples = 1000L,
  sample_size = 30L,
  seed = 42L,
  ...
)

Arguments

base_distribution

Distribution suffix passed to r<dist> (e.g. "unif", "exp", "pois").

n_samples

Number of sample means.

sample_size

Size of each sample.

seed

RNG seed.

...

Additional parameters passed to r<dist>.

Value

data.frame with sample_index, sample_mean, z_score.

Examples

# 1000 sample means of size 30 from an exponential(1) base;
# standardised z-scores converge to N(0,1):
res <- mrm_clt_demo(
  base_distribution = "exp",
  n_samples = 1000L,
  sample_size = 30L,
  seed = 42L, rate = 1
)
summary(res$z_score)
#>     Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
#> -2.71947 -0.67222 -0.05011  0.00000  0.62201  4.70571 
# mean ~ 0, sd ~ 1