时间，空间自相关会影响数据的独立性，也是审稿人问无可问时，经常顺道来一嘴的问题，当然更是时常困扰我们，一句话把我们打进万丈深渊的的问题之一。上周五（2月29号），Science上发表的Divergent responses of pelagic and benthic fish body-size structure to remoteness and protection from humans这篇文章，研究了两类群鱼类对不同因素相应的变化，其中就涉及到了时空自相关数据的处理问题(Letessier et al. 2024)。

该研究考虑了众多因素，以及一些因素之间的交互作用，对鱼类的体型大小的影响。其主模型可以用以下公式表示：

在这些模型中，作者采用nlme包下的gls函数，对数据存在的时空自相关进行了考虑，其文中讲到：

由于gls中提供了不同的自相关类型，作者便将所有的类型都跑了一遍，然后用AIC值选择哪种结构与其数据最为匹配，即按照以下代码操作：

### 定义模型因变量和自变量的关系

formula<-"betaslope ~ bruvs*protection_use*(logBathy + Slope + logDistP +

logDistSM + logTTM + logDistCR + logDistC+ poly(logSST, 2) + poly(logCHL,2))"

### 构建具有不同自相关结构的模型

mod1 <- nlme::gls(formular, data = data, correlation = nlme::corAR1(form=~1))

mod2 <- nlme::gls(formular, data = data, correlation = nlme::corExp(form=~1))

mod3 <- nlme::gls(formular, data = data, correlation = nlme::corGaus(form=~1))

mod4 <- nlme::gls(formular, data = data, correlation = nlme::corLin(form=~1))

mod5 <- nlme::gls(formular, data = data, correlation = nlme::corRatio(form=~1))

mod6 <- nlme::gls(formular, data = data, correlation = nlme::corSpher(form=~1))

mod7 <- nlme::gls(formular, data = data, correlation = nlme::corAR1(form=~mean_long+mean_lat))

mod8 <- nlme::gls(formular, data = data, correlation = nlme::corExp(form=~mean_long+mean_lat))

mod9 <- nlme::gls(formular, data = data, correlation = nlme::corGaus(form=~mean_long+mean_lat))

mod10 <- nlme::gls(formular, data = data, correlation = nlme::corLin(form=~mean_long+mean_lat))

mod11 <- nlme::gls(formular, data = data, correlation = nlme::corRatio(form=~mean_long+mean_lat))

mod12 <- nlme::gls(formular, data = data, correlation = nlme::corSpher(form=~mean_long+mean_lat))

###根据AIC值，选出与数据最为匹配的结构

a <- AIC(mod1, mod2, mod3, mod4, mod5, mod6, mod7, mod8, mod9, mod10, mod11, mod12)

print(a)

欢迎大家关注我的个人微信公众号：二傻统计

参考文献：

Letessier, T. B., D. Mouillot, L. Mannocci, H. Jabour Christ, E. M. Elamin, S. M. Elamin, A. M. Friedlander, A. Hearn, J.-B. Juhel, A. R. Kleiven, E. Moland, N. Mouquet, P. J. Nillos-Kleiven, E. Sala, C. D. H. Thompson, L. Velez, L. Vigliola, and J. J. Meeuwig. 2024. Divergent responses of pelagic and benthic fish body-size structure to remoteness and protection from humans. Science 383:976-982.