Skip to main content

Chapter 13 Comparing quantitative data within individuals

described in detail following the image
A diagram showing the six steps of research, arranged in a cycle: Ask, Design, Collect, Classify and summarise, Analyse, Report. Step 4 (Summarise) is highlighted, captioned β€œComparing quantitative data within individuals”.

Section 13.1 Introduction

Sometimes the same quantitative variable is measured on each individual more than once (i.e., within-individual changes for each unit of analysis) but only a small number of times. Examples of this type of data include:
  • measurements of weekly household water consumption for many households, before and after installing water-saving devices.
  • blood pressure recorded for people at \(8\)am, \(1\)pm and \(8\)pm each day.
In both cases, the same variable is measured multiple times for each individual. This chapter studies summarising within-individuals changes in quantitative variables.

Section 13.2 Numerical summary: mean differences

When each individual has two observations, the difference between the observations can be computed for each individual. Then, the appropriate numerical summary is to summarise these differences over all individuals; for example, using the mean of these differences.
When more than two observations are made for each individual, the changes from the first observation (the reference level) may be computed (e.g., if the first observation is pre-intervention, or a benchmark, observation); for example, using the mean change (see Sect. SectionΒ 13.5).

Example 13.1. Within-individual comparisons.

[111] studied children with atopic asthma, and measured the immunogoblin E concentrations (IgE) before and after an intervention for each child (Table TableΒ 13.2), plus the reduction in IgE for each child. The child is the individual.
For the IgE data, the numerical summary table is shown in Table TableΒ 13.3. The direction of the difference is implied by the word ’reduction’.
Table 13.2. The IgE before and after an intervention, and the change in IgE (in micrograms/L).
IgE (before) in \(\mu\)g/L IgE (after) in \(\mu\)g/L IgE reduction in \(\mu\)g/L
83 83 0
292 292 0
293 292 1
623 542 81
792 709 83
1543 1000 543
1668 1000 668
1960 1626 334
2877 2502 375
2961 2711 250
5504 4504 1000
Table 13.3. A numerical summary of the IgE data (in \(\mu\)g/L).
Mean Std dev. Sample size
Before 1690.5 1615.53 11
After 1387.4 1354.28 11
Reduction 303.2 325.28 11

Warning 13.4.

In the numerical summary table, the information for the differences is not found by subtracting the information in one row from the other. In Table TableΒ 13.3, for example, the number of differences is not \(11 - 11 = 0\text{;}\) the standard deviation of the differences is not \(1615.53 - 1354.28 = 261.25\,\mu\text{g}\,\text{L}^{-1}\text{.}\) These statistics are computed from the differences (i.e., the Reductions in Table TableΒ 13.2).

Section 13.3 Graphs for the differences

For within-individual changes for a quantitative variable, options for plotting include:

Subsection 13.3.1 Histogram of differences

Sometimes the same variable is measured on each unit of analysis twice, when the changes (or differences) for each individual can be produced, and a histogram of the changes or differences can be constructed. The direction of the differences should be clear (e.g., first measurement minus second, or second measurement minus first).
Issues relevant for constructing histograms (Sect. SubsectionΒ 11.3.1), such as bin widths and boundary values, also apply here.

Warning 13.5.

The axis displaying the counts (or percentages) should start from zero, since the height of the bars visually implies the frequency of those differences (see Example ExampleΒ 17.5).

Subsection 13.3.2 Case-profile plots

Sometimes the variable is measured or recorded more than twice, and so a single set of differences cannot be produced. In these cases, the values for each individual can be plotted using a case-profile plot: the measurements are shown on one axis (usually the vertical), and the various points at which measurements are taken are shown on the other axis. A case-profile plot is still useful for paired data, of course.

Warning 13.7.

The axis displaying the counts (or percentages) need not start from zero, since the distance from the axis to the lines do not visually imply any quantity of interest. Rather, the relative changes represented by the lines display the quantity of interest.

Example 13.8. Case-profile plot.

For the IgE data (Table TableΒ 13.2), the measurements of IgE for each child at both times can be shown in a case-profile plot (Fig. FigureΒ 13.9, right panel).
Each line corresponds to a unit of analysis (i.e., a child).
Figure 13.9. The IgE data. Left: a histogram of the differences (reductions); boundary observations are counted in the lower box. Right: a case-profile plot, where each line joins each person’s pre-intervention to their post-intervention measurement.

Example 13.10. Case-profile plot.

Runners use wearable devices to measure many performance indicators, including vertical oscillation (VO). VO contributes to running economy and injury risk, so reliable VO measurements are crucial. [105] compared four devices, and obtained data from video analysis for \(n = 150\) athletes; that is, each participant had the same runs measured using five methods. The case-profile plot (Fig. FigureΒ 13.23) shows the means for each method using a solid point; each line represents one runner. NOVA and Footpod give smaller VO measurements in general.
As in Example ExampleΒ 13.10, the case-profile plot is hard to read with large numbers of individuals, and so sometimes the mean (or median, as appropriate) is shown, with some measure of the variation of the observations (Fig. FigureΒ 13.23 shows the minimum and maximum values for each method, for instance).

Section 13.4 Example: invasive plants

Skypilot (Polemonium viscosum) is a native alpine wildflower growing in the Colorado Rocky Mountains (USA). In recent years, a willow shrub (Salix) has been encroaching on skypilot territory and, because willow often flowers early, researchers [117] are concerned that the willow may ’negatively affect pollination regimes of resident alpine wildflower species’ (p. \(6\,965\)). One RQ was:
In the Colorado Rocky Mountains, what is the mean difference between first-flowering day for the native skypilot and the encroaching willow?
Data for both species was collected at \(25\) different sites (Table TableΒ 13.11).
The site is the individual; the data are paired (Sect. SectionΒ 29.1), a form of blocking (Sect. SectionΒ 7.2). The ’first-flowering day’ is the number of days since the start of the year (e.g., January 12 is ’day 12’) when flowers were first observed.
Table 13.11. The day of the year of first flowering by encroaching willow and native skypilot.
Site Willow Skypilot
1 201 201
2 178 179
3 189 189
4 189 189
5 196 203
6 207 203
7 199 199
8 178 182
9 178 178
10 191 191
11 187 192
12 190 197
13 190 190
14 209 209
15 221 221
16 179 188
17 174 179
18 172 166
19 196 196
20 173 173
21 180 173
22 181 179
23 186 186
24 194 209
25 197 197
Since the data are available, the data should be summarised graphically (Fig. FigureΒ 13.14) and numerically (Table TableΒ 13.13), using software output (Fig. FigureΒ 13.12).
Figure 13.12. Software output for the flowering-day data.
Table 13.13. The day of first flowering for encroaching willow and native skypilot.
Mean Std dev. Sample size
Willow (encroaching) 189.40 12.20 25
Skypilot (native) 190.76 13.06 25
Differences 1.36 4.70 25
Figure 13.14. The flowering-day data. Left: a histogram of the difference between the first-flowering days (skypilot minus willow). Right: a case-profile plot of days of first flowering (unfilled points and dashed lines indicate earlier or same dates for willow).

Section 13.5 Example: pain-relieving tape

[118] studied the effect of using Kinesio Tape to alleviate pain in athletes. Pain was measured by applying a slow constant rate of pressure on the left arm, and subjects pressed a button when the sensation moved from pressure to pain. The pressure at which this occurred was recorded. This was repeated \(5\) minutes before applying the tape, \(5\) minutes after applying the tape, and again \(15\)–\(20\) minutes after applying the tape.
Figure FigureΒ 13.15 shows the reported pain for \(16\) subjects. A summary table is shown in Table TableΒ 13.16. The pain thresholds are increasing slightly as time progresses.
Figure 13.15. Pain threshold (left arm) at three time points when using Kinesio Tape, without applying tension, for \(n = 16\) subjects. The filled, black points represent the means for each time point.
Table 13.16. A numerical summary of the tape data: pain thresholds in kPa.
Mean (kPa) Std dev. (kPa) Sample size Mean increase from Pre SD of increase from Pre
Pre: 5 min 446.5 175.18 16
Post: 5 min 479.6 199.61 16 33.0 73.93
Post: 15–20 min 506.9 214.36 16 60.3 102.72

Section 13.6 Chapter summary

Quantitative data measured within individuals can be summarised using a histogram of differences when the variable is measured or observed twice, or a case-profile plot (with two or more measurement or observations per individual). A summary table should show the numerical summaries for the quantitative variable at each measurement or observation and for appropriate changes.

Reading Questions 13.7 Quick review questions

Are the following statements true or false?

1.

A histogram of the differences is only appropriate for within two within-individuals measurements or observations.
Answer.

2.

A case-profile plot is only appropriate for showing changes for two within-individuals measurements or observations.
Answer.
False; a case-profile plot can be used for two or more within-individual comparisons.

3.

The median and IQR are not appropriate for summarising differences.
Answer.
False; use whatever numerical summaries are appropriate.

Exercises 13.8 Exercises

Answers to odd-numbered exercises are given at the end of the book.

1.

[Dataset: Insulation] The Electricity Council in Bristol wanted to determine if a certain type of wall-cavity insulation reduced average energy consumption in winter [112][113]:
In Bristol homes, what is the mean reduction in energy consumption after adding home insulation?
  1. What are the individuals (units of analysis)?
  2. Explain why this study uses a within-individuals comparison.
  3. Use the collected data (shown below) to sketch a case-profile plot.
  4. Use the data to sketch a histogram of the differences.
  5. Use software or a calculator to prepare a summary table.
Table 13.17. The house insulation data: energy consumption before and after adding insulation, and the energy saving (all in MWh).
Home Before After Saving
A 12.1 12.0 0.1
B 11.0 10.6 0.4
C 14.1 13.4 0.7
D 13.8 11.2 2.6
E 15.5 15.3 0.2
F 12.2 13.6 \(-\)1.4
G 12.8 12.6 0.2
H 9.9 8.8 1.1
I 10.8 9.6 1.2
J 12.7 12.4 0.3

2.

[Dataset: Captopril] In a study of hypertension [113][114], \(15\) patients were given a drug (Captopril) and their systolic blood pressure measured (in mm Hg) immediately before and two hours after being given the drug (Table TableΒ 13.18).
  1. Explain why this study uses a within-individuals comparison.
  2. Construct a histogram of the differences.
  3. Construct a case-profile plot for the data.
Table 13.18. The Captopril data: before and after systolic blood pressures (in mm Hg).
Before After Differences
210 201 9
169 165 4
187 166 21
160 157 3
167 147 20
176 145 31
185 168 17
206 180 26
173 147 26
146 136 10
174 151 23
201 168 33
198 179 19
148 129 19
154 131 23

3.

[Dataset: PainRelief] [110] measured the reported pain of new mothers in Dodoma (Tanzania) at four times: near giving birth, then \(20\text{,}\) \(40\) and \(60\) minutes after giving birth. Mothers were administered either paracetamol or a cold pack as pain relief. Pain was recorded using a ’numeric rating scale represented by the horizontal line marked from zero to ten’, where higher scores mean greater pain.
Since the number of individuals is large (\(n = 912\)), use the summary data in Table TableΒ 13.19 to sketch a plot of the means and the range, like that in Figure FigureΒ 13.15.
Table 13.19. A summary table of reported pain for mothers after giving birth.
At birth After 20 min After 40 min After 60 min
Paracetamol Mean 7.44 6.89 4.69 2.84
(\(n = 456\)) Standard deviation 2.01 1.83 1.49 1.19
Minimum 2.00 2.00 2.00 0.00
Maximum 10.00 10.00 9.00 7.00
Cold pack Mean 8.63 5.67 3.19 0.99
(\(n = 455\)) Standard deviation 1.40 2.03 1.63 0.99
Minimum 4.00 0.00 0.00 0.00
Maximum 10.00 9.00 6.00 4.00

4.

[Dataset: Stress] The concentration of beta-endorphins in the blood is a sign of stress. One study ([113], Dataset 232; [116]) measured the beta-endorphin concentration for \(19\) patients about to undergo surgery.
Each patient had their beta-endorphin concentrations measured \(12\)–\(14\) hours before surgery, and also \(10\) minutes before surgery. A numerical summary (from software output) is in Table TableΒ 13.20.
  1. Explain why this study uses a within-individuals comparison.
  2. Explain why the standard deviation for the increase is not the difference between the two individuals time-point standard deviations.
  3. Using the data file and software, construct a histogram of the differences.
  4. Using the data file and software, construct a case-profile plot for the data.
Table 13.20. The numerical summary for the presurgical stress data.
Mean Std deviation Sample size
12–14 h before surgery 8.35 4.397 19
10 min before surgery 16.05 12.509 19
Increase 7.70 13.519 19

6.

[Dataset: WCTennis] [108] recorded the push-time for French wheelchair tennis players, while holding and not holding a racquet (Table TableΒ 13.21; [109]).
  1. What do the differences mean (as given in the table)?
  2. Create a plot of the data.
  3. Create a numerical summary table for the data.
Table 13.21. The wheelchair-tennis data. One observation is missing.
Person Push-time with racquet (s) Push-time without racquet (s) Difference (s)
01 0.2625 0.1833 0.0792
02 0.2375 0.2250 0.0125
03 0.2583 0.2042 0.0542
04 0.1917 0.1875 0.0042
05 0.1875 0.1708 0.0167
06 0.2542 0.1750 0.0792
07 0.2333 0.1917 0.0417
08 0.1917 0.1708 0.0208
09 0.2208 0.2208 0.0000
10 0.2583 0.2750 \(-\)0.0167
11 0.2083 0.1750 0.0333
12 --- 0.2042 ---
13 0.2208 0.2292 \(-\)0.0083

7.

[Dataset: Jumping] [115] recorded double-legged jumping distance for \(80\) healthy people, when they wore shoes and were barefoot (data below).
  1. What do the differences mean (as given in the table)?
  2. Create a plot of the data.
  3. Create a numerical summary table for the data.
Table 13.22. Jumping distances for athletes, with and without shoes (the first four and the last four observations).
With shoes (cm) Barefoot (cm) Difference (cm)
42.73 42.23 0.50
41.00 39.47 1.53
27.37 30.40 \(-\)3.03
46.80 36.60 10.20
\(\vdots\) \(\vdots\) \(\vdots\)
32.73 33.90 \(-\)1.17
56.50 55.10 1.40
33.57 32.07 1.50
27.77 33.57 \(-\)5.80

Section 13.9 Supplementary figures

Figure 13.23. Vertical oscillation (VO) measured using five methods for \(150\) runners. The solid black points represent the means for each method. Left: a line is plotted for each individual. Right: only the means are shown, with vertical lines from the minimum value to the maximum value for each method.