使用 nlme 包中的 gls 函数校正空间自相关的错误

Question

我正在尝试使用以下代码使用 nlme 包中的 gls 函数连接空间自相关：

modelExp<- gls(avg_diff ~ burnsev/Species, data= burn, na.action = na.omit,
                              correlation = corExp(form = ~Lat + Long, nugget = TRUE)

但是，我不断收到错误消息：

Error in getCovariate.corSpatial(object, data = data) : 
  cannot have zero distances in "corSpatial"

尽管已从我的数据框中删除了所有重复的坐标，但此错误仍然存​​在。有谁知道如何解决这个错误？可能存在空间分辨率的问题吗？我的一些点由

谢谢！

Answer 1

发生这种情况是因为您的 Lat Lon 在您的数据集中重复。通常会发生这种情况，因为您多年来一直从点位置采集样本。假设您的分组因子是年份，一个解决方案是执行以下操作：

modelExp<- gls(avg_diff ~ burnsev/Species, data= burn, na.action = na.omit,
                              correlation = corExp(form = ~Lat + Long|year, nugget = TRUE)