Fix outlier detection

[stevemussmann]

Addresses issue 37: #37

Briefly, run_bgc.sh now runs estpost to extract tau-alpha (ta) and tau-beta (tb) values from .hdf5 files.

combine_bgc_output.R now parses the ta and tb files.

get_bgc_outliers.R calculates mean values of ta and tb, and uses them to assign outliers. The square roots of the reciprocals of ta and tb are used for this purpose in the qnorm() function (e.g., qnorm(qn_lower,0,sqrt(1/ta)))

plot_traces.R now produces trace plots for ta and tb.

[codacy-production]

Not up to standards ⛔

🔴 Issues 2 medium

Alerts:
⚠ 2 issues (≤ 0 issues of at least minor severity)

Results:
2 new issues

Category	Results
BestPractice	2 medium

View in Codacy

_{TIP This summary will be updated as you push new changes. Give us feedback}