候補者がたくさんいて、網羅的に探すときはメンデリアンチェックで絞り込む
library(paramlink) tmpPed<-p[,1:5] tmpPed<-data.frame(ID=tmpPed[,1],FID=tmpPed[,3],MID=tmpPed[,2],SEX=tmpPed[,4]+1,AFF=tmpPed[,5]-1) tmpPed<-linkdat(tmpPed,model=1) # 候補者を当てはめて確定したジェノタイプがtmpG2 tttmpG<-tmpG2[,,1] for(i in 2:length(tmpG2[1,1,])){ tttmpG<-cbind(tttmpG,tmpG2[,,i]) } tttmpG tmpPed2=setMarkers(tmpPed, m=tttmpG) mendelianCheck(tmpPed2)