This project contains code for running an analysis of adverse drug events using the Multi-Item Gamma Poisson Shrinker (MGPS) model described in Empirical bayes screening for multi-item associations.
This analysis is designed to be small enough that you can run it on a single machine if you do not have access to a Hadoop cluster. You will need to have a version of CDH3 on your local machine, along with the version of Pig that is compatible with that version.
The input data for this analysis may be downloaded from the FDA's AERS website. You'll need to get the ASCII version of the data files for as many quarters as you would like to run over. For my own analysis, I used the data from 2008 through 2010.
The Pig scripts below assume that the input data is stored in three HDFS directories under the user's home directory: aers/drugs, aers/demos, and aers/reactions. All of the DRUG*.TXT files from the AERS website should go into aers/drugs, all of the DEMO*.TXT files should go into aers/demos, and all of the REAC*.TXT files should go into aers/reactions.
If you have not done so already, load the input data into the Hadoop cluster:
hdfs dfs -mkdir aers hdfs dfs -mkdir aers/drugs hdfs dfs -put DRUG*.TXT aers/drugs hdfs dfs -mkdir aers/demos hdfs dfs -put DEMO*.TXT aers/demos hdfs dfs -mkdir aers/reactions hdfs dfs -put REAC*.TXT aers/reactions
Each of these commands should be run from the project's top-level directory, i.e., the directory that contains this README file.
mvn package # Builds the Pig UDFs pig -f src/main/pig/step1_join_drugs_reactions.pig pig -f src/main/pig/step2_generate_drug_reaction_counts.pig pig -f src/main/pig/step3_generate_squashed_distribution.pig
At this point, you can optionally run the R code to solve the MGPS optimization problem. You will need to install the BB library in your local version of R using install.packages("BB") if you do not have it already.
hadoop fs -getmerge aers/drugs2_reacs_stats d2r_stats.csv Rscript src/main/R/ebgm.R d2r_stats.csv
The output from the optimization run may be plugged into the Pig script that scores the tuples, or you can just use the default parameters that are there now:
pig -f src/main/pig/step4_apply_ebgm.pig
The final output will be in aers/scored_drugs2_reacs. To generate the GEXF file of drug-drug interactions to load into Gephi, run:
hadoop fs -getmerge aers/scored_drugs2_reacs scored_d2r.csv ./src/main/python/gephi.py scored_d2r.csv > drugs.gexf