Publications by xi'an
the Kelly criterion and then some
The Kelly criterion is a way to optimise an unlimited sequence of bets under the following circumstances: a probability p of winning each bet, a loss of a fraction a of the sum bet, a gain of a fraction b of the sum bet, and a fraction f of the current fortune as the sum bet. Then is the fraction optimising the growth Here is a rendering of the...
1850 sym 8 img
bayess’ back! [on CRAN]
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...
396 sym 2 img
optimal leap year
A riddle about leap years: a solar year consists of approximately 365.24217 mean solar days, which is why there is a leap year approximately every four years. Approximately because the Gregorian calendar plans 97 and not 100 leap years over 400 years. Is this the optimal solution? No, since the Gregorian difference is 3.3 10⁻⁴ day per year, o...
1000 sym R (83 sym/1 pcs) 2 img
Bayes Factors for Forensic Decision Analyses with R [book review]
My friend EJ Wagenmaker pointed me towards an entire book on the BF by Bozza (from Ca’Foscari, Venezia), Taroni and Biederman. It is providing a sort of blueprint for using Bayes factors in forensics for both investigative and evaluative purposes. With R code and free access. I am of course unable to judge of the relevance of the approach for f...
4968 sym 2 img
shelled and riddled
Consider a shell game with three shells and a ball where only the location of the shell with the ball is exchanged with the location of an empty shell, randomly chosen. If one starts with the ball as rightmost, what is the distribution of the location of the ball after N steps? Running an exploratory R code like o=rep(0,3) for(n in 1:1e6){ b=...
935 sym R (113 sym/1 pcs) 2 img
optimal Gaussian zorbing
A zorbing puzzle from the Riddler: cover the plane with four non-intersecting disks of radius one towards getting the highest probability (under the standard bivariate Normal distribution). As I could not see a simple connection between the disks and the standard Normal, beyond the probability of a disk being given by a non-central chi-square cdf...
1443 sym R (602 sym/1 pcs) 4 img
fast track & slow lane
A riddle from the Riddler where N cars going at random (iid) speeds drive a road with a slow and a fast lane, each car choosing the fast lane iff any of the cars ahead in the slow lane is slowerthan them. With the question of the average number of car convoys. If there were one single lane, the problem would be to determine how many times a small...
1713 sym R (323 sym/1 pcs) 4 img
why is this algorithm simulating a Normal variate?
A backward question from X validated as to why the above is a valid Normal generator based on exponential generations. Which can be found in most textbooks (if not ours). And in The Bible, albeit as an exercise. The validation proceeds from the (standard) Exponential density dominating the (standard) Normal density and, according to Devroye, may ...
1613 sym R (456 sym/2 pcs) 2 img
riddle of the week
The Riddler of April 1 offered this simple question: start with the number 1 and then try to reach a target number through a series of steps. For each step, you can always choose to double the number you currently have. However, if the number happens to be one (1) more than an odd multiple of 3, you can choose to “reduce” — that is, subtrac...
1086 sym R (81 sym/1 pcs)
null recurrent = zero utility?
The stability result that the ratio converges holds for a Harris π-null-recurrent Markov chain for all functions f,g in L¹(π) [Meyn & Tweedie, 1993, Theorem 17.3.2] is rather fascinating. However, it is unclear it can be useful in simulation environments, as for the integral priors we have been studying over the years with Juan Antonio Cano a...
1585 sym 6 img