確率分布が正規分布に従うか調べたい、 二つの集団が同じ確率分布から得られたものか調べたい、 といった時に使うのが、コロモゴロフスミルノフ検定(Kolmogorov–Smirnov test) コルモゴロフ–スミルノフ検定（コルモゴロフ–スミルノフけんてい、英: Kolmogoro…

二項分布の試行回数を無限大に大きくしていくと正規分布に近づくことが知られている。 しかし、その証明は意外と知られていない。（中心極限定理でも証明は可能ではあるが、回りくどすぎて初学者には不親切）やが出てくるのがどうしても不思議で、自分なりに…

あるディレクトリの直下のファイルに対して巡回しながら処理したい時 シェルスクリプトとpythonでどうするか。 bash files="./hoge/*" for filepath in $files; do echo $filepath done ./hoge/1.txt ./hoge/2.txt のようにディレクトリごと表示される ファ…

以前の記事の続きとして、せっかくなので友愛数・完全数・婚約数を求めてみたい。 約数関数（再掲） divisor(n)として約数関数を使うので、再掲しておく。 sympyを使う。 import sympy def divisor(n): factors = sympy.factorint(n) rst = 1 for i,j in fac…

pythonのライブラリの一つsympyを使えば、簡単に素因数分解ができるということを知った。 import sympy sympy.factorint(1000) #{2: 3, 5: 3} ちなみに、因数分解も簡単にできる!! 凄い。 import sympy x = sympy.Symbol('x') eqn = x**2 - 3*x + 2 print(sy…

一元配置分散分析とは 「3つ以上の群があった時に、果たしてそれらの群の平均は等しいと言えるかどうか。」 という検定。 集団の分布が正規性を持つことが前提となっている。*1 ※すべての組み合わせペアでt検定を適用するのは間違いなので注意*2基本的な発想…

Flaskで開発していたり、http通信が必要となるjsのライブラリを使っているとlocalhostのwebサーバーを立てたくなることがある。 そこで苦労の末辿り着いた便利エイリアスを紹介したい。 それがこれ alias httpserver="python -m SimpleHTTPServer& start chr…

a = [1,2,3,2,3,4,5,4,3,2,1] こんなデータがあった時、3がある場所を取ってきたいとする。 a.index(3) # 2 とすると最初のインデックスだけが返ってくる。 希望は全てを取ってくることなので、これではダメ。 じゃあどうするか。 index = [] for i,j in enu…