diff --git a/README.md b/README.md
index 421aa14..7147ff4 100644
--- a/README.md
+++ b/README.md
@@ -5,7 +5,7 @@
 BreakoutDetection is an open-source R package that makes breakout detection simple and fast. The BreakoutDetection package can be used in wide variety of contexts. For example, detecting breakout in user engagement post an A/B test, detecting [behavioral change](http://wiki.cbr.washington.edu/qerm/index.php/Behavioral_Change_Point_Analysis), or for problems in econometrics, financial engineering, political and social sciences.
 
 ## How the package works
-The underlying algorithm – referred to as E-Divisive with Medians (EDM) – employs energy statistics to detect divergence in mean. Note that EDM can also be used detect change in distribution in a given time series. EDM uses [robust statistical metrics](http://www.wiley.com/WileyCDA/WileyTitle/productCd-0470129905.html), viz., median, and estimates the statistical significance of a breakout through a permutation test. 
+The underlying algorithm – referred to as E-Divisive with Medians (EDM) – employs energy statistics to detect divergence in mean. Note that EDM can also be used detect change in distribution in a given time series. EDM uses [robust statistical metrics](http://www.wiley.com/WileyCDA/WileyTitle/productCd-0470129905.html), viz., median, and estimates the statistical significance of a breakout through a permutation test.
 
 In addition, EDM is non-parametric. This is important since the distribution of production data seldom (if at all) follows the commonly assumed normal distribution or any other widely accepted model. Our experience has been that time series often contain more than one breakout. To this end, the package can also be used to detect multiple breakouts in a given time series.
 
@@ -32,4 +32,16 @@ res = breakout(Scribe, min.size=24, method='multi', beta=.001, degree=1, plot=TR
 res$plot
 ```
 
-From the plot, we observe that the input time series experiences two breakouts and also has quite a few anomalies. The two red vertical lines denote the locations of the breakouts detected by the EDM algorithm. 
+The above yields the following plot:
+
+![](https://blog.twitter.com/content/dam/blog-twitter/archive/breakout_detectioninthewild96.thumb.1280.1280.png)
+
+From the plot, we observe that the input time series experiences two breakouts and also has quite a few anomalies. The two red vertical lines denote the locations of the breakouts detected by the EDM algorithm. Unlike the existing approaches mentioned earlier, EDM is robust in the presence of anomalies. The change in mean in the time series can be better viewed with the following annotated plot:
+
+![](https://blog.twitter.com/content/dam/blog-twitter/archive/breakout_detectioninthewild97.thumb.1280.1280.png)
+
+The horizontal lines in the annotated plot above correspond to the approximate (i.e., filtering out the effect of anomalies) mean for each window.
+
+## Acknowledgements
+
+We thank James Tsiamis and Scott Wong for their support and Nicholas James as the primary researcher behind this research.