请教stata中怎么提取截断了的正态随机数

哈，不知道自己标题的表述对不对……
请教各位先辈啦~

比如说呢，我之前估计出来了一个正态分布的均值和方差，存在了r(m),r(s)里
那么，我现在要提取随机数，就是用rnormal(r(m),r(s))命令嘛，我知道这个

但是，现在考虑的是一个censor的情况
就是我想要提前出来的随机数，
不是从这个正态分布负无穷到正无穷区间提取的
而是想要从某个数a到正无穷提取    即 [a,+o0)
那么是否有这个命令呢？
还是根据 rnormal修改？
帮助文件如下：
    rnormal()

       Range:        c(mindouble) to c(maxdouble)    <----这个，我希望是 ·a‘ toc(maxdouble)
       Description:  returns standard normal (Gaussian) random variates,
                     i.e., variates from a normal distribution with a mean of
                     0 and a standard deviation of 1.

谢谢！！！~

并帖
https://bbs.pinggu.org/thread-1107862-1-1.html

自己顶一个~~期待高手help一下~~o(∩_∩)o

sungmoo 发表于 2011-10-14 12:35
并帖
https://bbs.pinggu.org/thread-1107862-1-1.html

哇~~~~~我都没搜索到~~！
~>.<~
好吧，谢谢版主~~
我把结果粘过来好了……

其中有例子（左截点-2，右截点2，方差1/4，期望0）：

set obs 10000
gen c = -2          /*left truncation point*/
gen d = 2           /*right truncation point*/
gen stdv = 0.5   /*standard deviation*/
gen v = uniform()
gen c_bar = c/stdv
gen d_bar = d/stdv
gen v_bar = (normal(c_bar) + (normal(d_bar) - normal(c_bar))*v)
gen u = stdv*invnorm(v_bar)


下载
http://fmwww.bc.edu/repec/bocode/g/gentrun.ado
http://fmwww.bc.edu/repec/bocode/g/gentrun.hlp

可生成左右截断的标准正态分布随机数。

*上述过程可以简化为一步：

g u=0.5*invnormal(normal(-2/0.5)+(normal(2/0.5)-normal(-2/0.5))*uniform())


以上解答来自于sungmoo

  *-----------------           
  *-截断型半正态分布（truncated standard normal distribution）
  *-----------------
   
        help gentrun
   
        *-这在模拟分析随机边界模型(Stochastic Frontier Model)时非常有用
    * 相关性质参见 Olive(2005)  Applied Robust Statistics, Chapter 4
      
       clear                  
       set obs 1000
       gentrun x, left(0)
       sum x
       histogram x
      
       gentrun y, right(-0.5)
       sum y
       histogram y