All functions |
|
|---|---|
An example dataset |
|
Generalized Hansen models |
|
Estimate ancestral states given a fitted OUwie model |
|
Parametric bootstrap function |
|
Generates data for contour plot of likelihood surface |
|
Generalized Detection of shifts in OU process |
|
Generalized Hansen model likelihood calculator |
|
Format data and tree for OUwie |
|
Generalized Hansen model simulator |
|
A test of regime identifiability |
|
Dents the likelihood surface This takes any values that are better (lower) than the desired negative log likelihood and reflects them across the best_neglnL + delta line, "denting" the likelihood surface. |
|
Propose new values This proposes new values using a normal distribution centered on the original parameter values, with desired standard deviation. If any proposed values are outside the bounds, it will propose again. |
|
Sample points from along a ridge This "dents" the likelihood surface by reflecting points better than a threshold back across the threshold (think of taking a hollow plastic model of a mountain and punching the top so it's a volcano). It then uses essentially a Metropolis-Hastings walk to wander around the new rim. It adjusts the proposal width so that it samples points around the desired likelihood. This is better than using the curvature at the maximum likelihood estimate since it can actually sample points in case the assumptions of the curvature method do not hold. It is better than varying one parameter at a time while holding others constant because that could miss ridges: if I am fitting 5=x+y, and get a point estimate of (3,2), the reality is that there are an infinite range of values of x and y that will sum to 5, but if I hold x constant it looks like y is estimated very precisely. Of course, one could just fully embrace the Metropolis-Hastings lifestyle and use a full Bayesian approach. |
|
Adjust tree for matrix condition |
|
Model average the parameter estimates over severl hOUwie fits. |
|
Generate a table from a set of hOUwie models describing their relative fit to data. |
|
Generate a continuous model parameter structure |
|
Fit a joint model of discrete and continuous characters via maximum-likelihood. |
|
Fit a joint model of discrete and continuous characters via maximum-likelihood with fixed regimes. |
|
Reconstruct the marginal probability of discrete node states under the hOUwie model. |
|
Simulate a discrete and continuous character following a Markov and Ornstein-Uhlenbeck model. |
|
Rerun a set of hOUwie models with the best mappings of the set. |
|
Sample points from along a ridge for a hOUwie model |
|
Contour plot |
|
Plot the dented samples This will show the univariate plots of the parameter values versus the likelihood as well as bivariate plots of pairs of parameters to look for ridges. |
|
Print dentist print summary of output from dent_walk |
|
Summarize dentist Display summary of output from dent_walk |
|