了解mcmc实现过程的一些部分

Question

我正在研究

MCMC in SAS

的一种实现 https://support.sas.com/documentation/cdl/en/statug/67523/HTML/default/viewer.htm#statug_mcmc_examples13.htm

它实现了

proc mcmc data=seeds outpost=postout seed=332786 nmc=20000;
   ods select PostSumInt;
   parms beta0 0 beta1 0 beta2 0 beta3 0 s2 1;
   prior s2 ~ igamma(0.01, s=0.01);
   prior beta: ~ general(0);
   w = beta0 + beta1*seed + beta2*extract + beta3*seed*extract;
   random delta ~ normal(w, var=s2) subject=ind;
   pi = logistic(delta);
   model r ~ binomial(n = n, p = pi);
run;

我不明白代码线

ods select PostSumInt;

和声明

parms beta0 0 beta1 0 beta2 0 beta3 0 s2 1;

即

beta0 0

的目的是什么？它们是起始值吗？

此外，

beta

系数似乎没有先验定义，但我看到有一个通用定义

prior beta: ~ general(0);

这里的

beta

如何与

beta0, beta1

等相关联，什么是

general(0)

？

我真的很感谢您对上述代码的任何澄清。

Answer 1

我只告诉你一件事，“ods”，它代表输出交付系统，用于在文件中以报告样式获取输出，例如pdf、rtf、Excel和文本文件。好吗？