各位大虾,本人初涉MCMCregress, 在R中敲入:
Model.posterior2 <- MCMCregress(log(s.p) ~ (kVp + AirGap + Phantom + FOV + mA)^2 , data = spFrame, b0 = 0, B0 = 0.1, c0 = 2, d0 = 0.11, marginal.likelihood = "Chib95")
然后想用这个模型去测试test dataset 看准确度如何, 请问在R中代码是怎么样的呢,多谢! google了很久也没看到例子啊。。。。