al<-matrix(0,36,2) al[,1]<-c(rep(2,9),rep(1,12),rep(0,15)) al[,2]<-c(2,2,1,1,1,0,0,0,0,2,2,2,1,1,1,1,0,0,0,0,0,2,2,2,2,1,1,1,1,1,0,0,0,0,0,0) sr<-al[,1]*1+al[,2]*1 P<-sr+rnorm(36,0,0.1) rgl.clear() rgl.bg(col="white") x<-cbind(P,al) d<-dis…

遺伝率0.5で、リスク０とリスクが±1くらい（正規分布でちょっとランダムに揺らす）にしたらこうなった。リスク0ではパフォーマンスがよくなくてほっとした。次は遺伝率を変えてみよう

リスクのあたえかたを正規分布にしてみた。右は平均0で、左は平均が1のものと-1のものを組み合わせたものになる。RVにリスクがある時は検出されることがわかったけれど、RVにリスクがない時検出されないかが分からない。今遺伝率１なのでリスク0にするとフェ…

200人、バリアント80個、ただし200人について観測されるRV(mafフェノタイプの値を、検定の時に、-0.5~0.5くらいの間になるようにした。最大値と最小値の絶対値の和で割る。