Rのimage()関数を使ってペアワイズＬＤプロットを簡単に描く方法を以前書いた(こちら)。 今日は、ハプロタイプの01表記ファイル"hoge.txt"から描くことにする。 横がマーカー、縦に染色体。こんな感じ。 0 0 0 0 1 0 1 1 0 1 0 1 1 0 0 1 1 1 0 0 1 0 0 0 1 …

数え上げ # 順列 permN<-function(N=10,k=3){ return(exp(lgamma(N+1)-lgamma((N-k)+1))) } 組み合わせ combN<-function(N=10,k=3){ return( exp(lgamma(N+1)-lgamma((N-k)+1)-lgamma(k+1)) ) } 重複順列 repPermN<-function(N=10,k=3){ return(N^k) } 重複…