$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The PBD mean, sum, and coverage uniquely provide information about pain responses not captured in other standardized pain scales. Between the two PBDs (Figure 2A,B), the mean pain intensity is identical (PBD mean = 79.6). An increased coverage and sum, however, reveals the greater spatial spread of pain and total pain intensity, respectively, that differentiate the two PBDs (Figure 2B). To accurately quantify pain using these metrics, researchers should avoid the following common PBD setup mistakes (Figure 2C). Excessively large pen thickness and extraneous elements outside the body outline, such as circling body regions or written descriptors will not be captured in the PBD processing. Similarly, a white pen used to remove color rather than the eraser tool will skew PBD metrics. Practice and reinforced instruction will empower patients to create accurate and quantifiable PBDs that reveal variability in pain intensity and distribution.
The PBD metrics were validated against the NRS, VAS, and MPQ (Figure 3B; Supplementary Figure 2) and scored high in usability (Supplementary Figure 1 and Supplementary Figure 2).
PBD metrics correlated to standard pain metrics
The PBD metrics were correlated with the NRS, VAS, and MPQ for most patients (Figure 3A, Supplementary Figure 1A,B). In four of five patients, the PBD sum, coverage, and mean were correlated to their VAS and NRS (Spearman's correlation, rs = 0.33-0.72, p < 0.004, Supplementary Table 1). For three out of five participants, PBD metrics were also significantly correlated with MPQ scores (Spearman's correlation, rs = 0.38-0.53, p < 0.004, Supplementary Table 1). However, patient 4 did not show significant correlations between the PBD metrics and standard pain scores. We further characterized non-linear relationships between PBD and standard metrics using information theory analyses (Supplementary Figure 2).
PBD metrics avoid response anchoring and share mutual information with standard pain metrics
PBD metrics contained more information (i.e., entropy) than the NRS. Across patients, NRS contained less information (2.32 ± 0.37 bits) compared to VAS intensity, VAS unpleasantness, MPQ total, PBD sum, PBD coverage, and PBD mean (3.21 ± 0.49 bits, 3.20 ± 0.31 bits, 3.16 ± 0.23 bits, 3.06 ± 0.32 bits, 3.34 ± 0.16 bits, 3.22 ± 0.39 bits, respectively; Supplementary Figure 2). This was confirmed with a one-way repeated measures ANOVA (F(4,1) = 12.10, p < 0.05) and a Tukey's t-test for individual comparisons (all p < 0.05). This shows PBD metrics had less response anchoring than the NRS.
The PBD was further validated against established metrics by mutual information analyses (permutation testing, α=0.05). In four of five patients, PBD metrics significantly shared MI with the NRS, VAS intensity, VAS unpleasantness, and the MPQ (p < 0.05, Figure 3B). In contrast, patient 4's PBD metrics did not significantly share MI with established metrics. Since their NRS contained the least information across patients' (Supplementary Figure 2), this suggests the NRS failed to capture nuances in pain experience that were captured by the PBD. In all patients, the NRS shared significant MI with VAS intensity, VAS unpleasantness, and MPQ while the PBD sum shared MI with PBD coverage and PBD mean (p < 0.05, Figure 3B). Altogether, for most patients, the PBD metrics shared MI with established pain metrics.
PBDs were easy to use for most participants
In the study, four of the five patients found the PBD easy to use and to accurately reflect their pain (Supplementary Table 2). However, patient 4 reported that the PBD was difficult to use (5 on a 5-point Likert Scale). This is primarily because they have deep, visceral pain-which is not well-captured in a 2-dimensional (2D) PBD. While patients varied in their familiarity with PBDs (2.8 ± 1.2, range 1-4, 5-point Likert Scale), they all used comparable electronics daily (5.0 ± 0.0, 5-point Likert Scale) and found the PBD to be user-friendly (5.2 ± 0.4, range 5-6, 6-point Likert Scale).

Figure 1. Pain body diagram (PBD) analysis workflow. Patients drew on blank PBD templates to represent the pain's location and intensity. Completed PBDs contained hues that ranged from green to blue to red, representing mild to moderate to severe pain regions, respectively. PBDs were masked to include only pixels within the body outline and then the template was removed to isolate only pixels containing hues. From the PBDs, PBD coverage (%), sum intensity (normalized to 0-100), and mean intensity (normalized to 0-100) were calculated. For PBD coverage, the number of colored pixels were first divided by the total number of pixels within the diagram (820,452 pixels for females, 724,608 pixels for males), then multiplied by 100. For PBD sum intensity, the hue values for all pixels in the body diagram were first summed (female range: 0-114,453,054; male range: 0-101,082,816). The sum was then divided by the maximum PBD sum intensity (females: 820,452 pixels multiplied by maximum hue value 139.5, males: 724,608 pixels by 139.5) and multiplied by 100. For PBD mean intensity, the sum of all hue values was divided by the total number of colored pixels, then normalized by dividing by the maximum hue value of 139.5. Please click here to view a larger version of this figure.

Figure 2. Representative PBDs showing examples of good and bad PBDs. (A,B) Good PBDs show the utility of calculating 3 pain metrics. (C) Bad PBD examples include excessively thick pen size, extraneous elements outside the body diagram, and inaccurate erasing. Please click here to view a larger version of this figure.

Figure 3. PBD metrics were validated against standard pain metrics via Spearman's correlation and mutual information analyses. (A) VAS intensity and PBD sum plotted with linear best-fit lines drawn for each patient. (B) Group-level data showing the mean mutual information (MI) between each pain metric, with MI indicated by color bar on the right. The text in each box represents the number of patients with statistically significant MI for a given pairwise comparison (e.g., 3/5 indicates 3 patients with significant values). MI is presented by the observed MI divided by the theoretical max MI. Abbreviations: NRS=numeric rating scale; VAS intensity = visual analog scale intensity; VAS unpl. = visual analog scale pain unpleasantness, MPQ=short form McGill pain questionnaire 2; PBD=pain body diagram; PBD cov. = PBD coverage, MI = mutual information, sig. = significant. Please click here to view a larger version of this figure.
Supplementary Figure 1. PBD mean (A) and PBD coverage (B) plotted against VAS intensity with linear best-fit lines drawn for each patient. Abbreviations: VAS=visual analog scale; PBD=pain body diagram. Please click here to download this File.
Supplementary Figure 2. Entropy per pain metric across patients. On the group-level, NRS intensity had lower entropy than every other pain metric as shown by a repeated measures one-way ANOVA followed by Tukey's test post-hoc for specific comparisons * = p < 0.05, ** = p < 0.001. Abbreviations: NRS=numeric rating scale; VAS=visual analog scale; MPQ=McGill pain questionnaire; PBD=pain body diagram. Please click here to download this File.
Supplementary Table 1. Spearman's correlations between PBD metrics and self-reported standard pain measures. Spearman's correlation coefficients (rho) for three extracted PBD metrics against NRS, VAS, and MPQ pain measures. Abbreviations: NRS=numeric rating scale; VAS=visual analog scale; MPQ=McGill pain questionnaire; PBD=pain body diagram. Please click here to download this File.
Supplementary Table 2. Patient impressions of completing a PBD were revealed through PBD-specific and system usability scale-modified questions. The modified usability scale questions alternated in positive and negative statements and were ranked on a 5-point scale (1=strongly agree, 5=strongly disagree). Abbreviation: PBD=pain body diagram. Please click here to download this File.
Supplementary Coding File 1: Python script for PBD metrics. The annotated python code processes a pain body diagram PNG file and outputs PBD mean, coverage, and sum values for each file. The script also includes import statements to download the required packages for the program to run. Please click here to download this File.
Supplementary File 1: Supplementary file for methodological details. Please click here to download this File.