Intra species factors calculation

There is variation between individuals concerning their individual sensitivities to experience health effects. In some scenarios the aim is to perform assessments for the sensitive individuals instead of the average individuals for which the points of departure are derived. If this is the case, then extrapolation is required to translate hazard characterisations derived for the average individual to hazard characterisations for a sensitive individual.

Traditionally a fixed safety factor describes the variation between individuals. Little information is available about the true distribution of human sensitivities, so there is also a large uncertainty. In MCRA, the intra-species variability is modelled explicitly using a lognormal distribution, characterised by a geometric mean (GM) equal to 1 and a geometric standard deviation (GSD). This distribution is used to sample individual hazard characterisations. This effectively converts the description of hazard characterisations to include variability, with an unbiased central value. GM is 1 by definition (50% of the population is assumed to be less sensitive and 50% more sensitive than the average individual) and has no uncertainty. On the other hand, there is uncertainty about the GSD or SD = log(GSD). This uncertainty is described by a chi-square distribution with df degrees of freedom:

dfSD2σ2χdf2

where σ is the true standard deviation and χdf2 denotes a chi-square distribution with df degrees of freedom. Thus df characterises the amount of uncertainty regarding the intra species variation quantified by (G)SD.

It is difficult to specify values for GSD and df. A practical way to do this is, define a p95-sensitive individual having an intra-species factor corresponding to the p95 percentile of the distribution describing the variability.

Assume that the p95-sensitive individuals are between e.g. 2 and 10 times more sensitive than the average human. The values 2 and 10 are to be interpreted as uncertainty bounds, the 2.5th and 97.5th percentiles of the uncertainty distribution for the Fintra,p95.

Then, p95-sensitive individuals are between Fintra,p95p2.5 and Fintra,p95p97.5 more sensitive than the average human, where the first index indicate the variability distribution and the second the uncertainty distribution.

For the lognormal distribution we have:

ln(Fintra,p95)=ln(GM)+1.645σ=1.645σ

And thus σ=ln(Fintra,p95)1.645. The uncertainty statement can be re-expressed as a 95% confidence interval for σ:

ln(Fintra,p95p2.5)1.645σln(Fintra,p95p97.5)1.645

A standard two-sided 95% confidence interval for σ derived from the χ2 distribution is:

SDdfχdf,0.9752σSDdfχdf,0.0252

Equating the two lower bounds as well as the two upper bounds we obtain two equations with two unknowns. Substituting for SD and rearranging we find that df can be found from:

χdf,0.9752χdf,0.0252=(ln(Fintra,p95p97.5)ln(Fintra,p95p2.5))2

And then GSD can be calculated as:

GSD=exp(SD)=(Fintra,p95p97.5)χdf,0.0252df/1.645

A bisection algorithm can be used to find the value for df for which the chi-square ratio in the example is

χdf,0.9752χdf,0.0252=(ln(10)ln(2))2=11.035

This value is attained for df = 6.25. See also van der Voet et al. (2009)