winBUGS操作小步骤
(2013-05-01 02:10:58)WinBUGS-GUI
我们先找一个例子来实际地运行一下WinBUGS,感受一下基本的操作流程,然后再考虑高级的操作。
第一步,打开WinBUGS,通过菜单File -> New新建一个空白的窗口
第二步,在第一步中新建的空白窗口中输入三部分内容:模型定义、数据定义、初始值定义(代码见附录)
第三步,点击菜单Model -> Specification,弹出一个Specification Tool面板。
第四步,在第二步中的提到的那个窗口中,将model这个关键字高亮起来,点击check model。你会看到WinBUGS的左下角状态栏上显示”model is syntactically correct.”
第五步,把定义的data前的关键字list也高亮起来,点Specification Tool面板上的load data
第六步,改Specification Tool面板上的马尔科夫链的数目,默认为1就好了
第七步,点击Specification Tool面板上的compile
第八步,把定义的初始值中的list关键字也高亮起来,再点击Specification Tool面板上的load inits
第九步,关了Specification Tool面板
第十步,点击菜单Inference -> Samples,弹出一个Sample Monitor Tool面板。
第十一步,在Sample Monitor Tool面板的node中填要估计的参数名,这里可以是tau, alpha0, alpha1, b, 把它们一个一个填在node中,逐一点set。
第十二步,关了Sample Monitor Tool面板
第十三步,点击菜单Model -> Update,弹出一个Update Tool面板。
第十四步,将Update Tool面板中的updates改大点,比如50000,点update按钮。
第十五步,运行完后,关了Update Tool面板
第十六步,点击菜单Inference -> Samples
第十七步,在弹出的Sample Monitor Tool面板上选一个node
第十八步,点history看所有迭代的时间序列图,点trace看最后一次迭代的时间序列图,点auto cor看correlogram时间序列图,点stat看参数估计的结果