zhang111111
如何用s-plus产生截尾正态分布?
请多指教!
yihui
这句话的宾语很奇怪啊,什么叫产生一个分布?产生服从某种分布的随机数?计算分布函数值?计算分位数?计算密度函数值?
zhang111111
如何用S-PLUS产生服从正态截尾分布的随机数?比如如何产生截尾于(-1,1)区间上的正态分布随机数
rtist
呵呵,一楼的问题提得很在理。
二楼的问题有很简单的解决方案,虽然computationally不够efficient,但是足够满足常见应用了,做rejection sampling啊。
yihui
如果是一维的连续型分布,那么随机数还是很好产生的,主要就是做一个分位变换(Quantile Transformation),这种随机数产生方法适用于任何一维连续分布。过程很简单,先产生均匀分布U(0, 1)的随机数x,然后假设invG(.)是分布G(.)的分位数函数(或者称之为广义逆Generalized inverse),定义为invG(y)=inf{x: G(x)>=y},那么invG(x)就服从G分布了。简单说来就是把广义逆函数作用在均匀分布上。
zhang111111
谢谢各位回答,我想问有没有简单的s-plus语句可以产生截尾于(-1,1)区间上的正态分布随机数,谢谢!
rtist
simply throw away any rnorm not in the range