# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add("Distributions") 
import Random  # 乱数を扱うためのモジュールをインポート
using Distributions
p0 = 0.65  # pの真のパラメータ値：p0
Random.seed!(1)  # 乱数のシード値を固定
# 幾何分布Ge(p0)に従う乱数を10000個生成
x = rand(Geometric(p0), 10000)
p_hat = 1/(1 + mean(x))  # 最尤推定値を計算
println("p_hat = ", p_hat)  # 最尤推定値を表示