广义最小二乘结果解释答案

【问题标题】：Generalized least squares results interpretation广义最小二乘结果解释
【发布时间】：2018-07-06 14:25:58
【问题描述】：

我检查了我的线性回归模型（WMAN = 物种，WDNE = 海面温度）并发现了自相关，因此我尝试使用以下脚本进行广义最小二乘；

library(nlme)
modelwa <- gls(WMAN ~WDNE, data=dat, 
               correlation = corAR1(form=~MONTH),
               na.action=na.omit)
summary(modelwa)

我比较了两个模型；

> library(MuMIn)
> model.sel(modelw,modelwa)
Model selection table 
        (Intrc)   WDNE class na.action  correlation df   logLik   AICc delta
modelwa   31.50 0.1874   gls   na.omit crAR1(MONTH)  4 -610.461 1229.2  0.00
modelw    11.31 0.7974    lm   na.excl               3 -658.741 1323.7 94.44
        weight
modelwa      1
modelw       0
Abbreviations:
na.action: na.excl = ‘na.exclude’
correlation: crAR1(MONTH) = ‘corAR1(~MONTH)’
Models ranked by AICc(x)

我相信结果表明我应该使用 gls，因为 AIC 较低。

我的问题是，我一直在报告 F-value/R²/p-value，但是 gls 的输出没有这些？

如果有人可以帮助我解释这些结果，我将不胜感激？

> summary(modelwa)
Generalized least squares fit by REML
  Model: WMAN ~ WDNE 
  Data: mp2017.dat 
       AIC      BIC    logLik
  1228.923 1240.661 -610.4614

Correlation Structure: ARMA(1,0)
 Formula: ~MONTH 
 Parameter estimate(s):
     Phi1 
0.4809973 

Coefficients:
                Value Std.Error  t-value p-value
(Intercept) 31.496911  8.052339 3.911524  0.0001
WDNE         0.187419  0.091495 2.048401  0.0424

 Correlation: 
     (Intr)
WDNE -0.339

Standardized residuals:
      Min        Q1       Med        Q3       Max 
-2.023362 -1.606329 -1.210127  1.427247  3.567186 

Residual standard error: 18.85341 
Degrees of freedom: 141 total; 139 residual
>

【问题讨论】：

我认为你跑题了，如果我没记错的话，这是关于统计分析的，应该经过交叉验证

标签： r linear-regression

【解决方案1】：

我现在已经克服了自相关的问题，所以我可以使用lm()

将残差的 lag1 作为 X 变量添加到原始模型中。这可以使用DataCombine 包中的slide 函数来完成。

library(DataCombine)
econ_data <- data.frame(economics, resid_mod1=lmMod$residuals)
econ_data_1 <- slide(econ_data, Var="resid_mod1", 
NewVar = "lag1", slideBy = -1)
econ_data_2 <- na.omit(econ_data_1)
lmMod2 <- lm(pce ~ pop + lag1, data=econ_data_2)

这个脚本可以在here找到

【讨论】：