F分布のヒストグラム

第一自由度3、第二自由度5のF分布で10000回抽出した場合のヒストグラムを示します。それぞれの母集団は、次の通りです。

1 6 1 8 0 3 3 9 8 8 7

2 4 1 4 2 1 3 5 6 2 3

ソースコード

using Plots
using StatsBase

x = [1 6 1 8 0 3 3 9 8 8 7]
y = [2 4 1 4 2 1 3 5 6 2 3]

values = []

for i ∈ 1:10000
    push!(values, (rand(x, 3) .^ 2 |> sum |> x -> x / 3) / (rand(y, 5) .^ 2 |> sum |> x -> x / 5))
end

histogram(
    values,
    title="F-distribution", label="F(3, 5)",
    linecolor="#4521A1", linewidth=2, fillcolor="#BFE35E", fillalpha=0.2,
    xlims=(0, maximum(values))
    )

savefig("f-distribution.png")

ヒストグラム