R/ST_Ordered.R

In ACSWR: A Companion Package for the Book "A Course in Statistics with R"

ST_Ordered <-
function(N,x,p_x){
  x <- x[order(p_x,decreasing=TRUE)]
  F_x <- cumsum(sort(p_x,decreasing=TRUE))
  disc_sim <- numeric(length=N)
  for(i in 1:N){
    temp <- runif(1)
    disc_sim[i] <- x[min(which(F_x>temp))]
  }
  return(disc_sim)
}