18.11.2012 Views

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

SHOW MORE
SHOW LESS

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 />

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!