■

フェノタイプとジェノタイプに関係があるかを、遺伝率で調整してもいいのだけど、リスクとアレル頻度に絞って考えることにする。たとえばリスク=0では絶対にフェノタイプとジェノタイプに関係性がなくなるので。

リスクを、risk=exp(-10000*af)にしてみると、検出力が弱くなった。どのくらいの大きさでどうなるのかを調べて、そのあとでアレル頻度の分布を変えつつやってみる。
→0/1の方と同じで、あまりにも極端な値があるとそれに引っ張られて検出力が落ちている気がする。アレル頻度の分布を考えながらやろう。

2変数のときはどうするのが分かりやすいだろうか。。。
シミュレーションを走らせつつ考える。

リスクを高くしたときこんな感じ