这几天要做这个hierachical model的MCMC,是我的毕业论文，都要被搞疯了
模型是这样的，有i=1:S个站点的从t=1年到t=50年的数据 {X(i,t)}
模型的第一层是 X(i,t)~GEV(Mu(i,t),Sigma(i,t),Xi(i)); i=1:S; t=1:50;
第二层是对Mu(i,t)和Sigma(i,t)的建模，他们既时间相关，又空间相关；
不妨令Mu(i,t)=a(i)+gg(i)+ mm(t),其中a(i)是固定效应，gg(i)用CAR模型，是个spatial random
effect, mm(t)是个时间效应。
用Bayesian的方法写出updating distribution的话，就可以用MCMC来做，因为有多重参数，所以在外面用到
Gibbs sampling,而在里面每次更新用到Metropolis-hasting算法(或其它)。
因为参数很多，更新的东西很多，最重要的是，GEV，Inverse Gammma...这些分布，在Winbugs里面都没有，所
以没法用WinBugs来做，只能用R写，但是我初次接触到MCMC就面对这么个复杂的问题，没有个示例，实在是无从
下手。不知道怎么控制更新，不知道怎么看何时收敛或者posterior之类的。
盼有用过MCMC的大神能解救小可，或者能告诉我R中用什么函数能够指定某个分布，让它以某种算法更新。
(我已经研究过Metropolis-hasting等算法，但我想，R上大能之士如此多，不会还要我来手动按算法写个
updating的函数吧，不然就哭死了)
我已经粗略研究过R中的MCMCpack, coda ,evdbayes等包，但是还是一团浆糊。示例都太过简单，或者特殊。
告诉我一些函数，或者一些包，或者一个范例都可，跪谢啊！！！
好心人啊，好人有好报啊！！！
