総数を固定しなくても… メモ は (スターリングの公式)より であって に注意したり、,に注意したりして変形すると この式の移行関係には、対数尤度比がカイ自乗分布に従うという話がある k <- 5 library(MCMCpack) f <- c(rdirichlet(1,runif(k))) N <- 10000…

２項分布の連続化：β分布 多項分布の連続化：ディリクレ分布 確率ベクトル、観測度数, 確率 に制約をつけつつ、確率を最大にするようなが知りたいとする Rのoptim()関数で最適化することにする 少し普通ではない変数の置き方をする 長さのベクトルを考える …

３つの多次元分布がある X,Y,Zとする 標本がたくさんあって、X,Y,Zが形成されるが、標本は３分布を横串で貫いているのでX,Y,Zの間には何かしらの関係がありえる X,Yの関係とX,Zの関係は既知として、Y,Z間の関係があるのかないのかをデータから読み取りたいと…

カテゴリカルな尺度が３軸ある X,Y,Zの３つ XとYが関係し、XとZが関係している YとZはXとの関係とは独立にさらに関係しているのかどうか知りたい ３軸が作る３次元表が観察されたとする X,Yの２軸が作る２次元カテゴリの比率は偏りがあるが、そのままにして …

CAST cohort allelic sum test collapses information on all rare variants within a region C-alpha test non-burden-based method robust to the direction and magnitude of effect not easy to adjust for covariates, incl. population stratification…