Example Regression
４次回帰プログラムの例
回帰計算をしたいデータファイルを表計算ソフトで作っていたとしましょう。サンプルですからパスを指定してプログラムからファイルを読み込んで、Row
Dataに表示します。"n th Poly fit"の入力では、１〜４次まで指定できます。"X-Y
Graph"には生データのポイントと回帰データ、回帰曲線が表示されます。"Coeff"の表示欄には０〜４次の係数が表示されます。
このプログラムの機能は（０）ファイルの読み込み、（１）回帰計算、（２）回帰曲線の計算、（３）グラフ表示の４個の機能ブロックに分けられ流れが直線的なので、シーケンス構造を使うことにしました。
（０）ファイルの読み込み
データファイルはこのプログラムと同じディレクトリに置いておくことにします。ファイル名も"GaAs
data"と言うことにしましょう。Current VI's Pathという機能を使うとそのプログラムのパス名が出力されます。プログラム名を切り取って"GaAs
data"をつなげるとよみたいファイルのパスができあがります。Read
From Spreadsheet File.viというアイコンにパスを入力（青緑の波線）すれば、ファイルからデータ(Raw
Data)を読み込んでデータが出力（オレンジ色の２重線）されます。データを次のシーケンスに渡します。（左上のボックス）
（１）回帰計算
左上のボックスから"Raw Data"が出てきます。左側に縦に並んだ２つの黄色いアイコンが配列の要素を取り出す機能を持っています。上のアイコンからはＸデータが、下のアイコンからはＹデータが取り出されます。青い
"n th Poly fit" と"Poly fit Alg" というアイコンから出ている線の先にあるのが"General
Polynomial Fit.vi" です。指定した次数で回帰を行い、入力されたＸデータに対応する回帰後のＹデータと回帰式の係数を出力します。係数は次数に対応する数だけ出力されるので、例えば２次回帰の場合は、３時、４次の係数としてゼロを結合します。（計算の次数にかかわらず大きさ５の１次元配列にします。大きさが変わらない方が後の処理は楽でしょ。）
（２）回帰曲線の計算
４次回帰はy = a0 + a1*x + a2*x^2 + a3*x^3 + a4*x^4 ですからx を次々入れてy
を求めます。すべて、ブロック図ではわかりにくくなる事もありますから、この例のように数式で演算することも可能です。Ｘの配列とＹの配列をバンドルという機能で、結合してｘ−ｙグラフに表示できる形式にします。
（３）グラフ表示
３本のプロットをビルドアレイと言う機能でまとめて、グラフに表示します。
戻る