geoR包中krige.bayes函数问题求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，向给位请教一个问题。本人在进行考虑空间关系的风险因子探索分析研究。
很多风险因子在空间上具有相关性，在进行回归分析的时候需要考虑这种相关性，故本人采用spatial linear gussian model. 该模型在可以通过geoR包中krige.bayes函数来实现，代码如下：

s<-read.geodata("D:/paper of mortality rate/EHP/R code/bayes.txt", header=TRUE,covar.col=c(5:12))
mod<-model.control(trend.d=trend.spatial(~ worker+density+income+slope+DEM+fault+MMI+house,s))
prio<-prior.control(phi.prior="r",phi.d=seq(0,100,l=51),tausq.rel.prior="u", tausq.rel.d=seq(0,10,l=51))
e<-krige.bayes(s,model=mod,prior=prio)

运行之后，出现了下列的warning:

警告信息：
In sample((1:n.points), n.posterior, replace = TRUE, prob = as.vector(phidist$probphitausq)) :
  Walker's alias method used: results are different from R < 2.2.0

不知何故？是不是phi和tausq的先验分布指定有误？请做过类似分析的朋友赐教！在此谢过了！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Bayes krige baye Rig bay 求助 函数 Bayes geoR krige

如果你正常出结果没有错误，只是这个warning的话，无视吧
它只是告诉你，版本不同，结果不同
仅此而已