利用附文中的数据,一元线性回归模型的方差的最大似然估计为12.2035179,而用winbugs做出来的均值或中位数结果都与这个相差较大(一般15左右或15以上,R中的MCMCpack包也是一样的结果),而两种方法估计出来的截距和斜率系数却相差不大,不知道为什么?怎么解释?
难道方差的MLE和贝叶斯估计就会有如此差异?
附:winbugs代码
model
{
alpha ~ dnorm( 0.0,1.0E-6)
beta ~ dnorm( 0.0,1.0E-6)
for( i in 1 : 15 ) {
mu <-alpha+beta*x
y ~ dnorm(mu,tau)
}
tau ~ dgamma(0.001,0.001)
sig2<-1/tau
}
list(x=c(135.1,139.9,163.6,146.5,156.2,156.4,167.8,149.7,145,148.5,165.5,135,153.3,152,160.5),
y=c(32,30.4,46.2,33.5,37.10,35.5,41.5,31,33,37.2,49.50,27.6,41,32,47.2))
list(alpha=0,beta=0,tau=1)