Publications by xi'an
composition versus inversion
While trying to convey to an OP on X validated why the inversion method was not always the panacea in pseudo-random generation, I took the example of a mixture of K exponential distributions when K is very large, in order to impress (?) upon said OP that solving F(x)=u for such a closed-form cdf F was very costly even when using a state-of-the-ar...
1621 sym R (263 sym/3 pcs) 2 img
plusquamperfect squares
A perfect riddle: For some perfect squares, when you remove the last digit, you get another perfect square. The first five perfect squares are 16, 49, 169, 256 and 361. What are the next three ones? Is there a more than perfect square other than 169 such that removing the last two digits returns a perfect square? Writing an R code for plusquamp...
893 sym R (244 sym/1 pcs)
ABC in Svalbard [#1]
It started a bit awkwardly for me as I ran late, having accidentally switched to UK time the previous evening (despite a record-breaking biking-time to the University!), then the welcome desk could not find the key to the webinar room and I ended up following the first session from my office, by myself (and my teapot)… Until we managed to reuni...
3414 sym 2 img
ABC in Svalbard [#2]
The second day of the ABC wwworkshop got a better start than yesterday [for me] as I managed to bike to Dauphine early enough to watch the end of Gael’s talk and Matias Quiroz’ in full on the Australian side (of zoom). With an interesting take on using frequency-domain (pseudo-)likelihoods in complex models. Followed by two talks by David Fra...
4261 sym 4 img
the new DIYABC-RF
My friends and co-authors from Montpellier have released last month the third version of the DIYABC software, DIYABC-RF, which includes and promotes the use of random forests for parameter inference and model selection, in connection with Louis Raynal’s thesis. Intended as the earlier versions of DIYABC for population genetic applications. Bien...
1897 sym 4 img
ABC in Svalbard [the day after]
The following and very kind email was sent to me the day after the workshop thanks once again to make the conference possible. It was full of interesting studies within a friendly environment, I really enjoyed it. I think it is not easy to make a comfortable and inspiring conference in a remote version and across two continents, but this has bee...
2066 sym 2 img
simulating Maxwell distribution
A question that came out on X validated a few days ago is how to efficiently simulate from a distribution with density x²φ(x). (Obviously properly normalised since the second moment of the standard Normal distribution is one.) The first solution that came out (by Jarle Tufto) exploits the fact that this density corresponds to a signed root of...
1531 sym 2 img
ten computer codes that transformed science
In a “Feature” article of 21 January 2021, Nature goes over a poll on “software tools that have had a big impact on the world of science”. Among those, the Fortran compiler (1957), which is one of the first symbolic languages, developed by IBM. This is the first computer language I learned (in 1982) and one of the two (with SAS) I ever co...
2515 sym 2 img
a common confusion between sample and population moments
Related To leave a comment for the author, please follow the link and comment on their blog: R – Xi'an's Og. R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job. Want to share your content on R-bloggers? click here i...
397 sym 2 img
one-way random walks
A rather puzzling riddle from The Riddler on an 3×3 directed grid and the probability to get from the North-West to the South-East nodes following the arrows. Puzzling because while the solution could be reasonably computed with an R code like sucz=0 for(i in 1:2^12){ path=intToBits(i)[1:12] sol=0 for(j in 1:12)sol=max(sol, pr...
976 sym R (205 sym/1 pcs) 2 img