高校数学の質問スレ Part434

725 ：１３２人目の素数さん：2024/04/30(火) 20:38:09.28 ID:VcpWQbIP.net

>>683
>レベルが下がる　そのまま　上がる　となりそれぞれに確率が設定されています
の確率に関しては情報がないため
下がる　そのまま　上がる　の確率は　形状パラメータ（１，１，１）のディリクレ分布に従って変動するとして計算する。

乱数発生させてWolfram言語でのシミュレーション（推敲希望）

sim[] :=(
item=0;
L=0;
While[L<10,
p1p2=RandomVariate[DirichletDistribution[{1,1,1}]];
p={p1p2[[1]],p1p2[[2]],1-Total[p1p2]};
L=L + RandomChoice[p -> {-1,0,1}];
item++];
item
)
試行回数に上限がないの算出までに時間がかかる。

出力例

In[20]:= items=Table[sim[],100]

Out[20]= {134, 1452, 108, 256, 427, 137, 258, 817, 38, 191, 33, 1340, 21084, 74730, 201, 106, 2523, 2909, 623, 2024,

> 26, 74, 246, 203, 5135, 4473, 536, 6742, 1341, 171, 22, 144, 115, 61, 32, 90, 88, 697, 105, 120, 21503, 355,

> 26018, 15051, 199, 18576, 936, 194, 531, 801, 1457, 90, 114, 104787, 3017, 434, 176, 1180, 494, 144, 1411, 358,

> 25, 1960, 429, 129997, 1960, 8345, 364, 1185, 356, 190, 139, 301, 149814, 547, 132, 458, 12, 231, 1351170, 17175,

> 981, 353, 136, 104657, 7607, 18538, 1621, 265, 923, 260, 58, 768, 1141, 180, 122, 197, 112, 78}


summary(items)
Min. 1st Qu. Median Mean 3rd Qu. Max.
12.0 136.8 361.0 21304.2 1498.0 1351170.0

uncertainty interval(分位数で算出）

In[25]:= Quantile[items,{0.025,0.975}]

Out[25]= {25, 129997}