# ケースとコントロールの人数 Nca<-Nco<-1000 # ケースコントロール別のレアアレル保有者数 nca<-20 nco<-0 # 計算するための、ケース・コントロールのアレル頻度(パラメタとして振ってある) dp<-0.001 pca<-seq(from=dp,to=0.05,length=100) pco<-pca # 確…