Spbayes 包-spLM 函数-spRecover( 中的错误，c++ 错误：dpotrf 失败

Question

我在 spBayes 包中的 spRecover 方面遇到问题。 spLM 运行得很好。 SpRecover 针对指数和球形模型运行，但是在高斯情况下出现“c++ 错误：dpotrf 失败”错误。我没有重复的坐标。有谁知道如何继续吗？

我尝试调整坐标舍入以消除所有重复项并添加抖动，但我仍然遇到同样的问题。

sp <- spLM(daily_data ~ data.type, data = data_tot,coords=coords,starting=list("phi"= 0.7249594 ,"sigma.sq"=30.88157, "tau.sq"= 7.122322), tuning=list("phi"=0.1,"sigma.sq"=0.05, "tau.sq"=0.05), priors=list("phi.Unif"=c(((0.7249594)/7.5),(9)), "sigma.sq.IG"=c((23.75925),(30.88157)), "tau.sq.IG"=c((23.75925),7.122322)), cov.model="gaussian",n.samples=n.samples) round(summary(mcmc(sp$p.theta.samples))$quantiles,3) 

#composition sampling to get posterior samples of marginalized regression coefficients and the spatial effects 
burn.in <- floor(0.75*n.samples) 
fin.sp <- spRecover(sp, start=burn.in)

Answer 1

你有没有弄清楚这一点？我遇到了同样的错误！