R语言acf()函数i、j参数含义与正确使用方法

i和j不是acf()计算函数本身的入参，而是acf类对象的子集提取方法（[.acf）的专属参数。
R 3.6.2版本的stats包中，这两个参数的说明被整合在acf函数的同一篇文档里，没有单独标注参数归属，这也是文档没有提供直接传参示例的原因——这种用法本身就不符合参数的设计逻辑。
你直接把i传给acf()时，因为acf()本身没有定义这两个形参，参数会通过函数的...通道被透传给底层的绘图逻辑plot.acf()，而绘图函数的合法参数列表里完全没有i、j，因此会抛出“不是合法图形参数”的警告，绘图和计算逻辑会直接忽略传入的值，最终结果不受任何影响。

结合官方文档的定义，两个参数的设计用途完全匹配描述，但只能在acf对象的子集提取环节使用：

• i：指定需要保留的滞后阶对应的存储位置索引，也就是acf结果矩阵里要保留的行。注意acf结果默认从滞后0阶（自相关系数恒为1）开始存储，第1行对应lag0、第2行对应lag1、第3行对应lag2，以此类推，滞后值和位置索引永远差1。
• j：指定需要保留的序列，可通过序列名称或列编号指定。仅在输入多变量时间序列、同时计算多序列自/互相关的场景下生效，单序列输入时该参数无实际作用。

错误写法（完全不生效）：

time_series = rnorm(100)
# 直接传i给acf()，参数会被透传给绘图逻辑，抛出警告且不生效
acf(time_series, i=c(1,2,4,7))

正确执行流程：

1. 先关闭自动绘图，计算完整的acf结果，拿到acf类对象
2. 用方括号子集方法传入i、j参数筛选结果
3. 对筛选后的对象绘图

# 第一步：计算完整acf，关闭自动绘图
acf_full <- acf(time_series, plot = FALSE)
# 第二步：筛选要保留的滞后阶：目标是lag1/2/4/7，对应位置索引为2/3/5/8（第1位是lag0）
acf_filtered <- acf_full[c(2,3,5,8), ]
# 第三步：绘制筛选后的结果，此时输出仅包含指定的滞后阶
plot(acf_filtered)

• 不要直接把滞后阶数值当成位置索引传给i。手动数索引很容易因为差1的固定偏移出错，更稳妥的写法是先匹配目标滞后值拿到对应索引：

target_lag <- c(1,2,4,7)
keep_pos <- which(acf_full$lag %in% target_lag)
acf_filtered <- acf_full[keep_pos, ]

• 多变量场景下j的用法示例：针对3个序列组成的时间序列对象，用acf_full[, "序列B"]即可筛选出仅和序列B相关的所有自/互相关结果，排除其他序列的交叉计算结果。
• 这套i、j的子集能力是R 3.6.2版本新增的，更早的R版本中acf类对象不支持子集提取，也不存在这两个参数。

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