Сгенерировать вторую половину нормального распределения, используя метод приемки-отклонения

Question

Следующий код R дает мне только половину нормального распределения; что я должен изменить на код, чтобы получить вторую половину?

halfnormal <- function(n){ 
    vector <- rep(0,n) 
    for(i in 1:n){ 
     uni_random <- runif(2) 
     y <- -log(uni_random) 
     while(y[2] < (y[1]-1)^2/2){ 
      uni_random <- runif(2) 
      y <- -log(uni_random) 
     } 
     vector[i] <- y[1] 
    } 
    vector 
} 

output <- halfnormal(1000) 
hist(output) 

Answer 1

Если вы настаиваете на использование этого кода для генерации стандартного нормального (не рекомендуются, так как rnorm будет быть намного быстрее и точнее), просто точечный продукт, что весь вектор вектор равной длины, состоящий из случайных значений (-1, +1).

Кстати, полунормальное явление также известно как Chi distribution (с 1 степенью свободы).

Answer 2

Это похоже на алгоритм Зиггурата с модификацией Марсалья, но это немного отличается? Если вы не хотите использовать любые гарантируемые к работе генераторов случайных чисел в R, возможно, это работает:

halfnormal <- function(n){ 
     vector <- rep(0,n) 
     for(i in 1:n){ 
      uni_random <- runif(2) 
      y <- -log(uni_random) 
      while(y[2] < (y[1]-1)^2/2){ 
       uni_random <- runif(2) 
       y <- -log(uni_random) 
      } 
      vector[i] <- sample(c(-1,1),size=1)*y[1] #randomly select the tail 
     } 
     vector 
    } 

    output <- halfnormal(1000) 
    hist(output)