FAQ
FAMOS
シーケンス / 関数top
└数値
ID.
Q. 決定係数を求めたい。
A.
直接求める関数は用意されてないため、下記の様なシーケンスで求めて下さい。 近似波形の係数を求めるために使用しているApproNonLin()関数はimc FAMOS Ver 7.2以降のProfessional版以上が必要です。
例)
LOAD sintest1
_data = sintest1 ; サンプルデータ
_A = ApproNonLin(_data,"A1*sin(A2*PI2*x)-y") ; 近似波形の係数
_approx = _A[1] * sin(_A[2] * PI2 * ramp(0, XDel?(_data), Leng?(_data))) ; 近似波形
_numerator = Sum((_data - _approx)^2) ; ①分子: 残差の二乗和
_denominator = Sum((_data - Mean(_data))^2) ; ②分母: 標本値の平均からの差の二乗和
R2 = 1 - _numerator / _denominator ; 決定係数: ①を②で割ったものを1から引いた値
; カーブウィンドウに表示
CwNewWindow("Cv1", "show")
CwSelectWindow("Cv1")
CwAction("reset")
CwPosition(0,0,640,480)
CwDisplaySet("title", "")
CwNewChannel("append last axis", _data)
CwNewChannel("append last axis", _approx)
CwNewElement("marker.text")
CwMarkerSet("x.type", 2)
CwMarkerSet("x", 73)
CwMarkerSet("y.type", 2)
CwMarkerSet("y", 80)
CwMarkerSet("font.size", 12)
CwMarkerSet("text", "決定係数＝" + TForm(R2, "F1.3"))