計算機科学のブログ

知識を広げる 確率と統計 平均(期待値)、中央値、3つのサイコロの目の合計、二項分布、リスト、線型フィット、直線、2次曲線

ハンズ・オン・スタートMathematica® -Wolfram言語™によるプログラミング (C・ヘイスティング(著)、K・ミショー(著)、M・モリソン(著)、ウルフラム・リサーチ(翻訳)、丸善出版)の第2部(知識を広げる)、第18章(確率と統計)の練習問題1、2、3、4、5、6、7、8、9、10の解答を求めてみる。

mean of 93, 86, 68, 94, 91, 88, 74

Mean[{93, 86, 68, 94, 91, 88, 74}]
Output
N[%,4]

84.86
N[%,2]

85.
steps1 = Table[i - 1, {i, 0, 5}];

Median[steps1]
Output
N[%]

1.5
Probability[x == 3,
            Distributed[x, BinomialDistribution[4, 8/10]]
]
Output
N[%]

0.4096
probs = Table[
    Probability[x == result,
                Distributed[x, BinomialDistribution[4, 8 / 10]]
    ],
    {result, 0, 4}
]
Output
BarChart[probs, ChartLabels -> Range[0, 4]]
Output
labresults1 = Table[{2i-1,i^2-5i+1}, {i, 1, 25}];

linearfit1 = LinearModelFit[labresults1, {1, x}, x]
Output
quadfit1 = LinearModelFit[labresults1, {1, x, x^2}, x]
Output
vislabresults1 = ListPlot[labresults1];

vislinearfit1 = Plot[linearfit1[x], {x, 0, 50}, PlotStyle -> Red];
visquadfit1 = Plot[quadfit1[x], {x, 0, 50}, PlotStyle->Green];

Show[
    visquadfit1,
    vislinearfit1,
    vislabresults1
]
Output
rolls = RandomChoice[Range[1, 6], 100]

{3, 1, 1, 2, 3, 4, 3, 4, 2, 3, 2, 3, 4, 5, 3, 2, 6, 4, 4, 3, 3, 6, 4, 3, 5, 5, 3, 6, 5, 
 
>   6, 4, 5, 5, 4, 3, 1, 2, 6, 3, 1, 1, 4, 4, 1, 6, 3, 1, 2, 1, 3, 6, 1, 3, 1, 6, 5, 3, 
 
>   4, 3, 1, 6, 2, 4, 5, 5, 2, 5, 3, 4, 4, 5, 6, 3, 4, 3, 4, 2, 4, 5, 4, 4, 3, 2, 2, 5, 
 
>   2, 3, 6, 2, 4, 1, 6, 1, 3, 5, 1, 4, 1, 5, 4}
Histogram[rolls]
Output
rolls = RandomChoice[Range[1, 6], 100];
Histogram[rolls]
Output
rolls = RandomChoice[Range[1, 6], 100];
Histogram[rolls]
Output
rolls = RandomChoice[Range[1, 6], 100];
Histogram[rolls]
Output
rolls = RandomChoice[Range[1, 6], 100];
Histogram[rolls]
Output
Probability[
    x + y + z == 11 || x + y + z == 12,
    Distributed[x, DiscreteUniformDistribution[{1, 6}]] &&
    Distributed[y, DiscreteUniformDistribution[{1, 6}]] &&
    Distributed[z, DiscreteUniformDistribution[{1, 6}]]
]
Output
N[%]

0.240741