Hypertension definitions in app versus PMID 35616029

[bcjaeger]

Comparing the number of participants in our data versus Figure S1 of PMID 35616029, I see some minor differences. In the reproducible example below, I find we have the same number of participants overall but different numbers with hypertension.

library(cardioStatsUSA)
library(magrittr)
library(data.table)

data_app <- nhanes_data %>% 
 .[svy_subpop_htn == 1] %>% 
 .[svy_year %in% c( "2009-2010", 
                    "2011-2012", 
                    "2013-2014",
                    "2015-2016", 
                    "2017-2020")] %>% 
 .[demo_pregnant == 'No' | is.na(demo_pregnant)]

# Number of observations (unweighted); comparison to PMID = 35616029
nrow(data_app) # = 30,311 in Figure S1 
#> [1] 30311
nrow(data_app[htn_jnc7 == 'Yes']) # = 11,007 in Figure S1 
#> [1] 11017
nrow(data_app[htn_accaha == 'Yes']) # = 15,100 in Figure S1
#> [1] 15027

^{Created on 2022-11-04 with reprex v2.0.2}

These differences are small and do not prevent us from replicating results in the paper, but I think it would be good to note why they occur. @pmuntner, do you have any thoughts on this?