按照特定要求生成随机分布

xingxf 发表于 2017-3-13 09:26
非常感谢你的思路。

至于你说在一个文件里一次生成1000次的问题能具体一点么？10000个观察值，1000个新 ...

set obs 100
gen g=mod(_n,10)
sort g
by g: gen t=_n
gen x=rnormal()
reshape wide x, i(t) j(g)

非常感谢，但是reshape实在不是提高效率的方法，这恐怕是Stata处理起来比较慢的函数之一。10个观察值，11个变量的long-wide转换用不了多长时间，但是10000个观察值，1000个变量的long-wide转换，reshape就得几十分钟吧。这还仅仅是reshape，再加上前面的simulation的代码。这种方法的效率远不如merge 1000次啊。

xingxf 发表于 2017-3-14 07:10
非常感谢，但是reshape实在不是提高效率的方法，这恐怕是Stata处理起来比较慢的函数之一。10个观察值，11个 ...

那你就merge咯。实际做了才知道哪个快。实在太慢的话改用eviews或者matlab吧。

夏目贵志 发表于 2017-3-14 07:18
那你就merge咯。实际做了才知道哪个快。实在太慢的话改用eviews或者matlab吧。

我觉得哪个统计软件并不是决定速度的决定性因素，核心还是怎么写代码。我以前用Stata做一个条件比较复杂的计数，循环程序跑了45小时，后来改用基本函数但是比较占内存，但是耗时只用了10分钟。所以关键是算法。

@夏目贵志你好 请问如何根据均值 标准误随机生成一个变量