changing the analysis script somewhat, factored out two more functions that take care of graphing and finding direct/indirect effects.  Changed the overall analysis script to "analysis_for_manuscript.R"