Experiment 1

Error Analysis and Graph Drawing

I. Introduction:
I.1 It is impossible to do an experimental measurement with perfect accuracy. There is
always an uncertainty associated with any measured quantity in an experiment even in
the most carefully done experiment and despite using the most sophisticated
instruments. This uncertainty in the measured value is known as the error in that
particular measured quantity. There is no way by which one can measure a quantity
with one hundred percent accuracy. In presenting experimental results it is very
important to objectively estimate the error in the measured result. Such an exercise is
very basic to experimental science. The importance of characterizing the accuracy and
reliability of an experimental result is difficult to understate when we keep in mind
that it is experimental evidence that validate scientific theories. Likewise, reliability
and accuracy of measurements are also deeply relevant to Engineering.
The complete science of error analysis involves the theory of statistics (see Ref. 1,2)
and is too involved to present here. This short presentation is intended to introduce the
student to some basic aspects of error analysis and graph drawing, which it is
expected that the student will then put into practice when presenting his/her results of
the coming experiments.

I.2 When a measurement of a physical quantity is repeated, the results of the various
measurements will, in general, spread over a range of values. This spread in the
measured results are due to the errors in the experiment. Errors are generally
classified into two types: systematic (or determinate) errors and random (or
indeterminate) errors. A systematic error is an error, which is constant throughout a
set of readings. Systematic errors lead to a clustering of the measured values around a
value displaced from the “true” value of the quantity. Random errors on the other
hand, can be either positive or negative and lead to a dispersion of the measurements
around a mean value. For example, in a time period measurement, errors in starting
and stopping the clock will lead to random errors, while a defect in the working of the
clock will lead to systematic error. A striking example of systematic error is the
measurement of the value of the electric charge of the electron ‘e’ by Millikan by his
 Oil Drop method. Millikan underestimated the viscosity of air, leading to a lower
value for his result
e = (1.591 ± 0.002) × 10 −19 C. (1)
Compare this with a more modern and accurate value (Cohen and Taylor 1973, Ref 3)
e = (1.602 189 ± 0.000 005) × 10 −19 C. (2)
Systematic errors need to be carefully uncovered for the particular experimental set-
up and eliminated by correcting the results of the measurements.

I.3 Random errors are handled using statistical analysis. Assume that a large number (N)
of measurements are taken of a quantity Q giving values Q1,Q2, Q3,…QN. Let Q be
the mean value of these measurements
1 N
Q = ∑ Qi ,i = 1,2 ,.....N
N 1
(3)

and let ‘d’ be the deviation in the measurements

1 N 2
d= ∑
N 1
Qi −(Q )
,i = 1,2 ,.....N . (4)

The result of the measurement is quoted (assuming systematic errors have been
eliminated) as
Q = Q ± d. (5)

The error ∆Q in the quantity Q is then taken to be the deviation d. (This is called the
standard error in Q).
In a single measurement of a physical quantity, the error can be estimated as the least
count (or its fraction) of the instrument being used.
As an example, the result of a measurement of the radius of curvature R, of a plaino-
convex could be quoted as
R = 140 ± 0.2 cm. (6)
This means that we expect that the value of R to be in the range 139.8 to 140.2 cm.
Note however, that this does not mean that the “true” value of R necessarily lies in
this range, only that there is a probability that it will do so.
The error in a measurement can also be quoted as a percent error,
∆Q d
× 100 = × 100 . (7)
Q Q
For example, the percent error in R is 0.143%.
I.4 Combination of errors:
Often the value of a quantity of interest may depend on other measured quantities. For
example we could have a quantity Q which is a function F of a number of
independent (actively controlled by us) variables say x,y and z i.e.,
Q = Q(x,y,z) (8)

In general, the error in Q is related to errors in x, y, z, …, as follows (for small errors)

(∆Q)2 = (∆Qx )2 + (∆Qy )2 + (∆Qz )2 + ... (9)

where

 ∂Q   ∂Q   ∂Q 
∆Qx =  ∆x ; ∆Qy =  ∆y ; ∆Qz =  ∆z etc.
 ∂x   ∂y   ∂z 

The following table summarizes the results for combining errors for some standard
functions. Try to derive some of these results.

Sr No Function Q(x,y) Error ∆Q or Fractional Error ∆Q/Q

1 Q=x+y ∆Q = (∆x )2 + (∆y )2

2 Q=x−y ∆Q = (∆x )2 + (∆y )2
2 2 2 2 2
 ∆Q   ∆x   ∆y   ∆Q   ∆x   ∆y 
3 Q=xy   =   +   ⇒   =   +  
 Q   x   y   Q   x   y 

2 2 2 2 2
 ∆Q   ∆x   ∆y   ∆Q   ∆x   ∆y 
4 Q = x/y   =   +   ⇒   =   +  
 Q   x   y   Q   x   y 

∆Q ∆x
5 Q = xn =n
Q x
∆x
6 Q = ln x ∆Q =
x
∆Q
7 Q = ex = ∆x
Q
II Drawing of best fit straight line graph:

Below we describe how to fit a straight line to a set of data. Relations that are not
linear can be transformed to a linear one by an appropriate transformation of the
variables (as you will learn from these assignments).

To draw the best fit straight line graph through a set of scattered experimental data
points we will follow a standard statistical method, known as least squares fit method.
Let us consider a set of N experimental data points (x1,y1), (x2,y2), (x3,y3)…(xN,yN). It
is well known that a straight line graph is described by the equation
y = mx + c. (10)
We ask the question: how are the slope m and the y-intercept c to be determined such
that a straight line best approximates the curve passing through the data points? Let
Si = yi - mxi – c be the deviation of any experimental point P(xi, yi), from the best fit
line. Then, the gradient ‘m’ and the intercept ‘c’ of the best fit straight line has to be
found such that the quantity
S = ∑ ( yi − mxi − c )
2

is a minimum. We thus require

∂S ∂S
= −2∑ xi ( yi − mxi − c ) = 0 and = −2∑ ( yi − mxi − c ) = 0 ,
∂m i ∂c i

which give,
m∑ xi + c∑ xi = ∑ xi yi m∑ xi + Nc = ∑ yi .
2
and

The second equation can be rewritten as y = mx + c , where

1 1 
y=
N
∑ yi and x =  ∑ xi  showing that the best fit straight line passes
N 
through the centroid ( x , y ) of the points (xi,yi). The required values of m and c can be
calculated from the above two equations to be

m=
∑ (x − x ) y
i i
and c = y − mx . (11)
∑ (x − x )
2
i

The best fit straight line can be drawn by calculating m and c from above. A graphical
method of obtaining the best fit line is to rotate a transparent ruler about the centroid
so that it passes through the clusters of points at the top right and at the bottom left.
 This line will give the maximum error in m, (∆m)1 on one side. Do the same to find
out the maximum error in m, (∆m)2 on the other side. Now bisect the angle between
these two lines and that will be the best fit line through the experimental data.
What are the errors in the gradient and intercept due to errors in the experimental data
points? The estimates of the standard errors in the slope and intercept are

1 ∑ Si  1 x 2  ∑ Si
2 2

(∆m) 2
≈ and (∆c )
2
≈  +  ,
D N −2  N D  N −2

where D = ∑ ( xi − x ) and Si is the deviation Si = yi - mxi – c.

2

II.1 Presentation of error associated with experimental data in a graph:

Let us consider a function, y = f(x), where x is an independent parameter which in the

hand of the experimentalist during performing the experiments and y is the
experimental data which is having a value depending upon the x and the instruments.
Let the error associated with x be ±∆x and that for y be ±∆y. One can represent ±∆x
and ±∆y with the experimental data point P(x,y) on the graph paper. To do that, first
plot P(x,y) on the graph paper, then draw a vertical line parallel to y axis about the
point P(x,y) of length 2∆y. So upper half of the line represents the error +∆y and the
lower half represents –∆y error. To present ±∆x, draw horizontal lines at the two ends
of the vertical line of length 2∆x each. The whole presentation is now giving the
errors associated with the experimental point P(x,y).
Figure 1 is an example of experimental
data of resonance absorption of γ-ray
experiment (Mössbauer spectroscopy)
with error associated with each
experimental data. The solid line gives
the fitted curve through the experimental
data. Note that the error in the variable
along horizontal axis is not shown.
II.3 Use of graphs in experimental physics:
Fig.1
In practical physics, the graph of the experimental data is most important in
improving the understanding of the experimental results. Moreover from the graphs
 one can calculate unknowns related to the experiments and one can compare the
experimental data with the theoretical curve when they are presented on same graph.
There are different types of graph papers available in market. So, one should choose
the appropriate type of graph paper to present their experimental results in the best
way depending upon the values of the experimental data and the theoretical
expression of the functions. To understand all those some of the assignments are
given below in addition to those we discussed before.

III Using calculator to fit curves

After doing these assignments you also have to verify the best fit parameters, you
have obtained, by using your calculator. Calculators do the kind of curve fitting we
want. Find out the procedure, specific to your calculator, to fit curves to a set of data.
You will find it in the calculator’s manual under heading like “statistical” and sub
heading “Regression”. If the formula of interest is one of the standard relations in the
calculator then you simply enter the data and read of the best fit parameters that the
calculator responds with. For other formulae you can use the trick of transforming to
a new set of variables. In all your future work you can do the regression analysis
using a calculator.
IV. Exercises and Viva Questions
1. What is the general classification of errors? Give an example of each. How are
they taken care of? Radius of curvature
(cm)
2. What is the meaning of standard 130.121
error? Calculate the standard error 130.136
130.139
for the hypothetical data given in the 130.148
130.155
adjacent table. Express the quantity 130.162
130.169
as in eq. (5), i.e. R = R ± d .
3. What is the percent error in Millikan’s measurement of the charge of the electron:
e = (1.591 ± 0.002) x 10−19 C?
4. What is the error in the volume of a cube V= L3 if the error in L is 0.01m? If L is
measured as L = 2 ± 0.01, express the value of V in a similar manner.
5. A small steel ball-bearing rests on top of a horizontal table. The radius (R) of the
ball is measured using a micrometer screw gauge (with vernier least count 0.05
mm) to be 2.15 mm. The height of the table is found using an ordinary meter scale
 to be 90 cm. What is the height of the center of the steel ball from the floor
(include the error)?
6. Let Q = x − y, where x = 100 ± 2 and y = 96 ± 2. Calculate Q (express the result
with the error included).
7. Consider the quantity Q = x/y. If x = 50 ± 1 and y = 3 ± 0.2. Calculate Q (express
the result with the error included).
8. In an experiment involving diffraction of sodium light using a diffraction grating,
the doublet lines are unresolved at first order and a single spectral line is seen at
an angle of 13o. If the least count of the vernier of the telescope is 1′, what will be
the error in the calculated value of the grating constant d? (Principal maxima of a
grating occur at angles θ such that dsinθ = mλ. The wavelength separation
between the sodium doublet lines is 6 Ao ).
9. Consider an experiment to measure the gravitational acceleration g by measuring
the time period of a simple pendulum. What are the possible sources of systematic
error in this experiment?
10. “If there are always errors in any measurement then there is nothing like the ‘true’
value of any measured quantity”. Comment on this statement. In what sense then
do you understand the values of ‘physical constants’ to be constants?
References:
1. “Practical Physics”, G.L. Squires, Cambridge University Press, Cambridge,
1985.
2. “Laboratory Experiments in College Physics”, C.H. Bernard and C.D. Epp,
John Wiley and Sons, Inc., New York, 1995.
3. Cohen, E.R. and Taylor, B.N., J. Phys. and Chem. Reference Data, Vol 2,
page 663, 1973.
 Experiment 1: Error analysis and graph drawing
Assignments

Assignment 1 Assignment 2 Assignment 3 Assignment 4 Assignment 5 Total

(4) (4) (4) (4) (4) (20)

Assignments:
1. Experimental data (in arbitrary units) of some experiment is given below:

X -15 -8 -1 5 10 15 22 29 36 47 55 68 79 85 100 115

Y -55 -42 -23 -18 -15 -10 3 9 17 22 33 40 70 84 92 100
X 134 145 164 180 195 210 215 223 230
Y 124 130 150 175 190 200 210 225 220

(a) Assuming 10% error in Y values, plot the data on preferred graph paper
showing the errors in terms of error bars.
(b) Calculate the slope and intercept of the best fit graph. Draw the best fit graph
on the above graph.
2. The expression of refractive index of a prism is given by the following relation:
§ A+ D·
sin¨ ¸
µ =
© 2 ¹ . Assuming the error of A and D as ∆A and ∆D, express the error
§ A·
sin¨ ¸
©2¹
of µ. Here A is the angle of the prism and D is the angle of deviation.

3. The relation between two independent variables X and Y is given as the empirical
expression Y = aX + bX3. The experimental data for X and Y are given below:

X: 0.088 0.121 0.185 0.234 0.284 0.333 0.367 0.387 0.421 0.447 0.469
Y: 0.149 0.255 0.338 0.474 0.584 0.678 0.821 0.912 0.966 1.131 1.104
X: 0.506 0.524 0.538 0.559
Y: 1.354 1.398 1.402 1.523
 Rearrange the equation to plot the graph in simpler form. (Hint: Plot Y/X vs. X2.
(Why ?)). Then find out the constants a and b from the graph. Try to co-relate the
expression with some practical experiment in physics and give your comments
about the constants.
4. Expression of some function is given by, Y = aXb, where, a and b are unknown.
Use the following experimental data to find out the constants by plotting an
appropriate graph of Y vs. X. Try to co-relate the above expression with some
practical experiment in physics and give your comments about the constants.

X: 250 365 489 623 796 956 1151 1358 1511 1778
Y: 458 2189 3680 13041 39576 58031 198234 277653 511257 582314

5. The ionic conductivity (C) of an crystal is given as a function of temperature (T)

by the equation, C = Co exp(E/kT) where k is Boltzmann constant.
(T is in Kelvin. And C is in C×107 cgs unit).

T: 789 824 845 873 910 918 945 968 986 999

C×107: 1.64 3.11 6.76 9.81 23 35.9 69 124 265 432

Plot the experimental data on suitable graph paper and find out the value of Co
and E.