TWSNMP FC:R言語の便利な関数STLを組み込むための調査

今朝も３時に猫が起こしにきてくれました。今朝はお世話した後、二度寝できました。４時にまた起こしにきてくれました。
まずは昨日の夜、寝る前に気になっていた外部コマンド実行のポーリングを修正しました。

外部コマンド実行ポーリングの改善 · twsnmp/twsnmpfc@77adeeb

です。
その後、前々から組み込みたかったR言語のSTL関数

stl function - RDocumentation

のような分析機能の調査を始めました。この関数はかなり便利で時系列データを

(Python編) 時系列データをサクッとSTLでトレンド・季節性に分解

のように分解して分析できるものです。R言語以外だとPythonの実装が見つかりました。JavaScriptかGO言語でできないか調べているとありました。

project-anomalia/stl

まずはサンプルコード試してみることにしました。元のサンプルコードだと結果が数値なのでちゃんと動いているかよくわりません。そこで、サンプルコードを

package main

import (
	"encoding/csv"
	"fmt"
	"os"
	"strconv"

	"github.com/chewxy/stl"
)

func main() {
	var data []float64
	f, _ := os.Open("co2.csv")
	r := csv.NewReader(f)

	r.Read() // read header
	for rec, err := r.Read(); err == nil; rec, err = r.Read() {
		// here we're ignoring errors because we know the file to be correct
		if co2, err := strconv.ParseFloat(rec[0], 64); err == nil {
			data = append(data, co2)
		}
	}

	// The main function:
	res := stl.Decompose(data, 12, 35, stl.Additive(), stl.WithRobustIter(2), stl.WithIter(2))
	printData(res.Data)
	printData(res.Seasonal)
	printData(res.Trend)
	printData(res.Resid)
}

func printData(data []float64) {
	fmt.Println()
	fmt.Printf("[")
	for i, d := range data {
		if i%20 == 0 {
			fmt.Println()
		}
		fmt.Printf("[%d,%f],", i, d)
	}
	fmt.Println("]")
}

のように書き換えてechartsで表示しやすいデータを出力するようにしました。このデータをechartsのサンプルページ

Examples - Apache ECharts

にコピペしてグラフに表示してみました。

それなりに分解できています。使えそうです。
ここで今朝は時間切れです。猫が後ろで「終了です。」と言っています。
明日は、このライブラリを使ってポーリング結果のSTL分析を組み込んでみようと思います。
明日に続く