Chapter 13 · Making sense of numbers
How do you summarise a set of numbers?
Mean, median, mode, SD, SE, coefficient of variation, skewness and the normal curve.
27 TNMGRMU questions answered hereset 35 times, 1990–2025
The problem first
You are a second-year MS General Surgery resident. Your thesis table reports hospital stay after emergency laparotomy as 6.2 ± 1.1 days for 40 patients. The software printed the standard error, and it looked tidy, so you used it.
Two of your patients had anastomotic leaks and stayed more than five weeks. At the viva the examiner reads the table and asks one question.
“Is 1.1 the standard deviation or the standard error, and should this be a mean at all?”
Start here
What does a column of numbers tell you?
A list of 40 hospital stays is hard to read. Descriptive statistics reduce it to a few numbers: a typical value, how much the values vary, and the shape of the distribution.
The typical value is a measure of central tendency (mean, median, mode). The variation is a measure of dispersion (range, interquartile range, variance, standard deviation). The shape is described by skewness and kurtosis. Draw the histogram before you choose any of them. How to draw it is in the data presentation chapter.
Central tendency
What is the typical value: mean, median or mode?
The mean is the sum of the values divided by their number. The median is the middle value when the values are put in order. The mode is the most frequent value.
| Measure | Uses every value? | Affected by extreme values? | Best for |
|---|---|---|---|
| Mean | Yes | Yes, strongly | Symmetrical measurements (haemoglobin, height) |
| Median | No, only their order | Very little | Skewed data, scores, data with extreme values (hospital stay, income) |
| Mode | No | No | Categories (the commonest blood group). There may be none, or more than one |
A distribution with two peaks is bimodal. It often means two different groups have been mixed, such as the heights of men and women together. Describe the groups separately, because a single mean or median falls between the two peaks and describes neither group. The mode is also unstable in small samples, where one or two extra values can move it.
Interactive
Move one patient's stay and watch the mean, the median and the SD
Teaching example: six patients stay 2, 3, 3, 3, 4 and 5 days. Move the seventh patient's stay and watch both averages.
Drag the red dot or move the slider
mean 3.7 · median 3 · SD 1.4 days
With a seventh stay of 6 days: mean 3.7, median 3, SD 1.4 days. With 40 days: mean 8.6, median still 3, SD 13.9 days. One extreme value changes the mean and the SD. The median does not change, because it depends only on which value is in the middle.
Central tendency
When is the ordinary mean the wrong average?
When values rise in multiples, when you average rates, and when you combine groups of different sizes.
Geometric mean: values that rise in multiples
GM = antilog of the mean of the logs = (x₁ × x₂ × … × xₙ)1/n
Teaching example: antibody titres in seven people after vaccination, 1:20, 1:40, 1:40, 1:80, 1:160, 1:320 and 1:640. Each dilution step doubles the titre, so the values are analysed on a log scale. The mean of the log₂ values is 6.608, and 26.608 = 97.5, a geometric mean titre (GMT) of about 1:98. The arithmetic mean, 185.7, is pulled up by the two highest titres. The median is 80. The geometric mean cannot be calculated if any value is zero or negative.
Harmonic mean: averaging rates
HM = n ÷ Σ(1 ÷ x)
Teaching example: two pharmacists each dispense 60 prescriptions, one at 20 per hour and one at 30 per hour. Together they take 3 + 2 = 5 hours for 120 prescriptions, which is 24 per hour. The harmonic mean gives the same answer, 2 ÷ (1/20 + 1/30) = 24, while the ordinary mean of the two rates, 25, is wrong. The harmonic mean suits rates when the same amount of work is done at each rate.
Combined mean of groups of different sizes
Combined mean = (n₁x̄₁ + n₂x̄₂) ÷ (n₁ + n₂)
Teaching example: ward A, 40 patients, mean stay 6.0 days. Ward B, 60 patients, mean stay 4.5 days. Combined mean = (240 + 270) ÷ 100 = 5.1 days. The average of the two means, 5.25, gives the smaller ward too much weight. Medians cannot be combined in this way.
Dispersion
How much do the values vary?
Two groups can have the same mean and be very different. Measures of dispersion describe how widely the values are spread.
| Measure | Definition | Strength | Weakness |
|---|---|---|---|
| Range | Largest − smallest value | Simple | Uses two values only, and grows with sample size and extreme values |
| Interquartile range (IQR) | 75th percentile − 25th percentile | Not affected by extreme values, and pairs with the median | Ignores the outer half of the data |
| Quartile deviation | IQR ÷ 2 | As for IQR | As for IQR |
| Mean deviation | Average of the absolute distances from the mean | Uses every value | Rarely used in further calculations |
| Variance | Average squared distance from the mean (divided by n − 1 in a sample) | Basis of ANOVA and many tests | In squared units (mmHg²) |
| Standard deviation (SD) | Square root of the variance | Same units as the data, and uses every value | Affected by extreme values |
| Coefficient of variation (CV) | SD ÷ mean × 100 | No units, so it compares variation between different measurements | Meaningless when the mean is near zero |
Why square the deviations? Deviations from the mean always add up to zero, so their signs must be removed. Absolute values do this, and give the mean deviation. Squares are used for the variance because variances add: for two independent measurements, the variance of their sum or difference is the sum of their variances. If two independent readings each have an SD of 3 mmHg, the difference between them has an SD of √(3² + 3²) = 4.24 mmHg, not 6 and not 0.
Notation: μ and σ for the population mean and SD, x̄ and s for the sample mean and SD. σ² and s² are the variances.
The formula
How do you calculate the mean, SD and CV?
Teaching example: haemoglobin in six women, 10, 11, 12, 12, 13 and 14 g/dL.
mean (x̄) = Σx ÷ n
SD = √[ Σ(x − x̄)² ÷ (n − 1) ]
CV = SD ÷ x̄ × 100 · SE = SD ÷ √n
- x
- each value
- n
- number of values
- Σ
- add up
- n − 1
- used in a sample, because deviations from the sample mean are slightly smaller than deviations from the population mean
| x (g/dL) | x − x̄ | (x − x̄)² |
|---|---|---|
| 10 | −2 | 4 |
| 11 | −1 | 1 |
| 12 | 0 | 0 |
| 12 | 0 | 0 |
| 13 | 1 | 1 |
| 14 | 2 | 4 |
| Σx = 72 | 0 | Σ = 10 |
x̄ = 72 ÷ 6 = 12.0 g/dL
variance = 10 ÷ 5 = 2.0 → SD = √2.0 = 1.41 g/dL
CV = 1.41 ÷ 12.0 × 100 = 11.8% · SE = 1.41 ÷ √6 = 0.58 g/dL
The same calculation in Python:
import statistics as st
hb = [10, 11, 12, 12, 13, 14] # g/dL, six women (teaching example)
mean = st.mean(hb)
sd = st.stdev(hb) # divides by n - 1
cv = sd / mean * 100
se = sd / len(hb) ** 0.5
print(f"mean {mean:.1f}, SD {sd:.2f}, CV {cv:.1f}%, SE {se:.2f}")
# mean 12.0, SD 1.41, CV 11.8%, SE 0.58Median 12, mode 12, range 10–14, IQR 11.25–12.75 (linear interpolation), mean deviation 1.0 g/dL. Dividing by n instead of n − 1 gives 1.29, which underestimates the population SD.
What happens when every value is changed in the same way?
| Change to every value | Mean | SD | CV |
|---|---|---|---|
| Add 5 (change of origin) | 12.0 → 17.0 | 1.41 → 1.41, unchanged | 11.8% → 8.3% |
| Multiply by 10 (g/dL to g/L, change of scale) | 12.0 → 120 | 1.41 → 14.1 | 11.8%, unchanged |
Adding a constant moves the centre but not the spread. Multiplying changes both by the same factor. The CV changes with the origin, which is why it is used only for measurements with a true zero.
The formula
How do you get the mean, median and mode from a frequency table?
When only grouped data are available, each value is assumed to lie spread evenly within its class interval. The results are approximations.
Teaching example (exam-style, made-up): systolic blood pressure of 150 adults. 100–110 means 100 to under 110 mmHg.
| Class (mmHg) | Midpoint x | Frequency f | f × x | Cumulative f |
|---|---|---|---|---|
| 100–110 | 105 | 8 | 840 | 8 |
| 110–120 | 115 | 22 | 2,530 | 30 |
| 120–130 | 125 | 41 | 5,125 | 71 |
| 130–140 | 135 | 38 | 5,130 | 109 |
| 140–150 | 145 | 24 | 3,480 | 133 |
| 150–160 | 155 | 12 | 1,860 | 145 |
| 160–170 | 165 | 5 | 825 | 150 |
| Total | 150 | 19,790 |
Mean = Σfx ÷ n = 19,790 ÷ 150 = 131.9 mmHg
Median = L + [(n/2 − cf) ÷ f] × h = 130 + [(75 − 71) ÷ 38] × 10 = 131.1 mmHg
Mode = L + [(f₁ − f₀) ÷ (2f₁ − f₀ − f₂)] × h
Mode = 120 + [(41 − 22) ÷ (82 − 22 − 38)] × 10 = 128.6 mmHg
- L
- lower limit of the class that contains the median (or the mode)
- cf
- cumulative frequency of the classes before it
- f
- frequency of that class. For the mode, f₁ is the largest frequency, f₀ and f₂ the frequencies of the classes before and after
- h
- width of the class, here 10 mmHg
Quartiles use the same formula as the median, with n/4 and 3n/4 in place of n/2. Q1 = 120 + [(37.5 − 30) ÷ 41] × 10 = 121.8, Q3 = 140 + [(112.5 − 109) ÷ 24] × 10 = 141.5 mmHg. The SD from the midpoints is 14.3 mmHg. With the individual values, always calculate from them instead.
Comparison
SD or SE: which one goes after the ± sign?
The SD describes how much patients differ from each other. The standard error (SE) describes how much the mean would differ from study to study if you repeated the study many times.
| Feature | Standard deviation (SD) | Standard error of the mean (SE, SEM) |
|---|---|---|
| Describes | Variation between individual patients | Precision of the sample mean as an estimate of the population mean |
| Formula | √[Σ(x − x̄)² ÷ (n − 1)] | SD ÷ √n |
| As n increases | Stays about the same | Becomes smaller |
| Type of statistic | Descriptive | Inferential (the basis of confidence intervals and t-tests) |
| Use it for | Table 1, baseline characteristics, normal ranges | Confidence intervals. A 95% CI is clearer than the SE itself |
| Haemoglobin example | SD 1.41 g/dL | SE 0.58 g/dL |
Back to the thesis table: an SE of 1.1 with 40 patients means SD = 1.1 × √40 = 7.0 days. Mean − 2 SD is 6.2 − 13.9 = −7.7 days, which is impossible for a hospital stay. The data are right-skewed, so the table should give the median and IQR.
Interactive
Change the number of patients and watch the SD and the SE
The SD describes how much patients differ from each other. The SE describes how much the mean differs between samples of the same size, and it becomes smaller as the sample size increases.
Try this: move the number of patients from 25 to 100 and compare the two red brackets. The SE becomes about half as large. The SD of the patients stays close to 15 mmHg.
Teaching simulation: systolic blood pressure in a population with mean 120 and SD 15 mmHg. The numbers change a little with each new draw.
More patients: the SE becomes smaller, the SD does not
SE = 15.39 ÷ √25 = 3.08 · SD of the 40 sample means = 2.97 mmHg
Laboratory imprecision
How do you compare spread when the means differ?
An SD is in the units of the measurement and depends on the size of the values. Divide it by the mean and it has no units.
CV (%) = SD ÷ mean × 100
In a laboratory, imprecision is measured by analysing the same control material repeatedly and calculating the SD and CV. The CV lets you compare analytes with different units, one analyte at different concentrations, and one method between laboratories. Laboratories then compare the CV with an analytical performance specification. A 2015 European Federation of Clinical Chemistry and Laboratory Medicine (EFLM) consensus described three models for setting these specifications.
Use the CV only for measurements with a true zero (concentrations, weights). It has no meaning for temperature in °C, and it is unstable when the mean is close to zero.
The normal curve
What is a normal distribution?
A symmetrical, bell-shaped distribution fixed completely by its mean and its SD. Mean, median and mode are equal, and the tails extend without ever reaching zero.
z = (x − mean) ÷ SD
The standard normal distribution has mean 0 and SD 1. The z-score tells you how many SDs a value lies from the mean, so values measured in different units can be placed on one scale. It is called the Gaussian distribution after Carl Friedrich Gauss. The 95% rule is used for reference ranges (mean ± 1.96 SD) only when the data are close to normal in shape.
For a distribution of any shape, a weaker statement holds: at least 1 − 1/k² of the values lie within k SDs of the mean. That is at least 75% within 2 SD and 88.9% within 3 SD. This is Chebyshev's inequality, also called the Bienaymé–Chebyshev rule. For a normal distribution the figures are 95.4% and 99.7%.
Percentiles
What does '3rd percentile' mean on a growth chart?
A percentile is the value below which a given percentage of values lie. The 3rd percentile has 3% of values below it.
Percentiles do not need a normal distribution. They are read from the ordered data, which is why growth charts, reference ranges for skewed tests and the IQR use them.
Percentiles belong to a family of quantiles, values that cut ordered data into groups of equal size: 2 tertiles make 3 groups, 3 quartiles make 4, 4 quintiles make 5, 9 deciles make 10 and 99 percentiles make 100. A percentile is not a percentage. A child on the 90th percentile of weight is heavier than 90% of children of the same age and sex on that chart, whatever the weight.
A median is a statement about a group. A median survival of 3 years means that half the patients lived 3 years or less and half lived longer. It does not predict 3 years for any one patient.
Interactive
Move a value along the normal curve and watch its z-score and percentile
On a normal curve each z-score corresponds to one percentile. The shaded area is the share of values below the value, above it, or no further from the mean than it.
Try this: choose "Within" and set the value to 130, then 140, then 150 mmHg. The shaded shares are 68.3%, 95.4% and 99.7%, the shares within 1, 2 and 3 SD of the mean.
Teaching example: systolic blood pressure with a normal shape, mean 120 and SD 10 mmHg.
z = (value − 120) ÷ 10
z = (130 − 120) ÷ 10 = +1.00 · percentile 84.1
Shape
What are skewness and kurtosis?
Skewness describes whether one tail of the distribution is longer than the other. Kurtosis describes the peak and tails compared with a normal curve of the same SD.
Pearson's skewness = (mean − mode) ÷ SD or 3 × (mean − median) ÷ SD
For the hospital-stay curve at the top of this chapter (SD 3.02 days): (5.67 − 3.89) ÷ 3.02 = 0.59, and 3 × (5.67 − 5.00) ÷ 3.02 = 0.66. Positive values mean right skew. Zero means symmetry.
Quartile coefficient of skewness = (Q3 + Q1 − 2 × median) ÷ (Q3 − Q1)
It uses only the quartiles, so extreme values do not affect it. For the same curve, Q1 = 3.57, median 5.00 and Q3 = 7.01 days: (7.01 + 3.57 − 10.00) ÷ (7.01 − 3.57) = 0.17. The value lies between −1 and +1. Software usually gives a moment coefficient instead, based on the average cube of the deviations (m₃ ÷ m₂1.5). The different coefficients are not interchangeable, so report which one you used.
The order mode < median < mean in a right-skewed distribution is a rule of thumb, not a law. It fails most often in discrete distributions, and can fail in distributions with more than one peak or with one long tail and one heavy tail (von Hippel, 2005).
β₂ = m₄ ÷ m₂² · excess kurtosis = β₂ − 3
m₂ and m₄ are the average second and fourth powers of the deviations from the mean. m₂ is the variance. Some programs report β₂ and others report excess kurtosis, β₂ − 3, so check which one your software gives before you interpret it.
Interactive
Lengthen the right tail and watch the mean, median and mode separate
In a right-skewed distribution the mode is smaller than the median and the median is smaller than the mean. The longer the tail, the larger the gaps between them.
Try this: move the slider to the right and watch the mean move away from the median. When the longest 5% of stays begin at 11.4 days, mean − 2 SD falls below zero.
Teaching example: hospital stay with a median of 5 days and a long right tail. The slider lengthens the tail.
This is the 95th percentile. Longer tail: the mean moves further from the median
mode 4.6 · median 5.0 · mean 5.2 · SD 1.61 days · 3 × (mean − median) ÷ SD = 0.43
Mean − 2 SD = 2.0 days, above zero, so this check alone does not show the skew. The mean is still larger than the median. Median (IQR): 5.0 (4.1 to 6.1) days.
Shape
How do you check whether your data are close to normal?
Look at the data first. Use a formal test as a supplement to the plots, not a replacement for them.
- Draw a histogram or a dot plot. Look for one peak and rough symmetry.
- Compare the mean and the median. For a value that cannot be negative, check whether mean − 2 SD falls below zero.
- Draw a box plot. The median should lie near the middle of the box, with whiskers of similar length.
- Draw a normal Q–Q plot. The points should lie close to the straight line.
- If needed, add a formal test such as the Shapiro–Wilk test. With small samples such tests often fail to detect a real departure, and with large samples they flag departures too small to matter.
Right-skewed data can often be analysed after a log transformation, with results summarised as geometric means, or with non-parametric methods (see choosing a test). Q–Q plots, residual checks and transformations are covered in detail in Checking assumptions.
Reference intervals
What is a normal range, and how is it set?
A reference interval is the range that contains the central 95% of values in healthy reference people. By definition, 5% of healthy people have a value outside it.
Two ways to calculate it. If the values, or their logarithms, are close to normal, use mean ± 1.96 SD. Otherwise use the 2.5th and 97.5th percentiles of the ordered values. The Clinical and Laboratory Standards Institute (CLSI) C28-A3 guideline recommends the non-parametric (percentile) method, and at least 120 reference individuals so that confidence intervals for the two limits can be calculated (Ozarda, 2016).
| Feature | Reference interval | Clinical decision limit |
|---|---|---|
| Based on | Values in healthy reference people | Outcomes, risk or the benefit of treatment in patients |
| Answers | Is this value unusual for a healthy person? | Should this patient be diagnosed or treated? |
| Set by | The central 95% of the reference distribution | Studies relating the value to disease or complications |
| Healthy people outside | 5%, by definition | Depends on where the limit is placed |
A result outside the reference interval is not proof of disease, and a result inside it is not proof of health. The more tests are done, the more healthy people have at least one result outside:
P(at least one outside) = 1 − 0.95k
With k independent tests: 5 tests 22.6%, 10 tests 40.1%, 20 tests 64.2%.
This assumes the tests are independent. Related tests, such as liver enzymes, are not, so the true figure differs. Establishing reference intervals in a laboratory, and comparing a result with the same patient's earlier results, are covered in Laboratory quality control and reference intervals. Choosing a diagnostic cut-off is covered in the diagnostic tests chapter.
Interactive
Lengthen the tail and add tests, and watch how many healthy people fall outside the normal range
Try this: move the right tail from its shortest to its longest setting and compare the share of healthy people above the upper limit for mean ± 1.96 SD and for the percentiles. Then raise the number of tests to 10.
Teaching example (made-up): a blood test in healthy adults, median 100 units, with a right-skewed distribution.
Longer tail → more skew
Each test flags 5% of healthy people
| Method | Interval | Healthy below | Healthy above |
|---|---|---|---|
| Mean ± 1.96 SD | 20 to 197 | under 0.1% | 4.5% |
| 2.5th to 97.5th percentile | 46 to 219 | 2.5% | 2.5% |
Mean 108.3, SD 45.1. Mean ± 1.96 SD puts under 0.1% of healthy people below and 4.5% above, instead of 2.5% on each side. The percentile interval always leaves 2.5% on each side.
Distributions
Which named distributions will you meet?
A distribution describes how often each value occurs. A few named distributions appear again and again in medical statistics.
| Distribution | Data | Fixed by | Where you meet it |
|---|---|---|---|
| Normal | Continuous, symmetrical | Mean and SD | Many measurements in healthy adults. Means of large samples. z-scores and confidence intervals |
| Log-normal | Continuous, positive, right-skewed | Mean and SD of the logs | Hospital stay, antibody titres, many laboratory values |
| t | Continuous, symmetrical, heavier tails than normal | Degrees of freedom | t-tests and confidence intervals for means of small samples |
| Chi-square | Continuous, zero upwards | Degrees of freedom | Chi-square tests for tables of counts |
| F | Continuous, zero upwards | Two degrees of freedom | ANOVA and comparison of two variances |
| Binomial | Whole-number counts out of a fixed n | n and the probability p | Number cured out of 50 treated patients |
| Poisson | Whole-number counts with no fixed upper limit | Mean number per unit of time or space | New cases per week, colonies per culture plate |
Binomial and Poisson distributions are discrete: only whole numbers are possible. The others are continuous. The t, chi-square and F distributions are used mainly to calculate p-values, in the hypothesis testing chapter. Counts and rates with a Poisson distribution are modelled in Poisson and log-linear models.
Categorical data
How do you summarise categories?
Give the number and the percentage in each category, and say what the percentage is out of.
| Measure | What it is | Teaching example |
|---|---|---|
| Proportion (percentage) | Part ÷ whole. The numerator is included in the denominator | 24 anaemic of 120 girls = 0.20, written 20% (24 of 120) |
| Ratio | One quantity ÷ another, separate quantity | 520 women and 480 men: 520 ÷ 480 × 1,000 = 1,083 women per 1,000 men |
| Rate | Events in a population over a stated period of time | 12 new cases per 1,000 people per year |
Percentage points or per cent? A rise from 20% to 25% is an increase of 5 percentage points, and a relative increase of 25%. Write which one you mean. Averages for categories: the mode is the only average for nominal data (blood group). For ordered categories (cancer stage) the median category can also be given.
Incidence, prevalence and the rates used in epidemiology are covered in the risk and survival chapter. Tables of questions that allow more than one answer, where the percentages add up to more than 100, are covered in the data presentation chapter.
Comparison
Are you describing your patients, or drawing conclusions about all patients?
Descriptive statistics summarise the data you collected. Inferential statistics use that sample to draw conclusions about the population, with a stated uncertainty.
| Feature | Descriptive statistics | Inferential statistics |
|---|---|---|
| Question | What do my data look like? | What can I conclude about the population? |
| Tools | Mean, median, mode, SD, IQR, percentages, tables, graphs | Confidence intervals, p-values, t-test, chi-square test, ANOVA, regression |
| Haemoglobin example | Mean 12.0 g/dL, SD 1.41 g/dL in these six women | 95% CI for the population mean: 10.5 to 13.5 g/dL |
| Uncertainty stated? | No | Yes |
| In a paper | Table 1, baseline characteristics | Results of the comparisons |
95% CI = 12.0 ± 2.571 × 0.58 (t with 5 degrees of freedom). Tests and confidence intervals are covered in the hypothesis testing chapter.
Decide
Which summary should you report?
Answer two or three questions about each variable.
What do you want to describe?
Which summary: all possible answers
- Numbers and percentages: Give the number and percentage in each category. The most frequent category is the mode.
- Median and interquartile range: Write median (IQR), for example 5 (3–8) days. Add the range if the extremes matter.
- Mean and standard deviation: Write mean (SD), for example 12.0 (1.41) g/dL, not mean ± SD.
- 95% confidence interval (or standard error): SE = SD ÷ √n. A 95% confidence interval is the clearer way to show precision.
- Coefficient of variation: CV = SD ÷ mean × 100. It has no units, so methods with different units or means can be compared.
By specialty
Descriptive statistics problems by specialty
Internal quality control: the same control serum is analysed every day for a month for glucose and for creatinine. The quality manager must say which method is more precise.
Compare the CVs. Each SD is in the analyte's own units and depends on the size of the values.
A diabetes clinic audit summarises HbA1c, whose histogram is roughly symmetrical, and serum triglycerides, whose histogram has a long right tail.
Mean (SD) for HbA1c, median (IQR) for triglycerides. Decide from the histogram, one variable at a time.
Hospital stay after laparotomy: most patients go home within a week, and a few with complications stay for several weeks.
Report median (IQR). A few long stays raise the mean and the SD but not the median.
A paper shows postoperative pain scores with small error bars labelled SEM. A review of four anaesthesia journals from 2001 found SEM used to describe sample variability in 198 of 860 articles (23%).
Variability between patients is shown with the SD. The SEM shows only how precisely the mean is estimated.
A newborn's weight plots below the 3rd percentile on a birth-weight chart.
97% of babies on that reference chart weigh more. On a normal-shaped chart the 3rd percentile is about 1.88 SD below the mean.
Gestational age at delivery in a labour ward: most births are at term, and a smaller number of preterm births are spread over many weeks before term.
This is a left (negative) skew. The mean is lower than the median. Report the median (IQR) and the percentage born preterm.
A household survey records monthly income. A few households earn many times more than the rest.
Income is right-skewed. The median describes the typical household better than the mean.
Plasma drug levels at 2 hours in 20 volunteers. Two slow metabolisers have levels several times higher than the others.
Plot the values first. Report median and range, or summarise on the log scale if the log values are symmetrical.
Common myths
Common mistaken beliefs
Myth: The SEM is fine for error bars and makes the data look tidy.
The SEM describes the precision of the mean, not the spread of patients, and it is always smaller than the SD. A review of four anaesthesia journals from 2001 found this misuse in 198 of 860 articles (23%) (Nagele, 2003).
Myth: The SD can be used only for normally distributed data.
The SD is a valid measure of variability whatever the distribution (Altman and Bland, 2005). The statement "about 95% of values lie within 2 SD" is true only for data close to normal in shape.
Myth: The mean is always the best average.
With skewed data or extreme values, use the median. One stay of 40 days raises the mean of seven stays from 3.7 to 8.6 days. The median stays at 3.
Myth: A "normal distribution" means the values of healthy people.
It is the name of one bell-shaped distribution. Values from healthy people can be skewed, and values from patients can be normal in shape.
Myth: Enrolling more patients makes the SD smaller.
The SD estimates the spread in the population, so it stays about the same. The SE becomes smaller (SE = SD ÷ √n).
Myth: A smaller SD always means a more precise test.
Only when the means are similar. Glucose SD 12 at mean 101 (CV 11.9%) is more precise than creatinine SD 0.5 at mean 1.2 (CV 41.7%).
Why this matters to you
Five consequences of getting descriptive statistics wrong
A protocol that does not match the data
The analysis plan says "all data will be expressed as mean ± SD", and then hospital stay, drain output and costs turn out to be skewed.
Thesis figures that hide the variation
SEM error bars look narrow and suggest that patients responded alike when they did not.
A common viva question
Examiners ask the difference between SD and SE, and why you chose the mean or the median for each variable.
Misreading papers
A stay reported as 4 (5) days with mean − 2 SD below zero is skewed. The mean overstates the typical stay.
Laboratory results that change with repetition
With creatinine SD 0.5 mg/dL at a mean of 1.2, about half of repeat pairs on the same sample differ by more than 0.5 mg/dL, enough to move a patient across a clinical cut-off.
History
History of descriptive statistics
- 1733
Abraham de Moivre
In a Latin pamphlet, gave an approximation to the binomial distribution for a large number of trials. The expression is the one now recognised as the normal curve. Historians disagree about whether this counts as discovering the normal distribution.
- 1809
Carl Friedrich Gauss
In Theoria motus corporum coelestium, a book on calculating planetary orbits, presented the method of least squares with the normal curve as the distribution of errors of observation. The normal distribution is also called the Gaussian distribution.
- 1893–94
Karl Pearson
Introduced the term "standard deviation" in 1893 and used it in print in 1894, where he equated it with Gauss's "mean error". In the same 1894 paper he called the error curve the "normal curve".
- 1895
Karl Pearson
In a paper on skew variation, used the term "mode" for the value of maximum frequency and measured skewness from the distance between the mean and the mode.
- 1896
Karl Pearson
Defined the coefficient of variation: the standard deviation divided by the mean, multiplied by 100.
- 1905
Karl Pearson
In Biometrika, called β₂ − 3 "the degree of kurtosis" and used the terms platykurtic, leptokurtic and mesokurtic.
- 1918
Ronald A. Fisher
Named the square of the standard deviation the "variance", in a paper on inheritance.
- 2003
Peter Nagele
Found the SEM used to describe sample variability in 198 of 860 articles (23%) in four anaesthesia journals published in 2001.
- 2005
Douglas Altman and Martin Bland
BMJ statistics note "Standard deviations and standard errors": the two are often confused, and the SD is a valid measure of variability regardless of the distribution.
Guidelines
Descriptive statistics in reporting guidelines
“Summarize data that are approximately normally distributed with means and standard deviations (SD). Use the form: mean (SD), not mean ± SD.” “Summarize data that are not normally distributed with medians and interpercentile ranges, ranges, or both.” “Do NOT use the standard error of the mean (SE) to indicate the variability of a data set.”
“We advise authors to summarize continuous variables for each study group by giving the mean and standard deviation, or when the data have an asymmetrical distribution, as is often the case, the median and percentile range (e.g., 25th and 75th percentiles).”
Continuous variables can be summarised by the mean and SD, or by the median and percentiles when the distribution is asymmetrical. “Standard errors and CIs are not appropriate for describing variability—they are inferential rather than descriptive statistics.”
Reading a paper
Five checks on any paper's descriptive statistics
- Does the paper say whether each number after ± or in brackets is an SD, an SE or a confidence interval?
- Are skewed variables (hospital stay, costs, many laboratory values) given as median and IQR?
- For a value that cannot be negative, is mean − 2 SD below zero? If so, the data are skewed.
- Are error bars in the figures labelled, and do they show SD when they describe patients?
- Are the summaries given for each group, with the number of patients in each?
Why is it called that?
Word surgery
“Normal”
Healthy, usual, nothing wrong.
The name of one bell-shaped, symmetrical distribution. Values from healthy people need not follow it.
“Mean”
Unkind or stingy.
The arithmetic average: the sum of the values divided by their number.
“Error”
A mistake.
The standard error is not a mistake. It is the SD of an estimate, such as a sample mean, across repeated samples.
SD describes the patients. SE describes the mean. SE = SD ÷ √n, so the SE is smaller than the SD whenever n is more than 1.
Mean
Late Latin medianus: "of the middle", through Old French meien
The arithmetic mean: the sum of the values divided by their number.
Median
Latin medianus: "of the middle"
The middle value of the ordered data. Used in this statistical sense in English from 1883.
Mode
Latin modus: "measure, manner, way"
The most frequent value. Karl Pearson used the term in 1895.
Deviation
Latin de via: "off the way"
The distance of a value from the mean. "Standard deviation" dates from 1894.
Variance
Latin variare: "to change"
The square of the SD, named by R.A. Fisher in 1918.
Normal
Latin norma: "carpenter's square", a rule or pattern
In statistics, the name of a distribution, not a judgement that the values are healthy.
Skew
Old North French eskiuer: "to shy away from, avoid"
A distribution with one tail longer than the other.
Kurtosis
Greek kyrtos: "curved, convex", with leptos "slender", platys "broad, flat", mesos "middle"
Leptokurtic: sharper peak. Platykurtic: flatter top. Mesokurtic: like the normal curve.
Statistics
German Statistik, from Latin status: "condition, state"
First meant data about the condition of a state. The meaning "numerical data of any sort" dates from 1829.
Exam answers
Model answer outlines
Measures of central tendency: mean, median and mode, their definitions and uses
Ph.D. RM 2017, 2019, 2022, 2023 · MD Community Medicine 2017Definition: a single value that represents the centre of a set of data.
| Mean | Median | Mode | |
|---|---|---|---|
| Definition | Σx ÷ n | Middle value of ordered data, or the mean of the two middle values if n is even | Most frequent value |
| Example (10, 11, 12, 12, 13, 14 g/dL) | 12.0 | 12 | 12 |
| Merits | Uses every value. Basis of SD, t-test and ANOVA | Not affected by extreme values. Usable for ordinal data | Simple. The only average for nominal data |
| Demerits | Affected by extreme values. Misleading in skewed data | Ignores the actual size of most values | May be absent or multiple. Unstable in small samples |
| Use | Symmetrical measurements (Hb, BP) | Skewed data (hospital stay, income), scores | Categories (commonest blood group) |
- Relation in a symmetrical distribution: mean = median = mode. Right skew: mode < median < mean. Left skew: the reverse. This is a rule of thumb that can fail, especially in discrete data.
- For moderately skewed data, mode ≈ 3 median − 2 mean (an approximation).
- Other means: geometric mean for ratios and data analysed on the log scale, harmonic mean for rates. Titres 1:20 to 1:640: GM 97.5 against an arithmetic mean of 185.7.
- Combined mean of groups = Σ n x̄ ÷ Σ n: 40 patients at 6.0 days and 60 at 4.5 days give 5.1 days, not 5.25.
- Grouped data: mean = Σfx ÷ n. Median = L + [(n/2 − cf) ÷ f] × h. Mode = L + [(f₁ − f₀) ÷ (2f₁ − f₀ − f₂)] × h.
- With the SD (Ph.D. RM 2019): the mean gives the centre and the SD the spread. Report mean (SD) for symmetrical data.
Measures of dispersion (with central tendency)
MD Biochemistry 2023 · Ph.D. RM 2021 · MD Community Medicine 2018 · DM Pulmonary Medicine 2025Definition: numbers that describe how widely values are spread around the centre. Two groups with the same mean can differ completely in spread.
- Range: maximum − minimum. Simple, but depends on two extreme values.
- Interquartile range: Q3 − Q1, the middle 50%. Pairs with the median and is resistant to extreme values. Quartile deviation = IQR ÷ 2.
- Mean deviation: Σ|x − x̄| ÷ n. In the Hb example, 6 ÷ 6 = 1.0 g/dL.
- Variance: Σ(x − x̄)² ÷ (n − 1). Hb example 2.0 (g/dL)².
- Standard deviation: √variance. Hb example 1.41 g/dL. The most used measure, in the same units as the data.
- Coefficient of variation: SD ÷ mean × 100. Hb example 11.8%. No units.
- Percentiles: describe position and spread without assuming a shape.
- Chebyshev's inequality: whatever the shape, at least 75% of values lie within 2 SD of the mean, and at least 88.9% within 3 SD.
Choice: symmetrical data, mean (SD). Skewed data, median (IQR). Comparing variation between measurements, CV. For central tendency, see the frame above.
Standard deviation: definition, calculation and significance
DM Medical Oncology 2000, 2023 · MD Biochemistry 2013, 2015 · MS OBSTETRICS AND GYNAECOLOGY 2015 · MD Community Medicine 1990, 2013 · MD Pharmacology 2011Definition: the square root of the mean of the squared deviations from the mean (divisor n − 1 in a sample). It measures the average distance of values from the mean, in the units of the data.
- Formula: SD = √[Σ(x − x̄)² ÷ (n − 1)].
- Steps: find the mean, subtract it from each value, square, add, divide by n − 1, take the square root.
- Worked example: 10, 11, 12, 12, 13, 14 g/dL → mean 12.0, Σ(x − x̄)² = 10, variance 2.0, SD 1.41 g/dL.
- Properties: never negative, and zero only if all values are equal. Unchanged when a constant is added to every value. Multiplied by k when every value is multiplied by k. Affected by extreme values.
Significance and uses:
- Summarises variation between individuals. Report as mean (SD).
- In a normal distribution, 68.3% of values lie within ±1 SD, 95.4% within ±2 SD, 99.7% within ±3 SD. Reference ranges use mean ± 1.96 SD when data are near normal.
- Gives z-scores, z = (x − mean) ÷ SD, to compare values in different units.
- Needed for the SE (SD ÷ √n), confidence intervals, t-tests and sample size calculations.
- Gives the CV (SD ÷ mean × 100) and control limits in laboratory quality control.
- Reference intervals: mean ± 1.96 SD only for near-normal data. Otherwise the 2.5th to 97.5th percentiles, with at least 120 reference individuals (CLSI C28-A3).
Standard deviation and standard error of the mean
MD Community Medicine 2025 · Ph.D. RM 2017, 2024, 2025| Standard deviation | Standard error of the mean | |
|---|---|---|
| Measures | Variation of individual values around the sample mean | Variation of sample means around the population mean |
| Formula | √[Σ(x − x̄)² ÷ (n − 1)] | SD ÷ √n |
| Effect of larger n | Stays about the same | Becomes smaller |
| Type | Descriptive | Inferential |
| Use | Describe the sample as mean (SD) | 95% CI = mean ± 1.96 SE (large samples), t-tests |
| Example | Hb SD 1.41 g/dL (n = 6) | SE 1.41 ÷ √6 = 0.58 g/dL |
- Common error: reporting mean ± SEM to describe patients. SEM is always smaller than SD and hides the variation (found in 23% of 860 anaesthesia articles, Nagele 2003).
- SAMPL, STROBE and CONSORT 2025 advise SD (or median and percentiles) for describing variability, and not the SE.
Coefficient of variation: why imprecision is expressed as CV rather than SD (glucose vs creatinine), and mean, SD and CV in medical practice
MD Biochemistry 2011, 2014, 2023 · Ph.D. RM 2018Definition: CV = SD ÷ mean × 100. A relative measure of dispersion, without units.
Glucose vs creatinine:
- Glucose: CV = 12 ÷ 101 × 100 = 11.9%.
- Creatinine: CV = 0.5 ÷ 1.2 × 100 = 41.7%.
- The raw SDs (12 vs 0.5) cannot be compared: they are in the units of each analyte and depend on the size of the mean.
- Relative to its mean, creatinine varies about 3.5 times as much. ±2 SD spans 77–125 mg/dL for glucose but 0.2–2.2 mg/dL for creatinine.
- Conclusion: the quality manager is right. Glucose analytical performance is better.
Why CV is better than SD for imprecision: it has no units. It allows comparison between analytes, between concentration levels of one analyte, between methods and between laboratories, and it can be compared with performance specifications in quality control.
Mean, SD and CV in medical practice: mean (SD) summarises a symmetrical measurement such as Hb (12.0 (1.41) g/dL). Mean ± 1.96 SD gives a reference range for near-normal data. CV (11.8% in the Hb example) compares variability of different measurements, and expresses method imprecision in the laboratory.
Limits: meaningless when the mean is near zero or the scale has no true zero.
Correlation coefficient and coefficient of variation
MD Biochemistry 2009| Correlation coefficient (r) | Coefficient of variation | |
|---|---|---|
| Describes | Strength and direction of the linear relation between two variables | Spread of one variable relative to its mean |
| Formula | r = Σ(x − x̄)(y − ȳ) ÷ √[Σ(x − x̄)² × Σ(y − ȳ)²] | SD ÷ mean × 100 |
| Range | −1 to +1, where 0 = no linear relation | 0 upwards, in % |
| Units | None | None |
| Example | Age and systolic BP | Imprecision of a glucose method |
| Caution | Does not show cause, or agreement between two methods | Unstable when the mean is near zero |
More on r in the correlation and regression chapter.
Normal (Gaussian) distribution, skewed distribution, skewness and kurtosis
MD Community Medicine 2019 · DM Neonatology 2017 · Ph.D. RM 2017 · MD Pharmacology 2011Normal distribution, definition: a continuous, symmetrical, bell-shaped probability distribution fixed by its mean (μ) and SD (σ).
- Mean = median = mode. Symmetrical about the mean. Total area = 1. Tails never touch the axis.
- 68.3% within μ ± σ, 95.4% within μ ± 2σ, 99.7% within μ ± 3σ; 95% within μ ± 1.96σ.
- Standard normal: z = (x − μ) ÷ σ, mean 0, SD 1. Skewness 0, kurtosis 3.
- Examples: many biological measurements in healthy adults are close to normal, such as height. Check with a histogram or normal Q–Q plot.
Skewed distribution: one tail is longer.
- Positive (right) skew: mode < median < mean. Example: hospital stay, household income.
- Negative (left) skew: mean < median < mode. Example: gestational age at delivery.
- Pearson's coefficient = (mean − mode) ÷ SD, or 3(mean − median) ÷ SD. Example: (5.67 − 3.89) ÷ 3.02 = 0.59.
Kurtosis: β₂ = m₄ ÷ m₂². Normal = 3 (mesokurtic). Above 3, leptokurtic (sharper peak, heavier tails). Below 3, platykurtic (flatter top). Excess kurtosis = β₂ − 3.
Checking normality: histogram, mean compared with median, box plot, normal Q–Q plot, and if needed the Shapiro–Wilk test. Other distributions: log-normal, t, chi-square and F (continuous), binomial and Poisson (discrete).
Why it matters: skewed data are summarised with median (IQR) and analysed with non-parametric tests or after log transformation. See choosing a test.
Concept of standard deviation and percentiles
DM Neonatology 2017- SD: average distance of values from the mean, √[Σ(x − x̄)² ÷ (n − 1)]. See the SD frame above.
- Percentile: the value below which a stated percentage of observations fall. P50 = median. P25 and P75 = quartiles. IQR = P75 − P25.
- Link between them: in a normal distribution each z-score has a fixed percentile. z −2 = 2.3rd, z −1.88 = 3rd, z −1 = 15.9th, z 0 = 50th, z +1 = 84.1st, z +1.88 = 97th, z +2 = 97.7th.
- Neonatal use: birth weight, length and head circumference are plotted on percentile charts. A value below the 3rd or above the 97th percentile, or beyond ±2 SD, is flagged for review.
- Difference: SD-based cut-offs assume a near-normal shape. Percentiles are read from the ordered data and work for skewed measurements too.
- Other quantiles: tertiles, quartiles, quintiles and deciles. Reference intervals use the 2.5th and 97.5th percentiles.
Calculation of standard deviation and t-value of small experimental groups
MD Biochemistry 1994Teaching example: haemoglobin (g/dL) in a control group A and a treated group B, six in each.
- A: 10, 11, 12, 12, 13, 14 → mean 12.0, Σ(x − x̄)² = 10, variance 2.0, SD 1.41.
- B: 12, 13, 13, 14, 15, 17 → mean 14.0, Σ(x − x̄)² = 16, variance 3.2, SD 1.79.
- Pooled variance = (10 + 16) ÷ (6 + 6 − 2) = 2.6, and pooled SD = 1.61.
- SE of the difference = 1.61 × √(1/6 + 1/6) = 0.93.
- t = (14.0 − 12.0) ÷ 0.93 = 2.15, degrees of freedom = 10.
- Critical t (two-sided, 5%, 10 df) = 2.228. Since 2.15 < 2.228, p = 0.057: not significant at 5%.
- Small groups have wide uncertainty, so a 2 g/dL difference was not detected. Assumptions: independent groups, near-normal data, similar variances.
Principle and types of ANOVA (analysis of variance)
MD Biochemistry 2025Principle: ANOVA compares the means of three or more groups by dividing the total variation into variation between groups and variation within groups. If the group means really differ, the between-group variation is large compared with the within-group variation.
F = mean square between ÷ mean square within
Teaching example, three groups of three: 4, 5 and 6 (mean 5), then 6, 7 and 8 (mean 7), then 9, 10 and 11 (mean 10). Sum of squares between = 38 (df 2), within = 6 (df 6). Mean squares 19 and 1. F = 19.0, p = 0.003.
Types:
- One-way ANOVA: one factor (three drug doses).
- Two-way ANOVA: two factors, with or without interaction (drug and sex).
- Repeated-measures ANOVA: the same patients measured at several times.
- Factorial and multi-way ANOVA: several factors together.
- ANCOVA: ANOVA adjusted for a continuous covariate (baseline value).
- MANOVA: more than one outcome variable.
Assumptions: independent observations, near-normal residuals, similar variances. A significant F says at least one mean differs. Post hoc tests (Tukey, Bonferroni) show which. Non-parametric alternative: Kruskal–Wallis test. See choosing a test.
The one thing to remember
Draw the histogram first: report mean (SD) for symmetrical data, median (IQR) for skewed data, and never use the SE to describe patients.
Worked examples use made-up numbers for teaching. All values were calculated in Python.
Sources: de Moivre A, Latin pamphlet, 1733 (MacTutor History of Mathematics, biography of Abraham de Moivre). de Oliveira MJ. Rev Bras Ensino Fís 2022;44 (doi:10.1590/1806-9126-rbef-2022-0087). Gauss CF. Theoria motus corporum coelestium, 1809. Magnus JR. Gauss on least-squares and maximum-likelihood estimation. SSRN 2022 (doi:10.2139/ssrn.3990758). Pearson K. Phil Trans R Soc Lond A 1894;185:71–110 (doi:10.1098/rsta.1894.0003); 1895;186:343–414 (doi:10.1098/rsta.1895.0010); 1896;187:253–318 (doi:10.1098/rsta.1896.0007). Pearson K. Biometrika 1905;4:169–212 (doi:10.1093/biomet/4.1-2.169). Fisher RA. Trans R Soc Edinb 1918;52:399–433, printed volume dated 1919 (doi:10.1017/S0080456800012163). Nagele P. Br J Anaesth 2003;90:514–6 (doi:10.1093/bja/aeg087; PMID 12644429). Altman DG, Bland JM. BMJ 2005;331:903 (doi:10.1136/bmj.331.7521.903; PMID 16223828). Lang TA, Altman DG. SAMPL guidelines. In: Science Editors' Handbook, EASE 2013; Int J Nurs Stud 2015;52:5–9 (doi:10.1016/j.ijnurstu.2014.09.006; PMID 25441757). Vandenbroucke JP et al. STROBE explanation and elaboration. PLoS Med 2007;4:e297 (doi:10.1371/journal.pmed.0040297). Hopewell S et al. CONSORT 2025 explanation and elaboration. BMJ 2025;389:e081124 (doi:10.1136/bmj-2024-081124; PMID 40228832). Sandberg S et al. Clin Chem Lab Med 2015;53:833–5 (doi:10.1515/cclm-2015-0067; PMID 25719329).
Added sources: von Hippel PT. Mean, median, and skew: correcting a textbook rule. J Stat Educ 2005;13(2) (doi:10.1080/10691898.2005.11910556). Ozarda Y. Reference intervals: current status, recent developments and future considerations. Biochem Med (Zagreb) 2016;26:5–16 (doi:10.11613/BM.2016.001; PMID 26981015). NIST/SEMATECH e-Handbook of Statistical Methods, Bienaymé–Chebyshev rule (itl.nist.gov/div898/handbook).
Word histories: Miller J. Earliest Known Uses of Some of the Words of Mathematics (MacTutor). Online Etymology Dictionary.
Past papers
Every TNMGRMU question on descriptive statistics
Each one is answered in this chapter. Tap a sitting to open the university's own question paper.
- ×4Standard Deviation
- ×3Measures of Dispersion
- ×2Standard error and standard deviation
- ×2Standard Error & Standard Deviation
- ×2Laboratory results for glucose shows a standard deviation of 12, mean of 101 mg/dL, whereas for creatinine the observed standard deviation is 0.5 and mean is 1.2 mg/dL. The quality manager judges analytical performance for glucose to be better than that of creatinine. Justify
- ×1Principle and types of ANOVA (analysis of variance)
- ×1Write a note on the measures of central tendency and dispersion in biostatistics
- ×1Mode and median
- ×1Median and mean
- ×1Explain ‘Skewed and normal distribution of data’ with relevant examples
- ×1Define the common measures of central tendency for data set
- ×1Mean, mode, median and standard deviation
- ×1Explain mean, standard deviation and coefficient of variation in medical practice with appropriate examples
- ×1Measures of central tendency and their uses
- ×1Concept of standard deviation and percentiles
- ×1Normal distribution (Gaussian curve)
- ×1Write in brief about (a) Mean (b) median and (c) mode
- ×1Explain in brief about skewness and kurtosis
- ×1Explain standard deviation and standard error mean
- ×1It is better to express imprecision of an analytical method in terms of coefficient of variation rather than standard deviation. Justify with an example
- ×1Statistical significance of Standard deviation
- ×1Define normal distribution
- ×1Define standard deviation
- ×1Correlation coefficient and coefficient of variation
- ×1Standard deviations
- ×1Calculation of standard deviation and t-value of small experimental groups. ]
- ×1Write brlefly an : (a) Analytical epidemiology. (b) Standard Deviation. (c) Vitamln A Deficiency. (d) Dimehslons of growth and development
Further reading