まずは2次元の正規分布を描画する。簡便のために対角成分以外は0とする。 import numpy as np from scipy.stats import multivariate_normal import matplotlib.pyplot as plt np.random.seed(2) mu_x,mu_y = 5,10 sigma_x,sigma_y = 1,5 cov = np.array([[1…

正規分布について考えるため、以下の図で記号を定義する。 出典: PDG 信頼区間(confidence interval) を定義したとき、有意度(σ)と、信頼区間から外れる確率(両側の場合はα、片側の場合はα/2)は以下のように計算できる。 from scipy.stats import norm for s…