This method computes the proportional change of absolute (rate differences) and relative (rate ratios) inequalities of prevalence rates for two different status groups, as proposed by Mackenbach et al. (2015).
inequ_trend(data, prev.low, prev.hi)Arguments
- data
A data frame that contains the variables with prevalence rates for both low and high status groups (see 'Examples').
- prev.low
The name of the variable with the prevalence rates for the low status groups.
- prev.hi
The name of the variable with the prevalence rates for the hi status groups.
Value
A data frame with the prevalence rates as well as the values for the
proportional change in absolute (rd) and relative (rr)
ineqqualities.
Details
Given the time trend of prevalence rates of an outcome for two status groups (e.g. the mortality rates for people with lower and higher socioeconomic status over 40 years), this function computes the proportional change of absolute and relative inequalities, expressed in changes in rate differences and rate ratios. The function implements the algorithm proposed by Mackenbach et al. 2015.
References
Mackenbach JP, Martikainen P, Menvielle G, de Gelder R. 2015. The Arithmetic of Reducing Relative and Absolute Inequalities in Health: A Theoretical Analysis Illustrated with European Mortality Data. Journal of Epidemiology and Community Health 70(7): 730-36. doi:10.1136/jech-2015-207018
Examples
# This example reproduces Fig. 1 of Mackenbach et al. 2015, p.5
# 40 simulated time points, with an initial rate ratio of 2 and
# a rate difference of 100 (i.e. low status group starts with a
# prevalence rate of 200, the high status group with 100)
# annual decline of prevalence is 1% for the low, and 3% for the
# high status group
n <- 40
time <- seq(1, n, by = 1)
lo <- rep(200, times = n)
for (i in 2:n) lo[i] <- lo[i - 1] * .99
hi <- rep(100, times = n)
for (i in 2:n) hi[i] <- hi[i - 1] * .97
prev.data <- data.frame(lo, hi)
# print values
inequ_trend(prev.data, lo, hi)
#> $data
#> lo hi rr rd
#> 1 200.0000 100.00000 2.000000 100.0000
#> 2 198.0000 97.00000 2.041237 101.0000
#> 3 196.0200 94.09000 2.083324 101.9300
#> 4 194.0598 91.26730 2.126280 102.7925
#> 5 192.1192 88.52928 2.170120 103.5899
#> 6 190.1980 85.87340 2.214865 104.3246
#> 7 188.2960 83.29720 2.260533 104.9988
#> 8 186.4131 80.79828 2.307141 105.6148
#> 9 184.5489 78.37434 2.354711 106.1746
#> 10 182.7034 76.02311 2.403262 106.6803
#> 11 180.8764 73.74241 2.452814 107.1340
#> 12 179.0677 71.53014 2.503387 107.5375
#> 13 177.2770 69.38424 2.555004 107.8927
#> 14 175.5042 67.30271 2.607684 108.2015
#> 15 173.7492 65.28363 2.661451 108.4655
#> 16 172.0117 63.32512 2.716326 108.6866
#> 17 170.2916 61.42537 2.772333 108.8662
#> 18 168.5886 59.58260 2.829494 109.0060
#> 19 166.9028 57.79513 2.887834 109.1076
#> 20 165.2337 56.06127 2.947377 109.1725
#> 21 163.5814 54.37943 3.008148 109.2020
#> 22 161.9456 52.74805 3.070172 109.1975
#> 23 160.3261 51.16561 3.133474 109.1605
#> 24 158.7229 49.63064 3.198082 109.0922
#> 25 157.1356 48.14172 3.264022 108.9939
#> 26 155.5643 46.69747 3.331321 108.8668
#> 27 154.0086 45.29655 3.400008 108.7121
#> 28 152.4685 43.93765 3.470111 108.5309
#> 29 150.9439 42.61952 3.541660 108.3243
#> 30 149.4344 41.34093 3.614684 108.0935
#> 31 147.9401 40.10071 3.689214 107.8394
#> 32 146.4607 38.89769 3.765280 107.5630
#> 33 144.9961 37.73076 3.842915 107.2653
#> 34 143.5461 36.59883 3.922150 106.9473
#> 35 142.1106 35.50087 4.003019 106.6098
#> 36 140.6895 34.43584 4.085555 106.2537
#> 37 139.2826 33.40277 4.169794 105.8799
#> 38 137.8898 32.40068 4.255769 105.4891
#> 39 136.5109 31.42866 4.343517 105.0823
#> 40 135.1458 30.48580 4.433074 104.6600
#>
#> attr(,"class")
#> [1] "sj_inequ_trend"
# plot trends - here we see that the relative inequalities
# are increasing over time, while the absolute inequalities
# are first increasing as well, but later are decreasing
# (while rel. inequ. are still increasing)
plot(inequ_trend(prev.data, lo, hi))
#> Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
