# R SCRIPT FOR WEEK 3 DEMONSTRATION. source("wk3funcs.r") par(mfrow=c(2,2)) set.seed(1); plot_thetas(sample_prior(2000)) set.seed(1); plot_thetas(sample_prior(2000,ordered=TRUE)) set.seed(1); plot_thetas(sample_posterior(2000,10,4,7)) set.seed(1); plot_thetas(sample_posterior(2000,10,4,7,ordered=TRUE))