# \[ANN\] blog post benchmarking julia vs R Gaussian mixture models

**URL:** <https://discourse.julialang.org/t/ann-blog-post-benchmarking-julia-vs-r-gaussian-mixture-models/26071>\
**Category:** General Usage\
**Tags:** announcement\
**Created:** [July 6, 2019, 12:51pm UTC](https://discourse.julialang.org/t/ann-blog-post-benchmarking-julia-vs-r-gaussian-mixture-models/26071 "2019-07-06T12:51:36Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![floswald](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/floswald/32/195_2.png) [@floswald](https://discourse.julialang.org/u/floswald)\
**Post date:** [July 6, 2019, 12:51pm UTC](https://discourse.julialang.org/t/ann-blog-post-benchmarking-julia-vs-r-gaussian-mixture-models/26071/1 "2019-07-06T12:51:36Z")

</div>

hi all!

I wrote up a [quick blog post](https://floswald.github.io/post/em-benchmarks/) on my first experiences with the [GaussianMixtures.jl](https://github.com/davidavdav/GaussianMixtures.jl) package by @davidavdav. My main alternative would be an R implementation, hence the benchmark is comparing 2 R with 2 julia implementations. Julia outperforms R (surprise? 😄). Comments most welcome!

---

<div class="post-metadata">

**Author:** ![mkborregaard](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mkborregaard/32/556_2.png) [@mkborregaard](https://discourse.julialang.org/u/mkborregaard)\
**Post date:** [July 6, 2019, 9:39pm UTC](https://discourse.julialang.org/t/ann-blog-post-benchmarking-julia-vs-r-gaussian-mixture-models/26071/2 "2019-07-06T21:39:45Z")

</div>

wow, the R code is a monster compared to the julia bit

---

<div class="post-metadata">

**Author:** ![floswald](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/floswald/32/195_2.png) [@floswald](https://discourse.julialang.org/u/floswald)\
**Post date:** [July 7, 2019, 6:49am UTC](https://discourse.julialang.org/t/ann-blog-post-benchmarking-julia-vs-r-gaussian-mixture-models/26071/3 "2019-07-07T06:49:09Z")

</div>

yes. well, I couldn’t claim that’s the best possible R code, but I couldn’t figure out how to do the `repeat` without this `spread` function. There may be better ways. Also, to be fair, I had another 3 lines of not shown julia code defining the `logsumexp` function 😉