﻿# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add(["Distributions","Statistics"])
using Distributions, Statistics
import Random  # 乱数を扱うためのモジュールをインポート
m0 = -2.5  # μの真のパラメータ値：μ0
s0 = 1.5  # σの真のパラメータ値：σ0
Random.seed!(1)  # 乱数のシード値を固定
# 正規分布N(μ0,σ0*σ0)に従う乱数を10000個生成
x = rand(Normal(m0, s0), 10000)
# 最尤推定の実行
normdist = fit(Normal, x)
mean_hat = mean(normdist)
sd_hat = std(normdist)
# 推定結果の表示
println("mean = ", mean_hat)
println("sd = ", sd_hat)