rmarkello / pypmi

A very small toolbox for working with data from the Parkinson's Progression Markers Initiative (PPMI)
https://pypmi.readthedocs.io
BSD 3-Clause "New" or "Revised" License
14 stars 8 forks source link

PyPMI

This package provides a Python interface for working with data from the Parkinson's Progression Markers Initiative (PPMI).

Build Status Codecov Documentation Status License PyPi

Installation and setup

This package requires Python >= 3.6. If you have the correct version of Python installed, you can install this package by opening a terminal and running the following:

git clone https://github.com/rmarkello/pypmi.git
cd pypmi
python setup.py install

Some of the functionality of this package—primarily, the functions to work with neuroimaging data—also requires the use of Docker.

Overview

The PPMI is an ongoing longitudinal study that begin in early 2010 with the primary goal of identifying biomarkers of Parkinson's disease (PD) progression. To date, the PPMI has collected data from over 400 individuals with de novo PD and nearly 200 age-matched healthy participants, in addition to large cohorts of individuals genetically at-risk for PD. Data, made available on the PPMI website, include comphrensive clinical-behavioral assessments, biological assays, single-photon emission computed tomography (SPECT) images, and magnetic resonance imaging (MRI) scans.

While accessing this data is straightforward (researchers must simply sign a data usage agreement and provide information on the purpose of their research), the sheer amount of data made available can be quite overwhelming to work with. Thus, the primary goal of this package is to provide a Python interface to making working with the data provided by the PPMI easier. While this project is still very much under development, it is neverthless functional! Most useful may be the functions accesible by directly importing pypmi, which help wrangle the litany of raw CSV files provided by the PPMI, and in pypmi.bids, which helps convert raw neuroimaging data from the PPMI into BIDS format.

I hope to continue adding useful features to this package as I keep working with the data, but take a look below at development and getting involved if you're interested in contributing, yourself!

Usage

Once you have access to the PPMI database, log in to the database and follow these instructions:

  1. Select "Download" from the navigation bar at the top
  2. Select "Study Data" from the options that appear in the navigation bar
  3. Select "ALL" at the bottom of the left-hand navigation bar on the new page
  4. Click "Select ALL tabular data (csv) format" and then press "Download>>" in the top right hand corner of the page
  5. Unzip the downloaded directory and save it somewhere on your computer

Alternatively, you can use the pypmi API to download the data programatically:

>>> import pypmi
>>> pypmi.fetch_studydata('all', user='username', password='password')

By default, the data will be downloaded to your current directory making it easy to load them in the future, but you can optionally provide a path argument to pypmi.fetch_studydata() to specify where you would like the data to go. (Alternatively, you can set an environmental variable $PPMI_PATH to specify where they should be downloaded to; this will be used if path is not specified.)

Once you have the data downloaded you can load some of it as a tidy pandas.DataFrame:

>>> behavior = pypmi.load_behavior()
>>> behavior.columns
Index(['participant', 'visit', 'date', 'benton', 'epworth', 'gds',
       'hvlt_recall', 'hvlt_recognition', 'hvlt_retention', 'lns', 'moca',
       'pigd', 'quip', 'rbd', 'scopa_aut', 'se_adl', 'semantic_fluency',
       'stai_state', 'stai_trait', 'symbol_digit', 'systolic_bp_drop',
       'tremor', 'updrs_i', 'updrs_ii', 'updrs_iii', 'updrs_iii_a', 'updrs_iv',
       'upsit'],
      dtype='object')

The call to pypmi.load_behavior() may take a few seconds to run—there's a lot of data to import!

If we want to query the data with regards to diagnosis it might be useful to load in some demographic information:

>>> demographics = pypmi.load_demographics()
>>> demographics.columns
Index(['participant', 'diagnosis', 'date_birth', 'date_diagnosis',
       'date_enroll', 'status', 'family_history', 'age', 'gender', 'race',
       'site', 'handedness', 'education'],
      dtype='object')

Now we can perform some interesting queries! As an example, let's just ask how many individuals with Parkinson's disease have a baseline UPDRS III score. We'll have to use information from both data frames to answer the question:

>>> import pandas as pd
>>> updrs = (behavior.query('visit == "BL" & ~updrs_iii.isna()')
...                  .get(['participant', 'updrs_iii']))
>>> parkinsons = demographics.query('diagnosis == "pd"').get('participant')
>>> len(pd.merge(parkinsons, updrs, on='participant'))
423

And the same for healthy individuals:

>>> healthy = demographics.query('diagnosis == "hc"').get('participant')
>>> len(pd.merge(healthy, updrs))
195

There's a lot of power gained in leveraging the pandas DataFrame objects, so take a look at the pandas documentation to see what more you can do! Also be sure to check out the pypmi documentation for more information.

Development and getting involved

This package has largely been developed in the spare time of a single graduate student (@rmarkello). While it would be :sparkles: amazing :sparkles: if anyone else finds it helpful, given the limited time constraints of graduate school this package is not currently accepting requests for new features.

However, if you're interested in getting involved in the project, we're thrilled to welcome new contributors! You should start by reading our contributing guidelines and code of conduct. Once you're done with that, take a look at our issues to see if there's anything you might like to work on. Alternatively, if you've found a bug, are experiencing a problem, or have a question, create a new issue with some information about it!

Acknowledgments

This package relies on data that can be obtained from the Parkinson's Progression Markers Initiative (PPMI) database https://www.ppmi-info.org/data. For up-to-date information on the study, visit https://www.ppmi-info.org/

The PPMI—a public-private partnership—is funded by the Michael J. Fox Foundation for Parkinson’s Research and funding partners, including AbbVie, Avid Radiopharmaceuticals, Biogen, BioLegend, Bristol-Myers Squibb, GE Healthcare, Genentech, GlaxoSmithKline (GSK), Eli Lilly and Company, Lundbeck, Merck, Meso Scale Discovery (MSD), Pfizer, Piramal Imaging, Roche, Sanofi Genzyme, Servier, Takeda, Teva, and UCB www.ppmi-info.org/fundingpartners.

License information

This codebase is licensed under the 3-clause BSD license. The full license can be found in the LICENSE file in the pypmi distribution.

All trademarks referenced herein are property of their respective holders.