https://wwwn.cdc.gov/nchs/nhanes/default.aspx

• Scrape .DOC files to Pandas DataFrame
• Parse .XPT and mortality .DAT files and convert to Pandas DataFrame
• Parse accelerometry .XPT files for 2003-2006 and 2011-2014 surveys to NumPy arrays

pip install pynhanes

NHANES website has hierarchical organization of data:

• Surveys (e.g. "2011-2012") ->

  • Components (e.g. "Questionnaire") ->

    • Categories (e.g. "Occupation") ->

      • Data variables (.DOC and .XPT files)

It is conveninet to have all data in Pandas DataFrame of NumPy arrays for data analysis. This repo is here to help you make it.

NOTE: Please, keep in mind, that some NHANES data fields have been recoded since 1999. Make sure you have reviewed the NHANES website and understand how the code processed and changed the data. Especially pay attention to categorical data. This may have effect on data analysis results.

NHANES Parser lib offers tool to get data in Pandas and NumPy:

1. Create a working folder, e.g. ~/work/NHANES/, copy notebooks from the repository folder sripts to the working folder and create subfolders XPT, CSV, NPZ

2. Copy nhanes_variables.json from the repository folder sripts to your CSV subfolder

3. Run parse_codebook.ipynb to scrape hierarchical structure of NHANES website to Pandas DataFrame (saves data to CSV subfolder)

4. Use pywgetxpt to download needed .XPT category files for all survey years (pywgetxpt DEMO -o XPT saves DEMO data to XPT subfolder)

5. Run parse_userdata.ipynb to get a list of selected data variable fields and converts .XPT and mortality .DAT files to Pandas DataFrame (saves data to CSV subfolder)

6. Optionally run parse_activity.ipynb to convert 2003-2006 and 2011-2014 accelerometry data to NumPy arrays (saves data in NPZ subfolder)

7. Run load_and_plot.ipynb to see an example of how to load and hadle parsed data

* parse_codebook.ipynb produces a codebook DataFrame which is a handy tool to convert numerically-encoded values to human-readable labels

** parse_userdata.ipynb may combine several variables into a sinle variable. Normally you would like to do that if:

a) Same data field has alternative names in diffrenet survey years (but be careful since the range of values may have changed -see the codebook):

SMD090, SMD650 - Avg # cigarettes/day during past 30 days

b) It is more reasonable to treat data fields together:

SMQ020, SMQ120, SMQ150 - Smoked at least 100 cigarettes in life / a pipe / cigars at least 20 times in life