PAF Analysis

Setting up

Set seed

Packages

Warning: package 'ggplot2' was built under R version 4.4.3
Loading required package: grid
Loading required package: checkmate
Loading required package: abind
Warning: package 'dplyr' was built under R version 4.4.3

Attaching package: 'dplyr'
The following objects are masked from 'package:stats':

    filter, lag
The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union
Warning: package 'patchwork' was built under R version 4.4.3
Warning: package 'scales' was built under R version 4.4.3

Attaching package: 'scales'
The following object is masked from 'package:purrr':

    discard

Attaching package: 'cowplot'
The following object is masked from 'package:patchwork':

    align_plots
The following object is masked from 'package:gt':

    as_gtable

Attaching package: 'gridExtra'
The following object is masked from 'package:dplyr':

    combine
Warning: package 'ggh4x' was built under R version 4.4.2
Registering fonts with R

Setting working directory

Read in Functions

Read in Data

Warning: package 'slider' was built under R version 4.4.2
Warning: There was 1 warning in `mutate()`.
ℹ In argument: `value = case_when(...)`.
ℹ In group 1: `age_group = "20-49"`, `sex = "Men"`, `level = "(0,10]"`,
  `variable = "redmeat_consumption"`, `chain = 1`.
Caused by warning:
! Calling `case_when()` with size 1 LHS inputs and size >1 RHS inputs was
  deprecated in dplyr 1.2.0.
ℹ This `case_when()` statement can result in subtle silent bugs and is very inefficient.

  Please use a series of if statements instead:

  ```
  # Previously
  case_when(scalar_lhs1 ~ rhs1, scalar_lhs2 ~ rhs2, .default = default)

  # Now
  if (scalar_lhs1) {
    rhs1
  } else if (scalar_lhs2) {
    rhs2
  } else {
    default
  }
  ```

Cleaning childhood BMI data

Ensuring the Joinpoint input is up to date

Main Figures

Figure 3 - Aggregated PAFs by cancer site

`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
ℹ Summaries were computed grouped by year, variable, and sex.
ℹ Output is grouped by year and variable.
ℹ Use `summarise(.groups = "drop_last")` to silence this message.
ℹ Use `summarise(.by = c(year, variable, sex))` for per-operation grouping
  (`?dplyr::dplyr_by`) instead.

Supplementary Figures

Supplementary Table 1 - Data availability and sourcing information

This table is generated in excel

Supplementary Table 2 - Behavioural risk factors

This table is generated in excel

Supplementary Table 3 - Relative risks

This table is generated in excel

Supplementary Table 4 - PubMed search terms

This figure has been created externally

Supplementary Table 5 - AAPC/APC by cancer site

Joinpoint Cancer Incidence Trends in Women
Age Group AAPC Most Recent APC
Time Period AAPC1 CI 95% P Value P Value Difference Time Period APC2 CI 95% P Value
Breast
20-49 (2001, 2019) 0.88 (0.56, 1.21) <0.001 - (2013, 2019) 0.21 (-0.67, 1.09) 0.62
50+ (2001, 2019) 0.70 (0.17, 1.23) 0.01 0.28 (2003, 2019) 0.42 (0.22, 0.61) <0.001
Colorectum
20-49 (2001, 2019) 3.20 (2.78, 3.62) <0.001 - (2001, 2019) 3.20 (2.78, 3.62) <0.001
50+ (2001, 2019) 0.36 (-0.13, 0.85) 0.15 <0.001 (2017, 2019) 2.01 (-1.88, 6.06) 0.28
Endometrium
20-49 (2001, 2019) 2.77 (2.22, 3.32) <0.001 - (2001, 2019) 2.77 (2.22, 3.32) <0.001
50+ (2001, 2019) 1.80 (1.46, 2.14) <0.001 0.003 (2011, 2019) 0.29 (-0.33, 0.91) 0.34
Gallbladder
20-49 (2001, 2019) 3.21 (1.86, 4.58) <0.001 - (2001, 2019) 3.21 (1.86, 4.58) <0.001
50+ (2001, 2019) 3.38 (2.95, 3.81) <0.001 0.59 (2001, 2019) 3.38 (2.95, 3.81) <0.001
Kidney
20-49 (2001, 2019) 4.53 (3.82, 5.25) <0.001 - (2001, 2019) 4.53 (3.82, 5.25) <0.001
50+ (2001, 2019) 3.04 (2.42, 3.66) <0.001 <0.001 (2014, 2019) -0.65 (-2.64, 1.38) 0.50
Liver
20-49 (2001, 2019) 5.08 (1.61, 8.66) 0.004 - (2003, 2019) 2.89 (1.66, 4.15) <0.001
50+ (2001, 2019) 3.39 (1.72, 5.10) <0.001 0.20 (2013, 2019) 0.70 (-1.59, 3.05) 0.52
Multiple Myeloma
20-49 (2001, 2019) 6.02 (4.93, 7.13) <0.001 - (2001, 2019) 6.02 (4.93, 7.13) <0.001
50+ (2001, 2019) 2.62 (2.13, 3.11) <0.001 <0.001 (2001, 2019) 2.62 (2.13, 3.11) <0.001
Oral
20-49 (2001, 2019) 2.15 (1.50, 2.80) <0.001 - (2001, 2019) 2.15 (1.50, 2.80) <0.001
50+ (2001, 2019) 2.70 (1.48, 3.93) <0.001 0.78 (2010, 2019) 1.61 (0.93, 2.30) <0.001
Ovary
20-49 (2001, 2019) 0.72 (0.05, 1.39) 0.04 - (2013, 2019) -0.57 (-2.33, 1.23) 0.51
50+ (2001, 2019) -1.31 (-1.51, -1.10) <0.001 <0.001 (2001, 2019) -1.31 (-1.51, -1.10) <0.001
Pancreas
20-49 (2001, 2019) 2.63 (1.86, 3.41) <0.001 - (2001, 2019) 2.63 (1.86, 3.41) <0.001
50+ (2001, 2019) 0.98 (0.67, 1.28) <0.001 <0.001 (2013, 2019) -0.09 (-0.89, 0.72) 0.81
Thyroid
20-49 (2001, 2019) 6.12 (5.51, 6.73) <0.001 - (2013, 2019) 3.27 (1.69, 4.88) <0.001
50+ (2001, 2019) 4.77 (3.94, 5.60) <0.001 0.005 (2014, 2019) 2.36 (-0.32, 5.12) 0.08
1 Average Annual Percentage Change
2 Annual Percentage Change
Joinpoint Cancer Incidence Trends in Men
Age Group AAPC Most Recent APC
Time Period AAPC1 CI 95% P Value P Value Difference Time Period APC2 CI 95% P Value
Colorectum
20-49 (2001, 2019) 2.60 (2.26, 2.93) <0.001 - (2001, 2019) 2.60 (2.26, 2.93) <0.001
50+ (2001, 2019) -0.21 (-0.79, 0.37) 0.48 <0.001 (2014, 2019) -0.58 (-1.37, 0.21) 0.13
Gallbladder
20-49 (2001, 2019) 5.87 (-0.68, 12.86) 0.08 - (2003, 2019) 2.29 (-0.05, 4.68) 0.05
50+ (2001, 2019) 2.68 (2.08, 3.29) <0.001 0.17 (2001, 2019) 2.68 (2.08, 3.29) <0.001
Kidney
20-49 (2001, 2019) 3.69 (2.18, 5.22) <0.001 - (2016, 2019) -0.67 (-7.32, 6.46) 0.84
50+ (2001, 2019) 2.79 (2.10, 3.49) <0.001 0.15 (2014, 2019) 0.05 (-2.19, 2.34) 0.96
Liver
20-49 (2001, 2019) 1.47 (-0.93, 3.93) 0.23 - (2014, 2019) -3.52 (-10.98, 4.56) 0.36
50+ (2001, 2019) 4.52 (3.86, 5.20) <0.001 0.99 (2014, 2019) 1.49 (-0.68, 3.70) 0.16
Multiple Myeloma
20-49 (2001, 2019) 4.07 (2.99, 5.16) <0.001 - (2001, 2019) 4.07 (2.99, 5.16) <0.001
50+ (2001, 2019) 2.58 (2.15, 3.02) <0.001 0.009 (2001, 2019) 2.58 (2.15, 3.02) <0.001
Oral
20-49 (2001, 2019) 1.27 (0.85, 1.70) <0.001 - (2013, 2019) -1.07 (-2.17, 0.05) 0.06
50+ (2001, 2019) 2.92 (2.27, 3.58) <0.001 >0.999 (2013, 2019) 1.75 (0.03, 3.50) 0.05
Pancreas
20-49 (2001, 2019) 1.13 (-0.55, 2.84) 0.19 - (2004, 2019) 2.34 (1.47, 3.22) <0.001
50+ (2001, 2019) 1.06 (0.88, 1.24) <0.001 0.47 (2001, 2019) 1.06 (0.88, 1.24) <0.001
Thyroid
20-49 (2001, 2019) 6.04 (4.81, 7.27) <0.001 - (2012, 2019) 3.79 (1.17, 6.48) 0.008
50+ (2001, 2019) 4.91 (3.52, 6.33) <0.001 0.12 (2014, 2019) 1.98 (-2.51, 6.67) 0.37
1 Average Annual Percentage Change
2 Annual Percentage Change

Supplementary Table 6 - AAPC/APC by risk factor

Joinpoint Trends in Risk Factors Women
Age Group AAPC Most Recent APC
Time Period AAPC1 CI 95% P Value Time Period APC2 CI 95% P Value
Cigarette Smoking - Former
20-49 (1995, 2019) -0.04 (-0.69, 0.61) 0.90 (2017, 2019) -6.12 (-9.43, -2.69) 0.003
50+ (1995, 2019) 0.13 (-0.07, 0.32) 0.20 (2010, 2019) -0.25 (-0.68, 0.18) 0.24
Cigarette Smoking - Current
20-49 (1995, 2019) -2.33 (-2.60, -2.07) <0.001 (2011, 2019) -2.23 (-2.70, -1.76) <0.001
50+ (1995, 2019) -2.23 (-2.61, -1.86) <0.001 (2016, 2019) 0.34 (-0.88, 1.58) 0.55
Alcohol Consumption - Light Drinker3
20-49 (2011, 2019) 0.19 (-0.25, 0.64) 0.39 (2016, 2019) -1.03 (-2.35, 0.30) 0.10
50+ (2011, 2019) 0.44 (0.24, 0.64) 0.001 (2011, 2019) 0.44 (0.24, 0.64) 0.001
Alcohol Consumption - Moderate Drinker3
20-49 (2011, 2019) -3.02 (-4.25, -1.77) <0.001 (2017, 2019) -0.46 (-6.69, 6.17) 0.85
50+ (2011, 2019) -0.27 (-0.65, 0.11) 0.14 (2011, 2019) -0.27 (-0.65, 0.11) 0.14
Alcohol Consumption - Heavy Drinker3
20-49 (2011, 2019) -2.85 (-5.00, -0.65) 0.01 (2014, 2019) 2.97 (-0.10, 6.14) 0.06
50+ (2011, 2019) 0.81 (-0.61, 2.26) 0.22 (2011, 2019) 0.81 (-0.61, 2.26) 0.22
Body Mass Index (BMI) - Overweight4
20-49 (1995, 2019) 0.01 (-0.23, 0.26) 0.93 (2016, 2019) 0.82 (0.03, 1.61) 0.04
50+ (1995, 2019) -0.72 (-0.87, -0.56) <0.001 (2012, 2019) -1.29 (-1.75, -0.82) <0.001
Body Mass Index (BMI) - Obese4
20-49 (1995, 2019) 2.61 (2.33, 2.88) <0.001 (2013, 2019) 3.23 (2.48, 3.99) <0.001
50+ (1995, 2019) 1.61 (1.00, 2.22) <0.001 (2013, 2019) 1.90 (1.15, 2.66) <0.001
Physical Inactivity - Below UK Recommendations5
20-49 (2003, 2012) -1.44 (-2.43, -0.44) 0.02 (2003, 2012) -1.44 (-2.43, -0.44) 0.02
50+ (2003, 2012) -1.25 (-1.88, -0.61) 0.008 (2003, 2012) -1.25 (-1.88, -0.61) 0.008
Fibre Intake Deficiency - Below UK Recommendations6
20-49 (2008, 2018) -0.20 (-0.45, 0.06) 0.13 (2012, 2018) -0.55 (-0.88, -0.22) 0.007
50+ (2008, 2018) -0.43 (-0.68, -0.18) <0.001 (2014, 2018) -0.02 (-0.63, 0.59) 0.94
Red Meat Consumption - Median
20-49 (2008, 2018) -7.06 (-12.04, -1.81) 0.009 (2015, 2018) -17.80 (-32.31, -0.19) 0.05
50+ (2008, 2018) -7.39 (-8.64, -6.13) <0.001 (2008, 2018) -7.39 (-8.64, -6.13) <0.001
Processed Meat Consumption - Median
20-49 (2008, 2018) -5.36 (-7.88, -2.77) <0.001 (2015, 2018) -10.99 (-19.09, -2.09) 0.02
50+ (2008, 2018) -5.55 (-7.37, -3.70) <0.001 (2008, 2018) -5.55 (-7.37, -3.70) <0.001
1 Average Annual Percentage Change
2 Annual Percentage Change
3 Heavy (>50g/day), moderate (12-50 g/day), light (<12 g/day)
4 Overweight (>30 kg/m²), obese (25-30 kg/m²)
5 Recommendations are 150 min/week of moderate or 75 min/week or vigorous activity
6 Guidelines are 30g of fibre per day
Joinpoint Trends in Risk Factors Men
Age Group AAPC Most Recent APC
Time Period AAPC1 CI 95% P Value Time Period APC2 CI 95% P Value
Cigarette Smoking - Former
20-49 (1995, 2019) -0.01 (-0.33, 0.32) 0.96 (2017, 2019) -3.54 (-5.93, -1.09) 0.009
50+ (1995, 2019) -1.07 (-1.31, -0.82) <0.001 (2015, 2019) -0.11 (-1.07, 0.85) 0.81
Cigarette Smoking - Current
20-49 (1995, 2019) -1.63 (-2.04, -1.22) <0.001 (2016, 2019) -1.26 (-2.43, -0.09) 0.04
50+ (1995, 2019) -2.25 (-2.62, -1.88) <0.001 (2013, 2019) -3.71 (-4.68, -2.74) <0.001
Alcohol Consumption - Light Drinker3
20-49 (2011, 2019) 1.30 (0.70, 1.91) <0.001 (2016, 2019) 0.35 (-1.45, 2.18) 0.62
50+ (2011, 2019) 0.61 (0.06, 1.16) 0.03 (2017, 2019) -1.31 (-3.99, 1.44) 0.25
Alcohol Consumption - Moderate Drinker3
20-49 (2011, 2019) -2.17 (-2.90, -1.44) <0.001 (2017, 2019) -0.75 (-4.44, 3.08) 0.61
50+ (2011, 2019) -1.10 (-1.56, -0.64) <0.001 (2017, 2019) 0.69 (-1.66, 3.09) 0.47
Alcohol Consumption - Heavy Drinker3
20-49 (2011, 2019) -3.56 (-5.03, -2.06) <0.001 (2016, 2019) -1.10 (-5.61, 3.63) 0.55
50+ (2011, 2019) 0.25 (-0.27, 0.77) 0.29 (2011, 2019) 0.25 (-0.27, 0.77) 0.29
Body Mass Index (BMI) - Overweight4
20-49 (1995, 2019) -0.41 (-0.67, -0.14) 0.003 (2015, 2019) 0.39 (-0.49, 1.28) 0.37
50+ (1995, 2019) -0.39 (-0.57, -0.20) <0.001 (2007, 2019) -0.24 (-0.40, -0.09) 0.005
Body Mass Index (BMI) - Obese4
20-49 (1995, 2019) 2.19 (1.89, 2.50) <0.001 (2001, 2019) 0.79 (0.58, 1.01) <0.001
50+ (1995, 2019) 2.11 (1.69, 2.53) <0.001 (2009, 2019) 0.50 (0.21, 0.80) 0.002
Physical Inactivity - Below UK Recommendations5
20-49 (2003, 2012) -1.86 (-3.51, -0.18) 0.04 (2003, 2012) -1.86 (-3.51, -0.18) 0.04
50+ (2003, 2012) -0.73 (-0.93, -0.53) 0.001 (2003, 2012) -0.73 (-0.93, -0.53) 0.001
Fibre Intake Deficiency - Below UK Recommendations6
20-49 (2008, 2018) -0.38 (-1.34, 0.58) 0.43 (2016, 2018) 2.19 (-3.31, 8.00) 0.37
50+ (2008, 2018) -0.30 (-0.50, -0.10) 0.007 (2008, 2018) -0.30 (-0.50, -0.10) 0.007
Red Meat Consumption - Median
20-49 (2008, 2018) -7.41 (-8.92, -5.87) <0.001 (2016, 2018) -19.99 (-27.17, -12.10) 0.001
50+ (2008, 2018) -4.40 (-8.10, -0.54) 0.03 (2016, 2018) -13.92 (-31.34, 7.91) 0.16
Processed Meat Consumption - Median
20-49 (2008, 2018) -4.88 (-7.38, -2.32) <0.001 (2012, 2018) -7.77 (-10.91, -4.52) 0.001
50+ (2008, 2018) -2.03 (-5.33, 1.38) 0.24 (2013, 2018) -9.97 (-15.25, -4.36) 0.005
1 Average Annual Percentage Change
2 Annual Percentage Change
3 Heavy (>50g/day), moderate (12-50 g/day), light (<12 g/day)
4 Overweight (>30 kg/m²), obese (25-30 kg/m²)
5 Recommendations are 150 min/week of moderate or 75 min/week or vigorous activity
6 Guidelines are 30g of fibre per day

Supplementary Figure 1 - AAPC for all cancers

NOTE: Kaposi Sarcoma is removed from the womens APC plot as neither 20-49 or 50+ generated an AAPC due to zeroes in the data

Supplementary Figure 2 - Flowchart

This figure has been created externally

Supplementary Figure 4 - PAFs by risk factor

`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
ℹ Summaries were computed grouped by year, variable, and sex.
ℹ Output is grouped by year and variable.
ℹ Use `summarise(.groups = "drop_last")` to silence this message.
ℹ Use `summarise(.by = c(year, variable, sex))` for per-operation grouping
  (`?dplyr::dplyr_by`) instead.

Scale for y is already present.
Adding another scale for y, which will replace the existing scale.

Supplementary Figure 5 - Attributable Graphs for 50+

`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
`summarise()` has regrouped the output.
ℹ Summaries were computed grouped by year, variable, and sex.
ℹ Output is grouped by year and variable.
ℹ Use `summarise(.groups = "drop_last")` to silence this message.
ℹ Use `summarise(.by = c(year, variable, sex))` for per-operation grouping
  (`?dplyr::dplyr_by`) instead.

Supplementary Figure 6 - PAF

Warning: The `size` argument of `element_rect()` is deprecated as of ggplot2 3.4.0.
ℹ Please use the `linewidth` argument instead.
Scale for y is already present.
Adding another scale for y, which will replace the existing scale.

Supplementary Figure 7 - Burden of cancer rates by age group

Supplementary Figure 8 - Colorectal cancer and risk factor plotting

Ignoring unknown labels:
• fill : ""
• linetype : ""
Scale for colour is already present.
Adding another scale for colour, which will replace the existing scale.
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Scale for colour is already present.
Adding another scale for colour, which will replace the existing scale.

Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""
Ignoring unknown labels:
• fill : ""
• linetype : ""