﻿import numpy as np  # モジュールNumpyをnpという名前でインポート
lambda_0 = 1.75  # λの真のパラメータ値：λ0
np.random.seed(1)  # 乱数のシード値を固定
# ポアソン分布Po(λ0)に従う乱数を10000個生成し，xに代入
x = np.random.poisson(lam=lambda_0, size=10000)
# ポアソン分布の最尤推定値
lambda_hat = np.mean(x)
print("λの最尤推定値:", lambda_hat)