我想从一个用MCMCglmm估计的模型中模拟感兴趣的数量,大致就像Zelig包那样。在Zelig中,您可以为独立值设置所需的值,软件可以为结果变量(期望值、概率等)计算结果。z = 10simulation <- sim(m1.zelig, x = s1)正如我们所看到的,# Same model with MCMCglmm
library(MCMC
我有一个具有以下(近似)形式的多元模型: cbind(OFT1, MIS1, PC1, PC2) ~ thin=100, #interval at which the Markov chain is stored为了出版目的,我被要求报告Gelman的统计数据我一直想跑:但是,我知道这个错误:
Error i
我使用嵌套循环在R中执行多次MCMCglmm运行,将后验分布输出到数据帧中。外循环执行MCMC过程,生成后验分布,然后内循环将每一个后验分布提升到分离数据帧(其中F是链每次运行的固定效果中后验分布的数目,n是要执行的随机数);因此,将有F(在这种情况下,4)数据帧,每个数据帧都有R = list(V = 1, nu = 0.002))
# Cha