求《model-based geostatistics》书本数据

喝汤汤 发表于 2015-4-23 09:54
没有啦，贝塔的估计值完全不一样，最大似然估计的结果是有点不一样。就是钙含量分布的那个例子。如果你发 ...

R给出的beta值是这样的

> m2$beta
intercept    covar1    covar2
38.415685  8.542968 15.580809

因为area是个factor，有三个level，（1，2，3），R把第一个level作为参照值，所以area等于1的估计值就是intercept，area等于2的估计值就是intercept加上covar1，area等于3的估计值就是intercept加上covar2

结果有些小的差别，可能作者调整了一下数据，不过大体上是基本吻合的

qoiqpwqr 发表于 2015-4-23 13:06
R给出的beta值是这样的

> m2$beta

恩恩，您的结果和我运行的是一样的。十分感谢，您是第一个帮我解答的人，也是唯一一个