我目前正在使用MCMCglmm软件包进行随机拦截逻辑回归。我已经指定:
mcmc <- MCMCglmm(Event ~ 1+Wind+Humidity, random = ~Engineer, data=spares.data,
family="categorical",nitt=20000,burnin=3000, thin=24,prior=list(R=list(V=1, fix=1), G=list(G1=list(V=1, nu=0.002))),pr=T,pl=T)
我现在的问题是:
的规范
prior=list(R=list(V=1, fix=1), G=list(G1=list(V=1, nu=0.002)))
给我以前分布的以下(数学)规范:
beta_0,beta_1,beta_2~Normal(0,10 ^ 10)
tau ^ 2~InvGamma(0.001,0.001)
提前致谢: - )