Setting relative bounds

[David Bass]

Is there a way to set relative parameter bounds during inference, e.g. such that the time of one population split is always greater than that of another? There does not appear to be a way to do this in the core scipy optimization functions. I am aware that I can simply discard inference runs that return estimates with unsatisfactory relative values of estimates, but I'm hoping that there is a way to do so while speeding up computation by telling the optimization algorithm to ignore the region of parameter space in which a certain parameter is less than another parameter.

Thanks,

David

[Ryan Gutenkunst]

Hello David,

Using the dadi.Inference.opt method, which is a wrapper around NLopt, you can specify inequality constraints among your parameters (ineq_constraints) to define functional relationships between the parameters you would like to preserve during optimization.

Best,
Ryan

> --
> You received this message because you are subscribed to the Google Groups "dadi-user" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to dadi-user+...@googlegroups.com.
> To view this discussion on the web visit https://groups.google.com/d/msgid/dadi-user/908b7741-0ed4-4389-9147-ea335b88d76bn%40googlegroups.com.