Statistics (Introduction to Radiochemistry)

From mn/safe/nukwik
Jump to: navigation, search

In statistics data is often characterized by being exact and reproducible. An exact result is close to the true value, while a reproducible result gives information about how precise the measurements are. However a reproducible data set does not necessarily imply that the result is exact. The average of a data set is given by:

[math]\bar{x} = \frac{x_{1}+x_{2}+...+x_{N}}{N}[/math] eqn 1

Where N is the number of measurements.

To describe a data set the variance or standard deviation must be specified as well. This shows the spread of the values in the data set around the average value. The variance (σ2) is given by;

[math]\sigma ^{2} = \frac{(x_{1}- \bar{x} )^{2} +(x_{2}- \bar{x} )^{2} +...+(x_{N}- \bar{x} )^{2}}{N-1}[/math] eqn 2

and the standard deviation (σ) is given by:

[math]\sigma = \sqrt{\sigma ^{2}}[/math] eqn 3

A coin toss has two different outcomes. For instance the probability of 10 heads in a row can be calculated as (1/2)100.00098, where every coin toss is independent of each other. This is a binary process that can be described by a binomial probability distribution function (which gives the probability for each possible result of the measurement). 

Over a period of time a given radioactive nucleus has 2 outcomes. It can either disintegrate or remain unchanged and, as with a coin toss, can be classified as a binary process. Within the half-life of a nuclide, half of the observed atoms will have undergone disintegration, exactly as the distribution of an exceedingly large number of coin tosses would be 0.5 head and 0.5 tail. The binary process applies to whether or not a given radioactive particle (α,β,...) will be registered by a detector as well.

Radioactive decay usually involves systems with so many atoms, that calculation of the probability for individual atoms is rather impractical. When the number of atoms is large (N>>1) and the time of observation is short compared to the half-life (σt<<1) Poisson distribution is normally used. Poisson distribution describes rare random events. P(x) is the probability to get a certain amount of counts x when [math]\bar{x}[/math] is the expected average value.

[math]P(x) = \frac{\bar{x}^{x} e^{7-\bar{x}}}{x!}[/math] eqn 4

The number of events N is large and [math]\bar{x} - x \lt \lt \sqrt{N}[/math]

The middle value and standard deviation is given by: 

[math]\bar{x} = \sum_{x=0}^N xP(x) = pN[/math] eqn 5


[math]\sigma = \sqrt{pN(1-p)} \approx \sqrt{pN} = \sqrt{\bar{x}}[/math] eqn 6

p the time t.


Normal distribution (Gaussian distribution) is another common probability function. It can be used as a simplification when x is large and [math]\bar{x} - x \lt \lt \bar{x}[/math]. Gaussian distribution is symmetrical around [math]\bar{x}[/math]and [math]\sigma \approx \sqrt{\bar{x}}[/math]applies to the standard deviation.

In a series of measurements where the middle value x is given with uncertainty σ the interval of x - σ and x + σis called the confidence interval. The probability of the true value to be within the confidence interval [math]x\pm \sigma[/math] is 68.3%. The table below shows how the probability of the true value to be within the confidence interval increases with the number of σ.

Interval around the measurement number Probability (%)
[math]\pm 0.674\sigma[/math] 50.0
[math]\pm \sigma[/math] 68.3
[math]\pm 1.96 \sigma[/math] 95.0
[math]\pm 3 \sigma[/math] 99.7

 

Calculations with values that each have uncertainty is often conducted, e.g subtraction of the background from measurements. The following rules apply for calculations in the experimental part:

  • Multiplication/division with a number c without uncertainty:

[math] x\pm \sigma gives \left\{ \begin{matrix} x_{c}=x \cdot c \\ \sigma_{c} = \sigma \cdot c \end{matrix} \right. [/math]

  • Addition/subtraction of a number with uncertainty:

[math] \left\{ \begin{matrix} x_{1} \pm \sigma_{1} \\ x_{2} \pm \sigma_{2} \end{matrix} \right. \rightarrow \left\{ \begin{matrix} x=x_{1} + x_{2} \\ \sigma = \sqrt{\sigma_{1}^{2} + \sigma{2}^{2}} \end{matrix} \right. [/math]