# diversi.gof

##### Tests of Constant Diversification Rates

This function computes two tests of the distribution of branching
times using the `z`

): in this case, `x`

and `z`

must
be of the same length. See the examples for how to compute the latter
from a sample of expected branching times.

- Keywords
- univar

##### Usage

`diversi.gof(x, null = "exponential", z = NULL)`

##### Arguments

- x
- a numeric vector with the branching times.
- null
- a character string specifying the null distribution for
the branching times. Only two choices are possible: either
`"exponential"`

, or`"user"`

. - z
- used if
`null = "user"`

; gives the expected distribution under the model.

##### Details

The

The distributions of both test statistics depend on the null hypothesis, and on whether or not some parameters were estimated from the data. However, these distributions are not known precisely and critical values were determined by Stephens (1974) using simulations. These critical values were used for the present function.

##### Value

- A NULL value is returned, the results are simply printed.

##### encoding

utf8

##### References

Paradis, E. (1998) Testing for constant diversification rates using
molecular phylogenies: a general approach based on statistical tests
for goodness of fit. *Molecular Biology and Evolution*,
**15**, 476--479.

Stephens, M. A. (1974) EDF statistics for goodness of fit and some
comparisons. *Journal of the American Statistical Association*,
**69**, 730--737.

##### See Also

`branching.times`

, `diversi.time`

`ltt.plot`

, `birthdeath`

, `yule`

,
`yule.cov`

##### Examples

```
data(bird.families)
x <- branching.times(bird.families)
### suppose we have a sample of expected branching times `y';
### for simplicity, take them from a uniform distribution:
y <- runif(500, 0, max(x) + 1) # + 1 to avoid A2 = Inf
### now compute the expected cumulative distribution:
x <- sort(x)
N <- length(x)
ecdf <- numeric(N)
for (i in 1:N) ecdf[i] <- sum(y <= x[i])/500
### finally do the test:
diversi.gof(x, "user", z = ecdf)
```

*Documentation reproduced from package ape, version 3.2, License: GPL (>= 2)*