sample-median-dis.R

#
# sample-median-dis.R,  6 Mar 20
# Data from:
# Example
#
# Example from:
# Evidence-based Software Engineering: based on the publicly available data
# Derek M. Jones
#
# TAG example_median example_mean example_binomial-distribution


source("ESEUR_config.r")


# Need to get this plot to fit in the margin, along with the plot before it
plot_layout(2, 1, max_height=12)
par(mar=MAR_default-c(0.8, 0, 0.0, 0))


NUM_REPLICATE=1000

discrete_sample=function(sample_size)
{
#return(sample(1:10, sample_size, replace=TRUE, dnorm(1:10, mean=5, sd=2)))
return(sample(0:9, sample_size, replace=TRUE, dbinom(0:9, size=10, prob=6.5/10)))
}


samp_mean=replicate(NUM_REPLICATE, mean(discrete_sample(30)))
samp_median=replicate(NUM_REPLICATE, median(discrete_sample(30)))

plot(table(samp_median), yaxt="n", col=point_col,
	cex.axis=1.4, cex.lab=1.4,
	xlab="Median", ylab="Samples")

# Explicitly draw x-axis to cure spurious choice of axis values.
plot(table(samp_mean), yaxt="n", col=point_col,
	xaxt="n",
	cex.axis=1.4, cex.lab=1.4,
	xlab="Mean", ylab="Samples")
axis(1, at=c(5.5, 6, 6.5, 7), cex.axis=1.4)