Draw random number from Normal distribution past a threshold value
2 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Imner Renmi
el 22 de Jun. de 2016
Comentada: Roger Stafford
el 23 de Jun. de 2016
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 comentarios
Respuesta aceptada
Roger Stafford
el 22 de Jun. de 2016
Editada: Roger Stafford
el 22 de Jun. de 2016
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 comentarios
Roger Stafford
el 23 de Jun. de 2016
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.
Más respuestas (1)
Walter Roberson
el 22 de Jun. de 2016
The method of operation is not documented; it might use a rejection method internally.
0 comentarios
Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!