如何用SAS生成双峰分布的数据呢？

shiyiming

求助各位高手，如何生成双峰分布的数据呢？
看到文献上是这样写的 Bimodal:  N(u=0.5;sigma=0.4) N(u=-0.5;sigma=0.4)

如何实现这样分布的数据呢？

谢谢！

shiyiming

假设mcmc在你的sas版本可用；假设双峰的分布为0.5 和0.5：

[code:2138sk06]data have; run;
proc mcmc data=have outpost=MixtureNorm seed=11 nmc=10000 nbi =2000;
      parm x 0;
      logp = logpdf('normalmix', x, 2, 0.5, 0.5, 0.5, -0.5, 0.4,  0.4);
      prior x  ~ general(logp);
      model general(0);
   run;[/code:2138sk06]
假设一个空集的输入数据。此处设定似然值为常数(flat)。故此posterior等同于prior， a mixture of normal distribution.

shiyiming

to jingju
非常谢谢你的回复
不过程序中nmc是最后生成的个数吧，那 nbi 是什么意思呢？
需要根据nmc相应的调整吗？

shiyiming

that is burning number. people thint it may not be necessary when you simply generate a distribution. but it is no harm to get rid of those first ones in my understanding.

shiyiming

京剧兄，码品好，人品也好，佩服佩服