# Forest ErrAna StatDist

## Parent Distribution

Let represent our ith attempt to measurement the quantity

Due to the random errors present in any experiment we should not expect .

If we neglect systematic errors, then we should expect to, on average, follow some probability distribution around the correct value .

This probability distribution can be referred to as the "parent population".

## Average and Variance

### Average

The word "average" is used to describe a property of a "parent" probability distribution or a set of observations/measurements made in an experiment which gives an indication of a likely outcome of an experiment.

The symbol

is usually used to represent the "mean" of a known probability (parent) distribution (parent mean) while the "average" of a set of observations/measurements is denoted as

and is commonly referred to as the "sample" average or "sample mean".

Definition of the mean

Here the above average of a parent distribution is defined in terms of an infinite sum of observations () of an observable x divided by the number of observations.

is a calculation of the mean using a finite number of observations

This definition uses the assumption that the result of an experiment, measuring a sample average of , asymptotically approaches the "true" average of the parent distribution .

### Variance

The word "variance" is used to describe a property of a probability distribution or a set of observations/measurements made in an experiment which gives an indication how much an observation will deviate from and average value.

A deviation of any measurement from a parent distribution with a mean can be defined as

the deviations should average to ZERO for an infinite number of observations by definition of the mean.

Definition of the average

the deviation average would be given by

But the AVERAGE DEVIATION is given by an average of the magnitude of the deviations given by

= a measure of the dispersion of the expected observations about the mean

Taking the absolute value though is cumbersome when performing a statistical analysis so one may express this dispersion in terms of the variance

A typical variable used to denote the variance is

and is defined as

#### Standard Deviation

The standard deviation is defined as the square root of the variance

S.D. =

The mean should be thought of as a parameter which characterizes the observations we are making in an experiment. In general the mean specifies the probability distribution that is representative of the observable we are trying to measure through experimentation.

The variance characterizes the uncertainty associated with our experimental attempts to determine the "true" value. Although the mean and true value may not be equal, their difference should be less than the uncertainty given by the governing probability distribution.

#### Another Expression for Variance

Using the definition of variance (omitting the limit as )

Evaluating the definition of variance

In the case where your observables are weighted by the probability of the observation one has an expectations value given by

and a variance given by

to compare with the previous expression, let

## Average for an unknown probability distribution (parent population)

If the "Parent Population" is not known, you are just given a list of numbers with no indication of the probability distribution that they were drawn from, then the average and variance may be calculate as shown below.

### Arithmetic Mean and variance

If observables are made in an experiment then the arithmetic mean of those observables is defined as

The "unbiased" variance of the above sample is defined as

If you were told that the average is then you can calculate the

"true" variance of the above sample as

= RMS Error= Root Mean Squared Error
Note
RMS = Root Mean Square = =

#### Statistical Variance decreases with N

The repetition of an experiment can decrease the STATISTICAL error of the experiment

Consider the following:

The average value of the mean of a sample of n observations drawn from the parent population is the same as the average value of each observation. (The average of the averages is the same as one of the averages)

sample mean
if all means are the same

This is the reason why the sample mean is a measure of the population average ( )

Now consider the variance of the average of the averages (this is not the variance of the individual measurements but the variance of their means)

If the measurements are all independent

The above is the essence of counting statistics.

It says that the STATISTICAL error in an experiment decreases as a function of

#### Biased and Unbiased variance

Where does this idea of an unbiased variance come from?

Using the same procedure as the previous section let's look at the average variance of the variances.

A sample variance of measurements of is

To determine the "true" variance consider taking average of several sample variances (this is the same argument used above which let to )

: as shown previously
: also shown previously the universe average is the same as the sample average

Here

the sample variance
an average of all possible sample variance which should be equivalent to the "true" population variance.
: if all the variances are the same this would be equivalent
unbiased sample variance

## Probability Distributions

### Mean(Expectation value) and variance for a Given Probability Distribution

#### Mean of Discrete Probability Distribution

In the case that you know the probability distribution you can calculate the mean or expectation value E(x) and standard deviation as

For a Discrete probability distribution

where

number of observations

number of different possible observable variables

ith observable quantity

probability of observing = Probability Distribution for a discrete probability distribution

#### Mean of a continuous probability distibution

The average (mean) of a sample drawn from any probability distribution is defined in terms of the expectation value E(x) such that

The expectation value for a continuous probability distribution is calculated as

### Variance

#### Expectation of Arbitrary function

If is an arbitrary function of a variable governed by a probability distribution

then the expectation value of is

or if a continuous distribtion

### Uniform

The Uniform probability distribution function is a continuous probability function over a specified interval in which any value within the interval has the same probability of occurring.

Mathematically the uniform distribution over an interval from a to b is given by

#### Variance of Uniform PDF

Now use ROOT to generate uniform distributions. http://wiki.iac.isu.edu/index.php/TF_ErrAna_InClassLab#Day_3

### Binomial Distribution

Binomial random variable describes experiments in which the outcome has only 2 possibilities. The two possible outcomes can be labeled as "success" or "failure". The probabilities may be defined as

p
the probability of a success

and

q
the probability of a failure.

If we let represent the number of successes after repeating the experiment times

Experiments with are also known as Bernoulli trails.

Then is the Binomial random variable with parameters and .

The number of ways in which the successful outcomes can be organized in repeated trials is

where the denotes a factorial such that .

The expression is known as the binomial coefficient and is represented as

The probability of any one ordering of the success and failures is given by

This means the probability of getting exactly successes after trials is

#### Mean

It can be shown that the Expectation Value of the distribution is

:summation starts from x=1 and not x=0 now
:factor out  : replace n-1 with m everywhere and it looks like binomial distribution
:change summation index so y=x-1, now n become n-1
:
:definition of binomial expansion
:q+p =1

#### variance

Remember

To calculate the variance of the Binomial distribution I will just calculate and then subtract off .

: x=0 term is zero so no contribution

Let m=n-1 and y=x-1

### Examples

#### The number of times a coin toss is heads.

The probability of a coin landing with the head of the coin facing up is

= Uniform distribution with a=0 (tails) b=1 (heads).

Suppose you toss a coin 4 times.

Here are the possible outcomes

 order Number Trial # # of Heads 1 2 3 4 1 t t t t 0 2 h t t t 1 3 t h t t 1 4 t t h t 1 5 t t t h 1 6 h h t t 2 7 h t h t 2 8 h t t h 2 9 t h h t 2 10 t h t h 2 11 t t h h 2 12 t h h h 3 13 h t h h 3 14 h h t h 3 15 h h h t 3 16 h h h h 4

How many possible ways can you get all tails (x=4).

0! = 1

The probability of order #1 happening is

P( order #1) =

P( order #2) =

The probability of observing the coin land on heads 3 times out of 4 trials is.

#### A 6 sided die

A die is a 6 sided cube with dots on each side. Each side has a unique number of dots with at most 6 dots on any one side.

P=1/6 = probability of landing on any side of the cube.

Expectation value :

The expected (average) value for rolling a single die.

The variance:

If we roll the die 10 times what is the probability that X dice will show a 6?

A success will be that the die landed with 6 dots face up.

So the probability of this is 1/6 (p=1/6) , we toss it 10 times (n=10) so the binomial distribution function for a success/fail experiment says

So the probability the die will have 6 dots face up in 4/10 rolls is

Mean = np = Variance =

### Poisson Distribution

The Poisson distribution is an approximation to the binomial distribution in the event that the probability of a success is quite small . As the number of repeated observations (n) gets large, the binomial distribution becomes more difficult to evaluate because of the leading term

The poisson distribution overcomes this problem by defining the probability in terms of the average .

#### Poisson as approximation to Binomial

To drive home the idea that the Poisson distribution approximates a Binomial distribution at small p and large n consider the following derivation

The Binomial Probability Distriubtions is

The term

IFF we have x terms above
then
example:

This leave us with

For

#### Derivation of Poisson Distribution

The mean free path of a particle traversing a volume of material is a common problem in nuclear and particle physics. If you want to shield your apparatus or yourself from radiation you want to know how far the radiation travels through material.

The mean free path is the average distance a particle travels through a material before interacting with the material.

If we let represent the mean free path
Then the probability of having an interaction after a distance x is

as a result

= probability of getting no events after a length dx

When we consider ( we are looking for small distances such that the probability of no interactions is high)

Now we wish to find the probability of finding events over a distance given the mean free path.

This is calculated as a joint probability. If it were the case that we wanted to know the probability of only one interaction over a distance , then we would want to multiply the probability that an interaction happened after a distance by the probability that no more interactions happen by the time the particle reaches the distance .

For the case of interactions, we have a series of interactions happening over N intervals of with the probability

= probability of finding events within the length

The above expression represents the probability for a particular sequence of events in which an interaction occurs after a distance then a interaction after ,

So in essence the above expression is a "probability element" where another probability element may be

where the first interaction occurs after the distance .

So we can write a differential probability element which we need to add up as

The N! accounts for the degeneracy in which for every N! permutations there is really only one new combination. ie we are double counting when we integrate.

Using the integral formula

we end up with

#### Variance of Poisson Dist

For Homework you will show, in a manner similar to the above mean calculation, that the variance of the Poisson distribution is

### Gaussian

The Gaussian (Normal) distribution is an approximation of the Binomial distribution for the case of a large number of possible different observations. Poisson approximated the binomial distribution for the case when p<<1 ( the average number of successes is a lot smaller than the number of trials ).

The Gaussian distribution is accepted as one of the most likely distributions to describe measurements.

A Gaussian distribution which is normalized such that its integral is unity is refered to as the Normal distribution. You could mathematically construct a Gaussian distribution that is not normalized to unity (this is often done when fitting experimental data).

= probability of observing from a Gaussian parent distribution with a mean and standard deviation .

#### Half-Width (a.k.a. Full Width as Half Max)

The half width is used to describe the range of through which the distribution's amplitude decreases to half of its maximum value.

ie
Side note
the point of steepest descent is located at such that

#### Probable Error (P.E.)

The probable error is the range of in which half of the observations (values of ) are expected to fall.

#### Binomial with Large N becomes Gaussian

Consider the binomial distribution in which a fair coin is tossed a large number of times (N is very large and an EVEN number N=2n)

What is the probability you get exactly heads and tails where s is an integer?

The Binomial Probability distribution is given as

p = probability of success= 1/2

q= 1-p = 1/2

N = number of trials =2n

x= number of successes=n-s

Now let's cast this probability with respect to the probability that we get an even number of heads and tails by defining the following ratio R such that

Take the natural logarithm of both sides

Stirling's Approximation says

similarly

For Large

If

Then

or

as a result

In binomial distributions

for this problem

or

= probability of exactly heads AND tails after flipping the coin N times (N is and even number and s is an integer).

If we let and realize that for a binomial distributions

Then

So when N gets big the Gaussian distribution is a good approximation to the Binomianl

#### Gaussian approximation to Poisson when

= Poisson probability distribution

substitute

Stirling's Approximation when

: if Note:

another mathematical identity

if then

In the Poisson distribution

replacing dummy variable x with r -

= Gaussian distribution when

#### Integral Probability (Cumulative Distribution Function)

The Poisson and Binomial distributions are discrete probability distributions (integers).

The Gaussian distribution is our first continuous distribution as the variables are real numbers. It is not very meaningful to speak of the probability that the variate (x) assumes a specific value.

One could consider defining a probability element which is really an integral over a finite region such that

The advantage of this definition becomes apparent when you are interested in quantifying the probability that a measurement would fall outside a range .

The Cumulative Distribution Function (CDF), however, is defined in terms of the integral from the variates min value

Probability that you measure a value less than or equal to
##### discrete CDF example

The probability that a student fails this class is 7.3%.

What is the probability more than 5 student will fail in a class of 32 students?

Answ:

= 8%

There is an 8% probability that 5 or more students will fail the class

##### 2 SD rule of thumb for Gaussian PDF

In the above example you calculated the probability that more than 5 student will fail a class. You can extend this principle to calculate the probability of taking a measurement which exceeds the expected mean value.

One of the more common consistency checks you can make on a sample data set which you expect to be from a Gaussian distribution is to ask how many data points appear more than 2 S.D. () from the mean value.

The CDF for this is

Let

The above integral can be done numerically by expanding the exponential in a power series

The probability is one which you can compute.

Below is a table representing the cumulative probability for events to occur outside and interval of in a Gaussian distribution

 1 2 3 4

### Cauchy/Lorentzian/Breit-Wigner Distribution

In Mathematics, the Cauchy distribution is written as

= Cauchy-Lorentian Distribution
Note; The probability does not fall as rapidly to zero as the Gaussian. As a result, the Gaussian's central peak contributes more to the area than the Lorentzian's.

This distribution happens to be a solution to physics problems involving forced resonances (spring systems driven by a source, or a nuclear interaction which induces a metastable state).

= Breit-Wigner distribution
mass resonance
FWHM
= uncertainty principle

A Beit-Wigner function fit to cross section measured as a function of energy will allow one to evaluate the rate increases that are produced when the probing energy excites a resonant state that has a mass and lasts for the time derived from the Half Width .

#### mean

Mean is not defined

Mode = Median = or

#### Variance

The variance is also not defined but rather the distribution is parameterized in terms of the Half Width

Let

Then

The above integral does not converge for large z [large deviations ] . The width of the distribution is instead characterized by = FWHM

### Landau

where is any positive real number.

To simplify computation it is more convenient to use the equivalent expression

The above distribution was derived by Landau (L. Landau, "On the Energy Loss of Fast Particles by Ionization", J. Phys., vol 8 (1944), pg 201 ) to describe the energy loss by particles traveling through thin material ( materials with a thickness on the order of a few radiation lengths).

Bethe-Bloch derived an expression to determine the AVERAGE amount of energy lost by a particle traversing a given material assuming several collisions which span the physical limits of the interaction.

For the case of thin absorbers, the number of collisions is so small that the central limit theorem used to average over several collision doesn't apply and there is a finite possibility of observing large energy losses.

As a result one would expect a distribution which is Gaussian like but with a "tail" on the side of the distribution.

### Gamma

where

The distribution is used for "waiting time" models. How long do you need to wait for a rain storm, how long do you need to wait to die,...

Climatologists use this for predicting how rain fluctuates from season to season.

If integer then the above distribution is a sum of independent exponential distributions

### Exponential

The exponential distribution may be used to describe the processes that are in between Binomial and Poisson (exponential decay)

## Skewness and Kurtosis

Distributions may also be characterized by how they look in terms of Skewness and Kurtosis

### Skewness

Measures the symmetry of the distribution

Skewness =

where

The higher the number the more asymmetric (or skewed) the distribution is. The closer to zero the more symmetric.

A negative skewness indicates a tail on the left side of the distribution. Positive skewness indicates a tail on the right.

### Kurtosis

Measures the "pointyness" of the distribution

Kurtosis =

where

K=3 for Normal Distribution

In ROOT the Kurtosis entry in the statistics box is really the "excess kurtosis" which is the subtraction of the kurtosis by 3

Excess Kurtosis =

In this case a Positive excess Kurtosis will indicate a peak that is sharper than a gaussian while a negative value will indicate a peak that is flatter than a comparable Gaussian distribution.