<p><div><div><div><div data-width="100" data-show-faces="false" data-send="true" data-layout="button_count" data-height="21" data-href="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/"></div></div><div><a data-url="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/" data-counturl="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/" data-count="horizontal" data-text="astroABC: ABC SMC sampler for cosmological parameter estimation" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
(This article was first published on <strong><a href="https://xianblog.wordpress.com/2016/09/06/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/"> R – Xi'an's Og</a></strong>, and kindly contributed to <a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers)</a> </div>
<blockquote><p><span><em>“…the chosen statistic needs to be a so-called sufficient statistic in that any information about the parameter of interest which is contained in the data, is also contained in the summary statistic.”</em></span></p>
</blockquote>
<p><strong>E</strong>lise Jenningsa and Maeve Madigan <a target="_blank" href="http://arxiv.org/abs/1608.07606">arXived a paper</a> on a new Python code they developed for implementing ABC-SMC, towards astronomy or rather cosmology applications. They stress the parallelisation abilities of their approach which leads to “crucial speed enhancement” against the available competitors, <a target="_blank" href="https://xianblog.wordpress.com/2015/05/04/abc-and-cosmology/">abcpmc</a> and <a target="_blank" href="https://xianblog.wordpress.com/2015/05/04/abc-and-cosmology/">cosmoabc</a>. The version of ABC implemented there is “our” <a target="_blank" href="https://xianblog.wordpress.com/2009/02/10/abc-pmc-re-revised/">ABC PMC</a> where particle clouds are shifted according to mixtures of random walks, based on each and every point of the current cloud, with a scale equal to twice the estimated posterior variance. (The <a target="_blank" href="http://arxiv.org/abs/1608.07606">paper</a> curiously refers to non-astronomy papers through their arXiv version, even when they have been published. Like our <a target="_blank" href="https://xianblog.wordpress.com/2009/03/30/incredibly-ugly-squalid-pictures/">2008 Biometrika</a> paper.) A large part of the <a target="_blank" href="http://arxiv.org/abs/1608.07606">paper</a> is dedicated to computing aspects that escape me, like the constant reference to MPIs. The algorithm is partly automated, except for the choice of the summary statistics and of the distance. The tolerance is chosen as a (large) quantile of the previous set of simulated distances. <em>Getting comments from the designers of abcpmc and cosmoabc would be great.</em></p>
<blockquote><p><em><span>“It is clear that the simple Gaussian Likelihood assumption in this case, which neglects the effects of systematics yields biased cosmological constraints.”</span></em></p>
</blockquote>
<p>The last part of the <a target="_blank" href="http://arxiv.org/abs/1608.07606">paper</a> compares ABC and MCMC on a supernova simulated dataset. Which is somewhat a dubious comparison since the model used for producing the data and running ABC is not the same as the Gaussian version used with MCMC. Unsurprisingly, MCMC then misses the true value of the cosmological parameters and most likely and more importantly the true posterior HPD region. While ABC SMC (or PMC) proceeds to a concentration around the genuine parameter values. (There is no additional demonstration of how accelerated the approach is.)</p>
<p>Filed under: <a target="_blank" href="https://xianblog.wordpress.com/category/books/">Books</a>, <a target="_blank" href="https://xianblog.wordpress.com/category/statistics/r-statistics/">R</a>, <a target="_blank" href="https://xianblog.wordpress.com/category/statistics/">Statistics</a>, <a target="_blank" href="https://xianblog.wordpress.com/category/university-life/">University life</a> Tagged: <a target="_blank" href="https://xianblog.wordpress.com/tag/abc/">ABC</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/abc-pmc/">ABC-PMC</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/abc-smc/">ABC-SMC</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/abcpmc/">abcpmc</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/cosmoabc/">cosmoabc</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/cosmology/">cosmology</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/cosmopmc/">cosmoPMC</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/mcmc-algorithms/">MCMC algorithms</a>, <a target="_blank" href="https://xianblog.wordpress.com/tag/supernova/">supernova</a> <a target="_blank" rel="nofollow" href="http://feeds.wordpress.com/1.0/gocomments/xianblog.wordpress.com/31434/"><img alt="" src="http://feeds.wordpress.com/1.0/comments/xianblog.wordpress.com/31434/" border="0"></a> <img alt="" src="https://i1.wp.com/pixel.wp.com/b.gif?resize=1%2C1&ssl=1" border="0" data-recalc-dims="1"></p>
<p><div><div><div><div data-width="55" data-show-faces="false" data-send="true" data-layout="box_count" data-height="62" data-href="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/"></div></div><div><a data-url="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/" data-counturl="https://www.r-bloggers.com/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/" data-count="vertical" data-text="astroABC: ABC SMC sampler for cosmological parameter estimation" class="twitter-share-button" data-via="rbloggers" href="https://twitter.com/share"></a></div><div></div></div><div></div></div>
<div>
<div>To <strong>leave a comment</strong> for the author, please follow the link and comment on their blog: <strong><a href="https://xianblog.wordpress.com/2016/09/06/astroabc-abc-smc-sampler-for-cosmological-parameter-estimation/"> R – Xi'an's Og</a></strong>.</div>
<hr>
<a rel="nofollow" href="https://www.r-bloggers.com/">R-bloggers.com</a> offers <strong><a rel="nofollow" href="https://feedburner.google.com/fb/a/mailverify?uri=RBloggers">daily e-mail updates</a></strong> about <a title="The R Project for Statistical Computing" rel="nofollow" href="https://www.r-project.org/">R</a> news and <a title="R tutorials" rel="nofollow" href="https://www.r-bloggers.com/?s=tutorial">tutorials</a> on topics such as: <a title="Data science" rel="nofollow" href="https://www.r-bloggers.com/?s=data%20science">Data science</a>, <a title="Big Data" rel="nofollow" href="https://www.r-bloggers.com/?s=Big%20Data">Big Data, <a title="R jobs" rel="nofollow" href="https://www.r-users.com/">R jobs</a>, visualization (<a title="ggplot and ggplot2 tutorials" rel="nofollow" href="https://www.r-bloggers.com/?s=ggplot2">ggplot2</a>, <a title="Boxplots using lattice and ggplot2 tutorials" rel="nofollow" href="https://www.r-bloggers.com/?s=boxplot">Boxplots</a>, <a title="Maps and gis" rel="nofollow" href="https://www.r-bloggers.com/?s=map">maps</a>, <a title="Animation in R" rel="nofollow" href="https://www.r-bloggers.com/?s=animation">animation</a>), programming (<a title="RStudio IDE for R" rel="nofollow" href="https://www.r-bloggers.com/?s=RStudio">RStudio</a>, <a title="Sweave and literate programming" rel="nofollow" href="https://www.r-bloggers.com/?s=sweave">Sweave</a>, <a title="LaTeX in R" rel="nofollow" href="https://www.r-bloggers.com/?s=LaTeX">LaTeX</a>, <a title="SQL and databases" rel="nofollow" href="https://www.r-bloggers.com/?s=SQL">SQL</a>, <a title="Eclipse IDE for R" rel="nofollow" href="https://www.r-bloggers.com/?s=eclipse">Eclipse</a>, <a title="git and github, Version Control System" rel="nofollow" href="https://www.r-bloggers.com/?s=git">git</a>, <a title="Large data in R using Hadoop" rel="nofollow" href="https://www.r-bloggers.com/?s=hadoop">hadoop</a>, <a title="Web Scraping of google, facebook, yahoo, twitter and more using R" rel="nofollow" href="https://www.r-bloggers.com/?s=Web+Scraping">Web Scraping</a>) statistics (<a title="Regressions and ANOVA analysis tutorials" rel="nofollow" href="https://www.r-bloggers.com/?s=regression">regression</a>, <a title="principal component analysis tutorial" rel="nofollow" href="https://www.r-bloggers.com/?s=PCA">PCA</a>, <a title="Time series" rel="nofollow" href="https://www.r-bloggers.com/?s=time+series">time series</a>, <a title="finance trading" rel="nofollow" href="https://www.r-bloggers.com/?s=trading">trading</a>) and more...
<img height="1" alt="" width="1" src="http://feeds.feedburner.com/~r/RBloggers/~4/gGRkfshEdco"></a></div>