大家好，向给位请教一个问题。本人在进行考虑空间关系的风险因子探索分析研究。
很多风险因子在空间上具有相关性，在进行回归分析的时候需要考虑这种相关性，故本人采用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的先验分布指定有误？请做过类似分析的朋友赐教！在此谢过了！