Plutonium Biokinetics in Human Body A. Luciani - Kit-Bibliothek - FZK
Plutonium Biokinetics in Human Body A. Luciani - Kit-Bibliothek - FZK
Plutonium Biokinetics in Human Body A. Luciani - Kit-Bibliothek - FZK
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
calculat<strong>in</strong>g the amount of <strong>Plutonium</strong> activity <strong>in</strong> the blood compartment (often named transfer<br />
compartment too).<br />
3.1.1.2 Target function for data vs. model comparison<br />
The comparison of model predictions for a certa<strong>in</strong> variable (e.g; activity <strong>in</strong> an organ or<br />
excreted <strong>in</strong> bioassays) with available data or empirical curves can be carried out qualitatively,<br />
as a first attempt, by plott<strong>in</strong>g them. For a quantitative comparison it is necessary to select a<br />
target function that quantifies how good a certa<strong>in</strong> model fits a given data set best.<br />
An example of target function is given by Jones [111]. It is based on the assumption<br />
that bioassay data, specifically for ur<strong>in</strong>ary excretion, are distributed accord<strong>in</strong>g to a lognormal<br />
distribution [153]. It is def<strong>in</strong>ed as:<br />
where:<br />
d(t i) is the experimental datum at time t i;<br />
x(t i) is the value predicted by the model at time t i;<br />
n is the number of empirical data.<br />
[ ] 2<br />
n<br />
∑<br />
L = Logd(t i)− Logx(t i)<br />
i=1<br />
85<br />
equation 3.1.4<br />
Other functions can be designed for such purpose, consider<strong>in</strong>g that the more a target<br />
function po<strong>in</strong>ts out a disagreement between data and a curve, the more this function is<br />
appropriate. A target function was already used <strong>in</strong> a prelim<strong>in</strong>ary study to the present work<br />
[154]:<br />
F =<br />
n<br />
∑<br />
i=1<br />
⎡<br />
⎣<br />
⎢<br />
d(t i)−x(t i )<br />
x(t i )<br />
equation 3.1.5<br />
where the different elements have the same mean<strong>in</strong>g as <strong>in</strong> Jones’s target function. This target<br />
function has a similar structure to the χ 2 function. It is based on the deviation of each<br />
experimental datum from the relative value predicted by the model. The terms are then<br />
squared <strong>in</strong> order to elim<strong>in</strong>ate the sign effect of the deviation. For such purpose the absolute<br />
value operation, <strong>in</strong>stead of term squar<strong>in</strong>g, could be also adopted and the target function would<br />
be def<strong>in</strong>ed as:<br />
A =<br />
n<br />
∑<br />
i=1<br />
d(t i)−x(t i)<br />
x(t i )<br />
equation 3.1.6<br />
By compar<strong>in</strong>g the target function F and A <strong>in</strong> the equation 3.1.5 and equation 3.1.6, it<br />
can be seen that the F function will give more weight to large deviations. In fact if the<br />
experimental datum exceeds by two times the model’s prediction, the squared deviation is<br />
greater than the absolute value of the deviation. On the contrary if the experimental datum<br />
⎤<br />
⎦<br />
⎥<br />
2<br />
⎥