全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
2643 4
2011-10-14
悬赏 15 个论坛币 已解决
哈,不知道自己标题的表述对不对……
请教各位先辈啦~

比如说呢,我之前估计出来了一个正态分布的均值和方差,存在了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.

谢谢!!!~



最佳答案

sungmoo 查看完整内容

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

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2011-10-14 12:21:32
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-10-14 12:24:23
自己顶一个~~期待高手help一下~~o(∩_∩)o
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-10-14 12:44:07
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
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-10-14 16:56:08
  *-----------------           
  *-截断型半正态分布(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
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群