Learn R Programming

PhylogeneticEM (version 1.8.0)

sum_simplex: Sum on a simplex

Description

sum_simplex returns the sum on k_1+...+k_p=K, k_i>0 of the products of NN[i,k_i].

Usage

sum_simplex(NN, K, p)

Value

double : the result of the sum.

Arguments

NN

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

K

an integer. The number of groups wanted.

p

an integer. The number of daughters of a node.

Details

This function uses xsimplex to perform the product of NN[i,k_i] for all combination of k_i such that k_1+...+k_p=K, k_i>0, using function prod_index. Then sum all the products.