r - R中的STL趋势分量是如何计算的？

Question

我试图了解 STL 时间序列函数的作用，特别是趋势是如何定义的。当我阅读 STL 代码中发生的事情时，我会看到以下代码行：

z <- .Fortran(C_stl, as.double(x), as.integer(n), as.integer(period), 
as.integer(s.window), as.integer(t.window), as.integer(l.window), 
s.degree, t.degree, l.degree, nsjump = as.integer(s.jump), 
ntjump = as.integer(t.jump), nljump = as.integer(l.jump), 
ni = as.integer(inner), no = as.integer(outer), weights = double(n), 
seasonal = double(n), trend = double(n), double((n + 
2 * period) * 5))

这里有一些关于趋势的事情，但看起来实际输出是在一段外部代码 (C_stl) 中创建的。当我通过这段代码传递我的时间序列时，我知道如何全面了解正在发生的事情吗？

谢谢，约赫姆

Answer 1

是的，它是 Fortran 代码。

代码可以在源代码中找到：例如这里。