http://bekkoame.okwave.jp/qa2595561.htmlで質問させて頂いた者です。
y = γ / (1 + exp(α + β * x))という関数において,パラメータγ,α,βを推定したいのですが,初期値の設定がうまくいかなくて困っております。要するに初期値の探し方を教えて欲しいのです。
この関数からγ / y -1 = exp(α - β * x)という式が導けます。ここで両辺の自然対数をとると,
log(γ / y - 1) = α - β * x
となります。ここでY = log(γ / y - 1)とするとY = α - β * xとなって,これは線形モデルなので最小二乗法によって簡単にαとβの値が求められます。今回の場合だと得られているデータは,
x=[54,57,60,63,66,69,72]
y=[0,1,2,5,8,9,10]
なので,とりあえずγ=11と仮定して計算すると結局は
Y = [*, 2.30, 1.50, 0.18, -0.98, -1.50, -2.30]
となります(*はうまく計算できないので無視)。最終的にYのデータからY(ハット) = -0.949 * x + 3.1883という式が得られます。以前の質問でも補足しましたが,Rというデータ解析用のソフトウェアを使うと,
nls(y ~ gamma / ( 1 + exp(alpha + beta * x)),start = c(gamma = 11,alpha = 3.1883,beta = -0.949))
とするのですが,これは要するにy ~ gamma / ( 1 + exp(alpha + beta * x))の部分がモデル式,start = c(gamma = 11,alpha = 3.1883,beta = -0.949)の部分が初期値を設定する部分です。yとxという変数には上に記したように得られたデータがベクトルとして格納されています。
以上,説明が長くなってしまいましたが,今回のようにx=[54,57,60,63,66,69,72]となっていると,初期値が悪いとのエラーが出て計算ができません。これをx=[1,2,3,4,5,6,7]としてやるとうまくいきます。だから平行移動させるために前回のような質問をしたのですが、、、(ちょっと文字数に制限があるので意味不明だったら補足します)
No.1ベストアンサー
- 回答日時:
>y = γ / (1 + exp(α + β * x))という関数において
が、
>γ / y -1 = exp(α - β * x)
となってるのはなぜでしょうか?
γ / y -1 = exp(α + β * x)
では無いのですか?(勘違いだったら申し訳ないです)
早速の回答ありがとうございます。
おっしゃるとおり,γ / y -1 = exp(α - β * x)ではなくγ / y -1 = exp(α + β * x)です。私が間違えていました。
γ / y -1 = exp(α - β * x) ・・・×間違い
γ / y -1 = exp(α + β * x) ・・・○こちらが正しいです
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・14歳の自分に衝撃の事実を告げてください
- ・架空の映画のネタバレレビュー
- ・「お昼の放送」の思い出
- ・昨日見た夢を教えて下さい
- ・ちょっと先の未来クイズ第4問
- ・【大喜利】【投稿~10/21(月)】買ったばかりの自転車を分解してひと言
- ・メモのコツを教えてください!
- ・CDの保有枚数を教えてください
- ・ホテルを選ぶとき、これだけは譲れない条件TOP3は?
- ・家・車以外で、人生で一番奮発した買い物
- ・人生最悪の忘れ物
- ・【コナン30周年】嘘でしょ!?と思った○○周年を教えて【ハルヒ20周年】
- ・ハマっている「お菓子」を教えて!
- ・最近、いつ泣きましたか?
- ・夏が終わったと感じる瞬間って、どんな時?
- ・10秒目をつむったら…
- ・人生のプチ美学を教えてください!!
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・都道府県穴埋めゲーム
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
エクセルで60進法計算の仕方...
-
エクセル関数で源泉徴収額を計...
-
100以下の自然数のうち、次のよ...
-
公務員試験の問題
-
常用対数と自然対数の違い
-
整式を整理する・計算する
-
一般結合法則
-
「量子超越」とは、どういうこ...
-
文字を含む三角関数の定積分の...
-
経費率の計算方法を教えて下さい。
-
微分方程式について
-
工事の共通仮設費率の計算がで...
-
10分の1は「10/1 それとも1/10...
-
実績を積むという表現
-
50以下は“50”も入るのですか?
-
1億x1億はいくらでしょうか?
-
【機械図面】 最大値・最小値...
-
dBm/HzからdBm/MHzへの単位変換
-
5進法を10進法への直し方
-
言葉遣いについて ○○を取りに行...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
100以下の自然数のうち、次のよ...
-
エクセル関数で源泉徴収額を計...
-
エクセルで60進法計算の仕方...
-
この産み分けの計算でハズレの...
-
工事の共通仮設費率の計算がで...
-
Excelで勤務の過不足時間を計算...
-
経費率の計算方法を教えて下さい。
-
中3の問題です
-
3桁の自然数の中で、次の個数を...
-
最小公倍数と最大公約数の求め...
-
エクセルでのシグマ計算
-
ラプラス変換の「s」とは?
-
端数を習うのは小学何年生の頃...
-
関数電卓の使い方
-
問1.絶対値が3より小さい整数に...
-
99の10乗の下位5桁の数を求める...
-
中学受験の算数の問題
-
(3)、(4)の解説をお願いし...
-
乱数について
-
文字を含む三角関数の定積分の...
おすすめ情報