一个if循环的小问题

想编一个小程序，在0到2的范围内，随机生成两个标准正态分布的数，但不成功。

1. xstar=rnorm(2)
   
2. for(i in 1:2)
   
3. if((xstar[i]<=2)&(xstar[i]>0)){x[i]=xstar[i]}
   
4. x

复制代码

你这个是reject sampling， 要考虑sample不到的情况，用while循环

rnormtrunc <- function(low, high) {
        x <- rnorm(1)
        while (x < low || x > high) x <- rnorm(1)
        x
}

x <- numeric(2)
for (i in seq_along(x)) x[i] <- rnormtrunc(0, 2)

“在0到2的范围内，随机生成两个标准正态分布的数”这句话本身就是有问题的，既然是服从正态分布，又何来有界。

dasad

1. a=rnorm(20)
   
2. a=a[a<=2&a>0]
   
3. x=sample(a,2)
   
4. x

复制代码

求证1加1 发表于 2014-2-9 14:49
“在0到2的范围内，随机生成两个标准正态分布的数”这句话本身就是有问题的，既然是服从正态分布，又何来有 ...

  截尾分布，今天上午碰到的。

耕耘使者 发表于 2014-2-9 17:26
截尾分布，今天上午碰到的。

http://www.stat.ucla.edu/~dinov/ ... jectionSampling.pdf

建议用while循环：

1. normtrunc <- function(n, a, b, mu, sig){
   
2.   x <- rnorm(n, mu, sig)
   
3.   xvalue <- (x < a) | (x > b)
   
4.   while(sum(xvalue) > 0){
   
5.     x[xvalue] <- rnorm(sum(xvalue), mu, sig)
   
6.     xvalue <- (x < a) | (x > b)
   
7.   }
   
8.   return(x)
   
9. }
   
10. 
    
11. n <- 2
    
12. a <- 0
    
13. b <- 2
    
14. mu <- 0
    
15. sig <- 1
    
16. normtrunc(n, a, b, mu, sig)

复制代码

求证1加1 发表于 2014-2-9 17:53
http://www.stat.ucla.edu/~dinov/courses_students.dir/04/Spring/Stat233.dir/STAT233_notes.dir/Rejec ...

学习了

ntsean 发表于 2014-2-9 15:42
你这个是reject sampling， 要考虑sample不到的情况，用while循环

rnormtrunc

ntsean兄，能否详细解释一下每一步？
请原谅我是初学者，谢谢。