#Not run due to time constraints
# \donttest{
Rs <- ruars(20, rvmises, kappa = 10)
#Compare the region size of the moment based theory mean estimator to the
#Bayes region.
region(Rs, method = "direct", type = "theory", estimator = "mean", alp=0.1, m = 100)
bayesCR <- region(Rs, type = "Mises", method = "Bayes", estimator = "mean", S0 = mean(Rs),
kappa0 = 10, tuneS = 5000, tuneK = 1, burn_in = 1000, alp = .01, m = 5000)
bayesCR$Radius #Region size is give by "Radius"
bayesCR$Shat #The Bayes region is centered around the posterior mode: "Shat"# }
Run the code above in your browser using DataLab