- 目次はこちら
- こちらでm次元分割表の正確確率検定の話、そして、そのp値計算のために、与えられた周辺度数を満足するテーブルの集合があった
- そして、その集合を動き回るための基底としてのマルコフ基底が登場した
- ここで、「テーブル集合」はマルコフ基底によって、網羅的に動き回ることができる
- この閉じていること、動き回れることが「代数的」なので、Algebraic statisticsの内容としてこの話題が登場している
- この記事は1. Markov Basesの後半
- 2次元分割表の「動き」は
- 行の和がゼロ、列の和がゼロになるような動き
- それを「単位」に分解すると、2行、2列を取り出し、分配方法を取り換える、という動きが、最小単位
- すべての動きはその組合せで表せる
- ここで言う動きは、
が同じ周辺度数を持つ表(を1次元ベクトルにしたもの)の関係を定める行列
を使って考えれば、
を満足するような
の例になっている
- m次元分割表の「動き」は
を満足する
に属する。
(Aのカーネル・核空間・零空間)→Wiki
- ここで、きちんと書くために、表記法を書いておく
は、m次元
であるようなテーブルのすべてのセルの値の和が
であるようなテーブルに対応したベクトルの集合
- このようなテーブル集合の1要素
に対して
は「
によって
と定義される「モデル
」に関する
次元分割表
のfiber」
- 言い換えると、総サンプル数
なるm次元分割表の集合を扱う。その集合の1例
があったときに、
と
を用いて、ある関係になるようなm次元分割表の集合が
のfiber。ただし、
はなんでもよいわけではなくて、制約を入れておく。この場合には、生起確率の積でうまくいくような関係に対応させた定義を用いており、そのことの意味は、m次元が相互に独立な確率変数であることである。また、そのことは、定義を表現するときに、対数の1次線形和になるので、log-linearモデルと呼んでいる。
,
を使った「動き」の定義とは
もそのfiber
も表の形になったベクトル(の集合)
- これらを足し算で移動したい。その足し算に制約を持たせる
- 足し算の「足される側」と「足す側」と、その「答」のすべてがこの集合に含まれること(和について閉じている)
- 足し算を要素の積み重ねとして表したときに、その途中の経路も集合に含まれている
- これを数式で書くと
に関して、
があったとき
,
を満足するような
を要素とする有限部分集合
がMarkob basis と呼ばれ、
は
の(定める)「動き」の要素と呼ばれる
- その例として
次元分割表の場合には、2次元間の取り換え、3次元間のぐるり、…、m次元間のぐるり、を集めたものもある(これは、2,3,...,m次の置換群)になるが、これは格子基底(lattice bases)と呼ばれMarkov bases に含まれるらしい…
- こんな感じの「動きの基底」であるが、定義を広げたり縮めたりすることでいろいろある
- lattice basis
Markov basis
Grobner basis
universal Grobner basis
Graver basis
- lattice basis
数学・コンピュータ関連の姉妹ブログ『ryamadaのコンピュータ・数学メモ』
京都大学大学院医学研究科ゲノム医学センター統計遺伝学分野のWiki
講義・スライド
医学生物学と数学とプログラミングの三重学習を狙う学習ツール
駆け足で読む○○シリーズ
ぱらぱらめくるシリーズ
カシオの計算機
オンライン整数列大辞典