親子トリオ
# setf:父のディプロタイプの集合。アレルの集合(ホモなら要素数1、ヘテロなら要素数2) # setm:母のディプロタイプの集合 # seto:子のディプロタイプの集合 # 父母カップル LikeTrio<-function(setf,setm,seto){ couple<-setf*setm cnt<-0 for(i in couple){ if(as.set(i) == seto)cnt<-cnt+1 } cnt/length(couple) } setf<-set(1,3) setm<-set(2,3) seto<-set(3) LikeTrio(setf,setm,seto)
> setf {1, 3} > setm {2, 3} > seto {3} > LikeTrio(setf,setm,seto) [1] 0.25