有没有高手做过这方面呀？查了很多资料都没看到有用的程序。大概步骤都知道。可具体用r实现去用gibbs抽取参数我还是做不出来。或者有人用winbugs做过这方面也好呀？