This command can infer some average edge lengths and assign them from a (bootstrap/MCMC) sample.

## Arguments

- tree
a phylogenetic tree or splitnetwork where edge lengths are assigned to.

- trees
an object of class multiPhylo, where the average for the edges is computed from.

- fun
a function to compute the average (default is median).

- rooted
rooted logical, if FALSE edge lengths is a function of the observed splits, if TRUE edge lengths are estimated from height for the observed clades.

## Examples

```
data("Laurasiatherian")
set.seed(123)
bs <- bootstrap.phyDat(Laurasiatherian,
FUN=function(x)upgma(dist.ml(x)), bs=100)
tree_compat <- allCompat(bs, rooted=TRUE) |>
add_edge_length(bs)
plot(tree_compat)
add_boxplot(tree_compat, bs, boxwex=.7)
```