ある画像から輝度ヒストグラムを取ったところ3つのガウス分布が重なっているような形状になりました。
これをサンプルデータとし、それを3つのガウス分布の和で近似しようと考えています。
求めたいパラメータはガウス分布A・B・Cの高さh1・h2・h3、平均値μ1・μ2・μ3、分散σ1・σ2・σ3です。
最小二乗法により解こうと考えましたが9元連立方程式となり、指数関数の内外にパラメータが入っており、1つの方程式自体長くなっていることにより解くことが困難です。
パラメータを求める上での条件は3つです。
・サンプルデータを1つのガウス分布で近似しその差を取り、さらにその差をガウス分布で近似し差を取り、再びガウス分布で近似するといった方法を取らない
・h1・h2・h3は必ず正の値をとる
・パラメータはあらかじめおおよその値を予想できない
これらの条件の下で解くことは可能でしょうか。
もし可能ならその解き方を教えてください。
この時、いくつかの解法の仕方がある場合はできるだけ簡単かつ計算量の少ないものを教えてください。
よろしくお願いします。
No.1ベストアンサー
- 回答日時:
この問題は、典型的な非線形最適化問題となります。
解析的に解くことはできません。
3つのガウス分布の和の形で与えた式g(x)と実測値f(x)の間の
平均2乗誤差を評価関数とするとよいでしょう。
ただし、パラメータを変換しておきましょう。
(1)h1,h2,h3≧0のためには、代わりにs1,s2,s3を使い
h1 = (s1)^2, h2=(s2)^2, h3=(s3)^2
(2)σ1、・・・の方は、2乗するので、そのままにしましょう。
これで、拘束条件なしの最適化問題になりました。
評価関数は、
J(s1,s2,s3,μ1,μ2,μ3,σ1,σ2,σ3)
=∫|f(x)-g(x|s1,・・・,σ3)|^2 dx
≒Σ_i{|f(x_i)-g(x_i|s1,・・・,σ3)|^2}Δx
x_iはヒストグラムのサンプル点です。
Jの導関数を使う方法と使わない方法があります。
使わない方法では、Hooke&Jeeves法があります。
(Hooke Jeevesで検索すると、Cのプログラムが見つかる
かも)
これは、パラメータを1つずつ±に少し変化させ、減少する
方向へ移動していくことを原理とする方法です。パターン
サーチ法のひとつです。
導関数を使う場合は、パラメータベクトルをwとすると、
w(n+1)=w(n)-∇_w J(w(n))
いずれにしても、初期値は適当に決める必要があります。
まったく分からない場合は、ランダムに決めざるを得ません。
繰り返し近似法ですから、計算量の評価は難しい。
なお、パラメータ群を推定することが目的ですか
または、関数近似が目的ですか。
パラメータ群の推定が目的ならば、パラメータの唯一性
(Identifiability)の証明が必要と思います。
あなたの条件には合いませんが、ご参考までに
やはり解析的に解くことはできませんか・・・
目的としては3つのガウス分布に近似しベイズの決定則により最適な閾値を決定することです。
まずは導関数を使う方法でやってみたいと思います。
わかりやすいアドバイスありがとうございました。
No.2
- 回答日時:
混合正規分布のパラメータを決める方法として、EMアルゴリズムを使う方法もあります。
繰り返し計算により、パラメータを最尤推定する方法ですが、詳しいところは、参考URLをご覧ください。
この方法も初期パラメータを適当に決めないといけませんが、EMアルゴリズムの場合かなり柔軟に決めていいそうです。
参考URL:http://www.neurosci.aist.go.jp/~kurita/lecture/s …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
二項母集団の母比率の区間推定
-
精密法について教えてください...
-
3σについて
-
二項分布と正規分布の違い?
-
成績順位の推定
-
2つのサイコロの目の合計が偶数...
-
4つのサイコロを同時に振って...
-
重複順列
-
統計学の問題なんですが・・・
-
五分を6回連続で外すのはなん...
-
数学の質問です。 一枚の硬貨を...
-
確率0.02%って10000人に2人です...
-
一般常識を教えてください。1割...
-
30%の確率が5回連続で起きない...
-
数学Aの確率の問題です(´;ω;...
-
1個のサイコロを3回続けて投げ...
-
数学の問題で4C0の答えを教えて...
-
じゃんけんの問題
-
サイコロを三つ同時に振るとき ...
-
3つのサイコロを振って出た目の...
おすすめ情報