親子トリオ

# setf:父のディプロタイプの集合。アレルの集合(ホモなら要素数１、ヘテロなら要素数２)
# 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