compute_sum_var_diff computes sum_{e edge} ell_j * Var[X_j - X_pa(j) | Y]
compute_sum_var_diff
compute_sum_var_diff(phylo, var_diff)
phylo
a phylogenetic tree
var_diff
result of function compute_var_diff.BM
compute_var_diff.BM
matrix p x p