ハプロタイプシミュレーション

Dirichlet分布、ぼーず確率

ぼーず確率については、ryamadaの遺伝学・遺伝統計学メモ - Coalescence発生待ち時間の算出で書いたが・・・。その発展編。 Dirichlet分布一般については、関連記事(こちらとこちら)も参照。 ながらくイメージのつかめなかった、Dirichlet分布がようやくわか…

Coalescent発生待ち時間の算出理論

参考資料:19.6 Coalescent theory pp591 by C. Neuhauser in "Handbook of Statistical Genetics Second Edition" by D.J.Balding, M. Bishop and C. Cannings (2003) Wiley Wiley::Handbook of Statistical Genetics, Second Edition 仮定 ある世代時点Tに…

ボーズ(一匹も釣れない)になる確率の指数関数exp(-x)表現

「(釣り人が)ボーズ」になる確率の説明404 Not Found 単位時間あたり、m匹の魚を釣り上げる釣り人が、t単位時間の間0匹しか釣果がない確率prob_Bを「ボーズ確率」と呼んでいる。ボーズとは、一匹も釣果がないこと(のようだ)。今、mが時間によらず一定だとす…

GBM(geometric Brownian model)※

Log-normal growthとも呼ばれる 特徴 Rateがランダムである ランダムではあるが、ランダム加減の程度は「現在のrate」に比例する Recombination rateのをある染色体の場所におけるrecombination rateとのlog ratio で表したとき、その2点間の物理距離と相関…

ケース・コントロール用ハプロタイプデータのシミュレーション生成に関して盛り込むパラメタ

参考文献no title Coalescent, recombination, diffusion, 変異のselection力, 人口増減, recombination rateのばらつきを盛り込んで、疾患責任変異とその周辺ハプロタイプをシミュレートし、そのハプロタイプに対して、phenocopy, 低penetranceを考慮して、…