「遺伝統計学」の編集履歴（バックアップ）一覧に戻る
遺伝統計学 - (2007/07/24 (火) 19:54:24) のソース
=項目=
*[[アーミテージ検定]]
*[[ベイズ因子による検定]]
*[[確率統計]]
*[[確率分布]]
*[[Coalescent_theory]]
=その他=
==<math>\lambda_s</math>==
*Rybicki and Elston. AJHG 2000;66:593.
*<math>\lambda_s</math>はRisch 1990のレビュー三枚により使われるようになった。このとき、多重ローカスモデルではλは各ローカスのλの積で表される（<math>\lambda=\lambda_1 \lambda_2 ... \lambda_i</math>)ことを示した。
*またある状況では罹患同胞対連鎖解析の検出力を予測できることを示したが、これは否定されたよう
*欠点: 共有する遺伝子と同じくらい共有する環境因子の影響を受ける。
*ジェノタイプ相対リスクを<math>\gamma</math>とする。ジェノタイプ相対リスクは、そのジェノタイプの人の発症率を、それ以外のジェノタイプの発症率で割ったものである。したがって、g+をsusceptible genotypeありの状態、g-を無しの状態とすると、
::<math>\gamma=\frac{Pr(D|g+)}{Pr(D|g-)} \quad (1)</math>
:と書ける。一方<math>\lambda</math>はaffected relativeをARと書くことにより、
::<math>\lambda=\frac{Pr(D|AR)}{Pr(D)} \quad (2)</math>
:と書ける。ここで集団でのsusceptible genotypeの頻度を<math>p_G</math>、それ以外を<math>p_0</math>とすると、
::<math>Pr(D)=p_G Pr(D|G+)+p_0 Pr(D|G-) \quad (3)</math>
:である。(2)式と(3)式より、
::<math>Pr(D|G+)=\frac{Pr(D|AR)\div\lambda-p_0Pr(D|G-)}{p_G} \quad (4)</math>
:である。(4)式を(1)式に代入すれば、
::<math>\gamma=\frac{Pr(D|AR)\div\lambda-p_0Pr(D|G-)}{p_G Pr(D|g-)} \quad (5)</math>
:たとえば、<math>\lambda_s=5</math>で、集団の発症率が1/1,000である疾患と感受性ジェノタイプを考える。感受性ジェノタイプがない場合の発症率を1/10,000、集団の90%が感受性ジェノタイプを持っていないとするなら、(5)式より<math>\gamma</math>は91となり、<math>\lambda_s</math>の20倍近い。
==sequence coverage==
* raw value for coverage:
:<math>\frac{L+T}{R}</math>
ここでRはリファレンスであるSNPセット（たとえばHapMap SNP）、Tはtag SNPセット(tagSNPと、それによってまとめられるSNP)である。ここで、ゲノム全体のSNPセット<math>G>R</math>であるためこの式ではcoverageを過大に評価してしまう。したがって
:<math>\frac{\left ( \frac{L}{R-T} \right ) (G-T)+T}{G}</math>
ENCODEの予測から、<math>G \approx</math>7500万と示唆されているが、これが1億だからといって予測が大きく変わるわけではない。
==r<sup>2</sup>とpower==
*真のsusceptible locusであるlocus 1でcase-control N<sub>1</sub>サンプルをタイピング、近傍のマーカーlocus 2でN<sub>2</sub>サンプルをタイピングし関連研究をしたとする。
**locus 1のアレルはAとa、locus2のアレルはBとb。
**<math>\pi_{DA}</math>はケースのAアレル頻度、<math>\pi_{CA}</math>はコントロールのAアレル頻度。
**q<sub>AB</sub>はある染色体上でlocus1がアレルA、locus2がアレルBである確率。
*すると、（論文ではつづめてあったがバカな自分のために展開）
::<math>\pi_{DB}-\pi_{CB}=(\pi_{DA}-\pi_{CA})q_{AB}+(\pi_{Da}-\pi_{Ca})q_{aB}</math>
::<math>=(\pi_{DA}-\pi_{CA})q_{AB}+(1-\pi_{DA}-1+\pi_{CA})q_{aB}</math>
::<math>=(\pi_{DA}-\pi_{CA})(q_{AB}-q_{aB})</math>
*するとlocus1での関連解析の[[遺伝用語#カイ二乗|カイ二乗統計量]]は、ケースのサンプル数:コントロールのサンプル数=φ:1-φとすると、
::<math>\chi^2_1 = \frac{\{(\phi 2N_1+(1-\phi) 2N_1\}\{\phi 2N_1\hat{\pi}_{DA}(1-\phi)2N_1\hat{\pi}_{Ca}-\phi 2N_1\hat{\pi}_{Da}(1-\phi)2N_1\hat{\pi}_{CA}\}^2}{\phi 2N_1 (1-\phi) 2N_1 \hat{\pi}_A 2N_1 (1-\hat{\pi}_A)}</math>
::<math>= \frac{(\hat{\pi}_{DA}-\hat{\pi}_{CA})^2 2N_1 \phi (1-\phi)}{\hat{\pi}_A(1-\hat{\pi}_A)}</math>
::ここで、<math> \hat{\pi}_{DA}, \hat{\pi}_{CA}, \hat{\pi}_A</math>はそれぞれ患者、コントロール、全体でのアレルAのサンプル頻度である。AをBにしてN<sub>1</sub>をN<sub>2</sub>にすればlocus 2の<math>\chi^2_2</math>となる。
*<math>\chi^2_1</math>と<math>\chi^2_2</math>の分布はだいたい正規分布の確率変数の二乗である。正規分布の平均は
::<math>(\pi_{DA}-\pi_{CA})\left \[ \frac{2N_1 \phi(1-\phi)}{\bar{\pi}_A(1-\bar{\pi}_A)}\right \] ^{\frac{1}{2}}</math>
::と
::<math>(\pi_{DA}-\pi_{CA})(q_{AB}-q_{aB})\left \[ \frac{2N_2 \phi(1-\phi)}{\bar{\pi}_B(1-\bar{\pi}_B)}\right \] ^{\frac{1}{2}}</math>
::で
::<math>\bar{\pi}_A=\phi \pi_{DA}+(1-\phi) \pi_{CA} \approx \pi_A</math>
*ここで、（下の式は自分で導出できず）
::<math>r^2=\frac{(q_{AB}-q_{aB})^2\pi_A(1-\pi_A)}{\pi_B(1-\pi_B)}</math>
*したがって、N<sub>2</sub>=N<sub>1</sub>/r<sup>2</sup>のとき<math>\chi^2_2</math>は<math>\chi^2_1</math>とほぼ同じ分布となる。
=Journal abstracts=
*[[Altshuler_and_Daly_2007]]
