A 回答 (4件)
- 最新から表示
- 回答順に表示
No.1
- 回答日時:
この回答への補足
丁寧に答えて頂き、本当にありがとうございます。
教えていただいた方法で解こうとしたのですが,
下の形だと、微分したときにln中の分母のcやfが消えてしまうのですが、
どのようにしたらよいのでしょうか?
lnY((1/L)lnA) = aln((x-b)/c) - dln((x-e)/f) (Aは定数です。)
あと、解いて中でxで微分をして、その後にxで積分していますが、
これはどのような意味なのでしょうか?
また、最小二乗法をもっと勉強してみたいので、なにか良い参考書やURLを
ご存知でしたら教えていただきたいと思います。
大変自分勝手に質問してしまって申し訳ありません。
もしお時間がありましたらご返事を頂きたいとおもいます。
No.2
- 回答日時:
> 丁寧に答えて頂き、本当にありがとうございます。
きついなあ~。手抜きですいませんねえ。
ええとですね。(x,y)の観測値=サンプルが沢山得られているからこそ最小二乗法の対象になるわけですが、このご質問の場合の観測値(x,y)は式のどの部分に該当するんでしょうか。さらにどれが定数なのか。Pも定数?それともこれがy?どうもはっきりしません。
仮にLとPが定数だとすれば、これはxに関する方程式ですから最小二乗法とは縁もゆかりもない。ですからPが実はサンプル?
という訳で、どの変数についてサンプルが得られるのかを明記して戴かないと詳細なことは申し上げられません。
aln((x-b)/c) - dln((x-e)/f)=aln((x-b)) - dln((x-e)) +(-alnc+dlnf)
と書けますから、xで微分すると定数項(-alnc+dlnf)が消滅します。消滅しちゃって構いません。で、ご紹介したURLの方法でa,b,d,eが決定できてから、それを元の式
ln{1/L*ln(1/(1-P))} = aln((x-b)) - dln((x-e)) +(-alnc+dlnf)
に入れてみれば、未知数なのは(-alnc+dlnf)だけ。つまり線形最小二乗法で(-alnc+dlnf)が決定できますね。
これとは別の話として、cとfがそれぞれ幾らであるかは決められません、どんな方法を使おうとも。(-alnc+dlnf)が幾らであるか、だけしか決まらない。なぜなら、a,d,(-alnc+dlnf)=vが決まっているときcに好きな正の数を代入すれば、f = exp[(v+alnc)/d] と決まる。つまりcとfの一方は好きな値にできるからです。
> xで微分をして、その後にxで積分していますが、 これはどのような意味なのでしょうか?
原理的には微分しただけでも構わないんですが、それだと数値微分を計算しなくてはなりません。数値微分はサンプルのでこぼこの影響を受けやすくて、実用上は無理がある。そこで積分して滑らかになるようにしているんです。
この方式はstomachmanが開発し勝手に「超平面法」と称しているもので、一般的に良く知られた方法という訳ではありませんから、もしご質問の目的がレポートだったりする場合には、先生が目を白黒させるかも。
そのうち論文にしたろか、とも思うんですが、どういう場合には手が出ないのかなど、理論の方がまだ纏まりません。
過去のQ/Aから「最小二乗法」を検索なされば、バリエーションが幾つか見つかります。
この回答への補足
本当にありがとうございます。
ln{1/L*ln(1/(1-P))} = aln((X-c)/d)-bln((e-X)/f)
この式で(x,y)にあたるのはXとPなんです。
Pは確率で、0~1の間で変化するんだろうと思います。
これはレポートではないんです。
私は学生で、行なっている実験の解析にこの最小二乗法が必要なんです。
でも、手も足もでませんでした。
それにしても、stomachmanさんがどんなことをやってらっしゃるひとなのか
存じませんが、すごいですね。本当にすごいです。
No.3
- 回答日時:
もし差し支えなければ、ですけれど。
どんな実験で出てくるモデルなのか、それぞれの変数・定数が何を意味するのか、ちょっとご紹介願えないでしょうか。この回答への補足
はい。
私は今、繊維一本(ガラス繊維や炭素繊維など)の強度について
色々と研究しています。
繊維は長くなるほど、その強度は小さくなっていき、
この式はその長さと強さの関係を表している式なんです。
ln{1/L*ln(1/(1-P))} = aln((X-c)/d)-bln((e-X)/f)
Lは繊維の長さ、Xが実験での強度が入ります。
c,eがそれぞれある範囲での強度の最小値と最大値です。
このc,eは定数とみてもいいそうです。
a,bはshape parameterとなっており、形状母数にあたると思います。
d,fはscaling parameterとなっており、尺度母数にあたると思います。
Pは、その負荷で繊維が切れてしまう確立です。
この式から実験を使って、長さと強度の関係のグラフを求めたいのですが、
まだまだ勉強不足でわかりません。
もしお時間がありましたら、アドバイスを頂きたいと思います。
よろしくお願いいたします。
No.4
- 回答日時:
> Pは、その負荷で繊維が切れてしまう確立です。
なるほど。そういうことでしたら、ご紹介したURLのやり方よりも、まじめに非線形最小二乗法をやる方が良さそうです。
●c<X<eである。ln(X-c)とln(e-X)が両方定義できるためにはこれは必要十分ですね。
切れない確率をQ=1-Pとするとき
ln(Q) =-L {((X-c)/d)^a ((e-X)/f)^-b}
つまり
Q=exp[-KL]
これはなんとなく納得できます。深さがある分布に従うランダムな傷があって、或る程度以上のがLに含まれない確率。で係数Kが
K={((X-c)/d)^a} / {((e-X)/f)^b}
= g {(X-c)^a} / {(e-X)^b}
でやんすか。(ここにg=(f^b)/(d^a)としました。既に申し上げたように、このモデルからはf,dをそれぞれ幾らと決めることはできない。)
そうして、c<X<eである。このc,eは定数。つまりP(あるいはQ)とは関係なくc,eが分かっているんでないと、なんだか話がおかしいっすよね。
しかし、X→eになったらK→∞だからQ→0つまり絶対切れる。X→cになったらK→0、つまりQ→1ですから、断じても切れない。だからXは強度と言うより「弱度」あるいは負荷っすね。
●さて、これはPが確率ですから、一個のサンプル(X[k],P[k])を取るだけでも何度も実験を繰り返さないといけません。従って沢山のサンプルが期待できるものではない。そうなると、ご紹介したURLのように「データが非常によくモデルと合っていて、しかも沢山ある」という条件を前提にしたやり方は適当でないでしょう。
だからむしろ非線形最小二乗法をやる方が良さそうです。(教科書はいつも「最小二乗法による実験データ解析」(東京大学出版会)をお勧めしています。)しかし、非線形最小二乗法の場合、パラメータの大体の推定値が予め得られていることが必要です。
●この問題の場合、c,eを既知の定数とした時には、めちゃ簡単っす。
C=-a(ln d)+b(ln f)
として
ln{1/L*ln(1/(1-P))} = aln(X-c)-bln(e-X) +C
ここで負荷X[k]を与えたときの切れる確率をP[k]と書くとき、
y[k] = ln{1/L*ln(1/(1-P[k]))}
u[k] = ln(X[k]-c)
v[k] = -ln(e-X[k])
と定義すれば、
y[k] = a u[k] + b v[k] + C
ですからパラメータa,b,Cを含む単なる線形最小二乗法の問題。
●これでa,b,Cを決めれば、c,e はだいたいは分かっているんですから、
y[k] = aln(X[k]-c)-bln(e-X[k]) +C
という非線形最小二乗法を正攻法でやっつける時の初期値が決まったことになります。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
- ・ゆるやかでぃべーと タイムマシンを破壊すべきか。
- ・歩いた自慢大会
- ・許せない心理テスト
- ・字面がカッコいい英単語
- ・これ何て呼びますか Part2
- ・人生で一番思い出に残ってる靴
- ・ゆるやかでぃべーと すべての高校生はアルバイトをするべきだ。
- ・初めて自分の家と他人の家が違う、と意識した時
- ・単二電池
- ・チョコミントアイス
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
2階微分d^2y/dx^2を詳しく教え...
-
log(1+x)の微分
-
3階微分って何がわかるの??
-
y^2をxについて微分してください
-
分母が文字の分数を微分する方...
-
数学についての質問です。 1.f(...
-
sinx^2の微分って2xcosx^2であ...
-
サイン二乗xの微分を教えてく...
-
授業で「yをxで微分する」とい...
-
二次関数 y=x^2 を微分すると---
-
d^2y/dx^2は何と読めばいいので...
-
三角関数 微分
-
y=1-sinxを微分せよ
-
二回微分して 上に凸下に凸 が...
-
数学の微分の問題です。 xのsin...
-
微分
-
e^(logx+1)の微分
-
f(x)=0はxで微分可能か
-
これらの数式を声に出して読む...
-
y=1/(2x-1)を微分する方法につ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
2階微分d^2y/dx^2を詳しく教え...
-
サイン二乗xの微分を教えてく...
-
二回微分して 上に凸下に凸 が...
-
3階微分って何がわかるの??
-
授業で「yをxで微分する」とい...
-
分母が文字の分数を微分する方...
-
y^2をxについて微分してください
-
y=e^x^x 微分 問題
-
これらの数式を声に出して読む...
-
和積・積和の公式について質問...
-
log(1+x)の微分
-
三角関数の微分の問題なんです...
-
二次関数 y=x^2 を微分すると---
-
なぜ微分したら円の面積が円周...
-
不定積分の計算で出た定数は捨...
-
位置を微分したら速度?
-
lim[x→0](e^x - e^-x)/x
-
dxやdyの本当の意味は?
-
y=logxA(Aは定数)をxで微分
-
数3の「eのh乗引く1をでh割...
おすすめ情報