ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ R シミュレーション 家系 1. アレル、ジェノタイプ アレルとアレル頻度 # あるマーカー # アレルの数 Na<-4 # アレル頻度 Pa<-c(0.3,0.2,0.1,0.4) # Pa<-c(3,2,1,4)のようにして、あとで頻度の和が１になるようにしてもよい # 念のため、アレル頻度の和が１になるように補正する Pa<-Pa/sum(Pa) アレル頻度を適当に決めたいとき # 一様乱数からの発生 Pa<-runif(Na) # アレル頻度が1/Na付近に固まるので… # そうならないように道具をウェブ上から導入する install.packages("MCMCpack") library(MCMCpack) … 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20120116%2F1324171876" title="家系データをシミュレーションで作る１ - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/r/ryamada22/20120113/20120113091222.jpg Hatena Blog https://hatena.blog 2012-01-16 10:31:16 rich https://ryamada22.hatenablog.jp/entry/20120116/1324171876 1.0 100%