OLSシリーズ第4弾。 OLSは本当に重要なのでしつこくやる。今回は、scipyのstats.linregressを使う。いつも通り下記の回帰式を推定する。コードは以下の通り。 # -*- coding: utf-8 -*- import numpy as np from scipy import stats #サンプリング数 N=1000 #…

対数尤度が合わないとAICやBICといった値も必然的に合わなくなると思われるので、けっこう困った話になる。 今回はずれるという事象だけをメモしておく。気が向いたらなぜずれるのかを検証してみるかもしれない。回帰式は例によってこれ。Pythonで作ったデー…