核家族ごとにディプロタイプ総当たり条件の確率・尤度を計算する
LikeNucWithPrior<-function(cpnout,nucs,DiplotypePrior){ prob<-list() for(nn in 1:length(nucs)){ nucDimVector<-cpnout[[nn]][[2]] tmp<-DiplotypePrior[[nucDimVector[1]]] for(i in 2:length(nucDimVector)){ tmp<-tmp%o%DiplotypePrior[[nucDimVector[i]]] } #print(dim(cpnout[[nn]][[1]])) #print(dim(tmp)) #print("---") prob[[nn]]<-tmp * cpnout[[nn]][[1]] } prob }