Y=AX+Bという形で(Y,X)の組から係数A,Bを推定する方法が最小二乗法かと思います。(そのような整理が妥当であるか、意味があるか、についてはデータの分布形状や変動係数などで別途検証する必要がありますがここでは不問として)
さて、その次ですが、YY=AX+BY+CZという関数で表示しようと考えたとします。(YY,X,Y,Z)というデータ系列があり、係数(A,B,C)を求めるというような解析です。(前述のようにそれが妥当かどうかは不問として。)係数A,B,Cを最小二乗法のようにできるだけ妥当に求める方法を教えていただきたいのですが。
1)たとえば、YY=AX+CCとして最小二乗法でAを求めたあと、残差CCについてCC=BY+CZとおいて再度最小二乗法でBを求めるというような手法を何度も行うとか? これだと解く順番に依存するようなので全部の順番でやって平均値を取るとか? 下手な考え休むに似たりのようですが。
2)カルマンフィルタ、ニューラルネットワークの手法で解析するということもあるでしょうか。その場合、YY=AX+BY+CZという枠組みではなくなりますが、それでもよいと言えばよいです。
このような問題を解く方法について解説されている本とか分野の名称(数理統計学とか?)などを教えて頂けると探索する手がかりとなります。今は手がかりすらはっきりしていないのです。
最小二乗法についてはその程度だったら専門分野ではなく、各分野で個別に解説されているようで込み入った問題になったときにどこを捜索したらよいかわからないもので質問しました。よろしくお願いします。
A 回答 (4件)
- 最新から表示
- 回答順に表示
No.4
- 回答日時:
解き方の計算法は
#3さんの方法で良いですね。
実際に最小二乗法を解くのは手計算では大変ですので
エクセルやfree-softの「gnuplot」の中の関数fitting機能
を使います。
僕は後者のソフトを利用しています。データのファイルを読み込んで簡単にfittingによりA,B,Cを計算してくれます。
参考URL
http://ayapin.film.s.dendai.ac.jp/~matuda/Gnuplo …
gnuplotダウンロード先
↓
参考URL:http://www.gnuplot.info/
回答ありがとうございます。
gnuplotはインストールしてあります。これはグラフ作成ツールという風に位置づけておりますが、このような計算機能もあるということですね。関数フィッティング機能ということは裏でそのような作業をやっているということ推察できます。エクセルもそうだと思いますが。
他にもたとえばRとかmaximaなどのフリーソフトにもそのような機能がありそうですね。
このようなフリー系の数理解析ソフト群はかなり進んできたように思います。ただ、使いこなすのが大変で1つに秀でたらかなり役に立つとは思いますが、目先があちこちに動いてなかなか習得できないでいます。
No.3
- 回答日時:
Y(i)とA*X(i)+Bの差の2乗和を最小にするのがY=AX+Bという形での最小二乗法ですよね。
そのためにf(A,B)=Σ(Y(i)-A*X(i)+B)^2のA,Bに関する偏微分を計算し、
df/dA=0とdf/dB=0を連立させてA,Bを求めるとOKというものです。
YY=AX+BY+CZという形でもやることは同じです。
g(A,B,C)=Σ(YY(i)-A*X(i)-B*Y(i)-C*Z(i))^2として、A,B,Cに関する偏微分が0になる点を求めます。
なるほど。
gをA,B,Cでそれぞれ偏微分してゼロとなるようなA,B,Cを求めるということですね。いわれて見ればそういうことかと思いますが。
これは非線形となると話が違ってくるでしょうか。x^2, xy, y^2というような項が出てくるということですが。展開する関数に直交性があれば、それでも良いのかなとは思いますが。
回答ありがとうございました。
No.2
- 回答日時:
参考URLを書き忘れていました
ソルバーはすべてのエクセルにもとから実装してある機能ではないので、使用されるときは参考URLを参考にしてください
参考URL:http://homepage1.nifty.com/gfk/square_solver.htm
No.1
- 回答日時:
私は以前はYY=AX+BY+CZなどの関数の係数を求める解析にはエクセルのソルバー機能を使って求めていました。
ソルバー機能はある目的のセルの値が最小になるように変数をエクセルが自動的に変化させ最適な値を求めてくれるという機能です。
f(x)=AX+BY+CZとしてA,B,Cをまず適当に定めて、各X,Y,Zにおいて(YY-F(X))^2を求めその合計値が最小になるようにすると最適なA,B,Cが求まったことになりますよね
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
0から1になった時の増加率を教...
-
自然対数Ln(x)からxを求める方...
-
9X2乗-6X+1 はどうやった...
-
素数定理では満足できない
-
z=(-1-√3i)/2 とするとき、
-
数学 ∑(1からnまで)1/k2乗...
-
数学の大学入試問題の解き方
-
中学数学 a※b=1/3(a+b)とする...
-
バスケの得点、フィボナッチ数...
-
1/3乗などの計算方法
-
漸化式の問題です。
-
3次関数f(x)=x~3-9x~2+12x-1で...
-
(x^4-2x^3-4x^2+13x-2)/(x-1)^5...
-
反復計算で指数方程式の解を求...
-
指数計算の確認
-
20〜200までの自然数の和
-
逆関数の求め方
-
複素関数ー留数についての質問です
-
不定積分の答えをどこまで出す...
-
1-1+1-1+…=sqrt(2)って証明でき...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報