2012-01-28から1日間の記事一覧

カットオフ

こちらの企画 親子鑑定 親子かそうじゃないか(親子ではなくて、『ただの関係』) 尤度の比をとる ルール・インする尤度比の基準がある ルール・アウトする尤度比の基準がある 尤度比検定 パラメタがあって、仮説がパラメタで表現されている 仮説ごとに尤度が…

DiagnosisMedパッケージを使ってROCカーブを描いてみる

こちらを参考 library(DiagnosisMed) Ncase<-100 Ncont<-100 mean.case<-40 sd.case<-3 mean.cont<-35 sd.cont<-2 Case<-rnorm(Ncase,mean.case,sd.case) Cont<-rnorm(Ncont,mean.cont,sd.cont) boxplot(cbind(Case,Cont)) ra<-range(c(Case,Cont)) breaks<-…

無性生殖集団の親子鑑定

倍々で増える無性生殖生物集団を考える すべての個体は不死であって、増殖を続けると仮定する DNA配列はポアッソン乱数で変異が入る Infinite modelを仮定し、変異は必ず新規なものとする グラフオブジェクトに納めることで グラフ上の距離(パスの長さ)が…