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)で作成されているようです。時系列をこのコードに通したときに何が起こっているのかを完全に理解する方法はありますか?
ありがとう、Jochem