﻿# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add(["Distributions","Statistics"])
using Distributions, Statistics
import Random  # 乱数を扱うためのモジュールをインポート
lambda_0 = 1.75  # λの真のパラメータ値：λ0
Random.seed!(1)  # 乱数のシード値を固定
# ポアソン分布Po(λ0)に従う乱数を10000個生成し，xに代入
x = rand(Poisson(lambda_0), 10000)
# ポアソン分布の最尤推定値
lambda_hat = mean(x)
println("λの最尤推定値：", lambda_hat)