# NOT RUN {
n <- 500
theta <- rnorm(n)
ses <- sqrt(rgamma(n,shape=1,scale=1))
XX <- theta + ses*rnorm(n)
dd <- cbind(XX,ses)
rvs <- rvalues(dd, family = gaussian)
TopList(rvs, topnum = 12)
TopList(rvs, topnum = 15, sorted.by = "MLE")
# }
Run the code above in your browser using DataLab