如何在Julia中执行Yule-Walker分析？是否有类Python statsmodels的等效方法？

等效实现说明 #

Julia的StatsBase.jl包提供了和Python statsmodels里statsmodels.regression.linear_model.yule_walker功能完全等效的yule_walker函数，能完成AR模型系数、残差方差等核心计算。

使用步骤 #

1. 先安装StatsBase包（如果还没装）：

   using Pkg
   Pkg.add("StatsBase")
   

2. 导入包并调用函数：
   拿一个随机生成的时间序列示例，计算2阶AR模型的Yule-Walker估计：

   using StatsBase
   
   # 构造示例时间序列数据
   x = randn(100)
   # 求解p=2的Yule-Walker方程
   ar_coeffs, residual_var = yule_walker(x, 2)
   

   • ar_coeffs会返回AR模型的系数数组
   • residual_var返回模型残差的方差值

可选参数说明 #

和statsmodels的实现对齐，这个函数也支持自定义参数：

• method: 指定求解方程的算法，可选:cholesky（默认）、:levinson
• demean: 布尔值，控制是否先对输入数据去均值，默认true

示例用法：

# 使用Levinson算法求解，且不对数据去均值
ar_coeffs, residual_var = yule_walker(x, 2; method=:levinson, demean=false)

内容的提问来源于stack exchange，提问作者Jose Manuel de Frutos