【メモ書き】pythonでロジット変換・オメガ変換・プロビット変換

この記事で各種線形変換を行ったので、忘れないようにpythonのコードも記載しておく。

 

#データの作成
data = np.random.normal(loc=25, scale = 5, size=100000) #平均25,標準偏差5のサンプル100000個
x = np.sort(data)#
y = 1. * np.arange(len(data)) / (len(data) - 1) #正規分布の累積分布関数

y_p = stats.norm.ppf(y) #プロビット変換(正規分布の逆関数)

y_l = np.log(y/(1-y)) #ロジット変換

y_o = 10 * np.log10(y/(1-y)) #オメガ変換

 

f:id:Chemstat:20210227102007j:plain