As part of a Bayesian method to estimate the divergence times of species, priors have to be set with values based on previous literature or known fossil dates. These priors can have different distributions selected to model the uncertainity.
For example, if I had prior information that a particular species diverged between 5 and 10 million years ago, I might want to model this with a normal distribution with mean 7.5, the 5% quantile at 5, and the 95% quantile at 10.
However, I need to input the mean and standard deviation (SD) to describe this. Is there a way, other than trial and error, to find the SD of a normal distribution using only the mean and values for the 5% and 95% quantiles?
q.95then in R
sd=(q.95-q.05)/(2*qnorm(0.95))calculates the standard deviation; similar commands would apply in other languages. $\endgroup$