我使用mcmc package中的metrop函数来估计线性回归系数的后验分布（13个解释变量）。
系数的先前分布都是gammas（选择形状和速率以使分布看起来一致）。初始状态下的后验对数似然是有限的。初始状态接近于伽马分布的均值。
但是，启动mcmc函数时，至少9次中有10次出现以下错误消息：
system.time出错（out＆lt; - .Call（＆＃34; metrop＆＃34;，func1，initial，nbatch，blen，：
logh：func返回NA或NaN
另外：
警告讯息： 在dnorm（y，mean = pred，sd = ecart_type，log = T）：产生NaNs
当算法工作时，接受率非常弱（~0.001）。
有谁知道可能是什么问题？这可能是因为提案功能的支持包括负值，这与伽马先验不一致吗？