Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Coding the experiment is straightforward:
g=sample(1:N,N,rep=TRUE) o=0*g for(i in 1:N){ if(min(o[g[i]:N])){f=f+1;break() }else{ o[min(which(!o[g[i]:N]))+g[i]-1]=1 }}}
returning an estimated probability of approximately 0.764.
As I had some free time during the early mornings at ISBA 2022, I tried to reformulate the question as a continuous event on uniform order statistics, turning to be at most one uniform larger than (N-1)/N, at most two larger than (N-2)/N, and so on… Asking the question on math.stackexchange quickly produced an answer that reversed engineered my formulation back to the goats (or parking lot), with a generic probability of
which of course coincides with the Monte Carlo approximation!
As an aside, I once drank South-African wines named Goats-do-Roam and Goat-Roti at my friends Jim and Maria’s place, and they were quite enjoyable!
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 if you have a blog, or here if you don't.