昔 SAS/IML で離散型 HMM にチャレンジしたときに できたもの (Dropbox)、対数取る処理 (アンダーフロー防止のための Scaling) とかする前に遅くてやめてしまった。 シンボル行列の次元が 10000×c ぐらいになると激遅でもうだめですね。 IML は計算速度が・…