在R中获取prais.winsten时间序列回归的观测数(N)

因为prais.winsten包返回的模型对象没有实现基础R的nobs()方法，所以直接调用会失效，你可以用以下两种简单方法获取模型实际使用的观测数：

• 提取模型残差的长度：

  length(pw_model$residuals)
  

  模型的残差向量只包含实际参与回归计算的有效观测，它的长度就是最终用的观测数。

• 提取模型拟合值的长度：

  length(pw_model$fitted.values)
  

  拟合值和残差是一一对应的，长度完全一致，结果和上面的方法相同。

举个完整的示例：

library(prais)
# 构造面板数据示例
panel_data <- data.frame(
  group = rep(1:8, each = 6),
  time = rep(2015:2020, 8),
  dep_var = rnorm(48),
  ind_var1 = rnorm(48),
  ind_var2 = rnorm(48)
)
# 拟合Prais-Winsten模型
pw_fit <- prais_winsten(dep_var ~ ind_var1 + ind_var2, data = panel_data, index = c("group", "time"))

# 获取观测数
length(pw_fit$residuals)

内容的提问来源于stack exchange，提问作者Dale