Draw random number from Normal distribution past a threshold value

Question

Imner Renmi 2016-6-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/291745-draw-random-number-from-normal-distribution-past-a-threshold-value

评论： Roger Stafford 2016-6-23

Hi,

I want to draw values that are larger than 1.96 from a normal distribution. Anyone know how I can do this in matlab. Is there a built-in function to do this?

Note, of course I could use rejection sampling, but I want to draw a vector of 10000 values larger than 1.96 and I want to do this 5000 times. Using rejection sampling could be extremely time-consuming and I want to avoid that.

Thank you very much for your help.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Roger Stafford 2016-6-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/291745-draw-random-number-from-normal-distribution-past-a-threshold-value#answer_226403

编辑：Roger Stafford 2016-6-22

在 MATLAB Online 中打开

If your normal distribution without the truncation would have been mean = mu and standard deviation = sig, then do this:

   r0 = normcdf(1.96,mu,sig); % The probability that sample is below 1.96
   x = norminv((1-r0)*rand(1,n)+r0,mu,sig);

The value n is the number of desired samples. There is no rejection in this method. You should realize that this is no longer a normal distribution and its mean and standard deviation are different from 'mu' and 'sig'.

2 个评论
显示无隐藏无

Imner Renmi 2016-6-23

Thanks. One question, why do you have the "+r0" in the second line?

Roger Stafford 2016-6-23

This is to provide a statistically uniform distribution between r0 and 1 as input to the 'norminv' function so that it will generate only values starting at 1.96 and on higher in accordance with that end of the normal distribution. Try it and see.

请先登录，再进行评论。