高校数学の質問スレ(医者・東大卒専用) Part438 (882レス)
上
下
前
次
1-
新
848
: 2025/05/25(日) 05:57:52.88
ID:P4nhnL8B(2/6)
調
AA×
[240|
320
|
480
|
600
|
100%
|
JPG
|
べ
|
レス栞
|
レス消
]
848: [sage] 2025/05/25(日) 05:57:52.88 ID:P4nhnL8B >847のRのコードをChatGPTで Mathematicaにコメント付きで移植 (* betaParameter 関数: 指定された信頼区間 [L, U] に、指定された信頼度 credMass(例: 95%)の確率質量を持つ ベータ分布のパラメータ α, β を算出する。 *) betaParameter[L_: 1/7, U_: 1/5, credMass_: 0.95] := Module[ {α, β}, (* f[α] は、PDF[BetaDistribution[α, β], L] == PDF[BetaDistribution[α, β], U] を満たすように β を α に基づいて算出する関数。 *) f[α_] := 1 + ((α - 1) * Log[U / L]) / Log[(1 - L) / (1 - U)]; (* g[α] は、ベータ分布 Beta[α, f[α]] の区間 [L, U] に 含まれる確率(CDFの差)を返す関数。 *) g[α_] := CDF[BetaDistribution[α, f[α]], U] - CDF[BetaDistribution[α, f[α]], L]; (* g[α] = credMass を満たす α を数値的に求める *) α = α /. FindRoot[g[α] == credMass, {α, 1, 1*^5}]; (* 対応する β を算出 *) β = f[α]; (* 結果を返す *) {α, β} ] (* 関数を実行して α, β を取得 *) {α, β} = betaParameter[] (* g[α] を評価して、[L, U] に credMass の質量があることを確認 *) g[α] http://rio2016.5ch.net/test/read.cgi/math/1723152147/848
ののコードをで にコメント付きで移植 関数 指定された信頼区間 に指定された信頼度 例 の確率質量を持つ ベータ分布のパラメータ を算出する は を満たすように を に基づいて算出する関数 はベータ分布 の区間 に 含まれる確率の差を返す関数 を満たす を数値的に求める 対応する を算出 結果を返す 関数を実行して を取得 を評価して に の質量があることを確認
上
下
前
次
1-
新
書
関
写
板
覧
索
設
栞
歴
あと 34 レスあります
スレ情報
赤レス抽出
画像レス抽出
歴の未読スレ
ぬこの手
ぬこTOP
0.039s