prod_index {PhylogeneticEM} | R Documentation |
prod_index
return the product of chosen elements of a matrix.
prod_index(X, Id)
X |
a matrix with p rows and K column. Each row contains the number of partition in 1<=k<=K groups for one of the p children of a given node |
Id |
a vector of length p, result of the function |
This function is to be used in sum_simplex
to be applied to all the
elements given by xsimplex
.
Performs the product : X[1,Id[1]]*X[2,Id[2]]*...*X[p,Id[p]] if all the
elements of Id are positive. Otherwise, just return 0.
double : the result of the product.