10. The following R code is an R function to simulate realisations from a distribution with...

 

Transcribed Image Text:

10. The following R code is an R function to simulate realisations from a distribution with density function (x) x h(x) = 1/(1+x). The method used by the R code is the ratio-of-uniform algorithm. ru.sim=function(n) { r = NULL for (i in 1:n) { t = 1 while (t > 0) { } u = runif(1, 0, 1) v = runif(1, -1, 1) t = u + v - 1 r[i] = u/v r (a) Explain why the while loop in the above ratio of uniform algorithm simulates (u, v) from Ch = {(u, v) 0 10. The following R code is an R function to simulate realisations from a distribution with density function (x) x h(x) = 1/(1+x). The method used by the R code is the ratio-of-uniform algorithm. ru.sim=function(n) { r = NULL for (i in 1:n) { t = 1 while (t > 0) { } u = runif(1, 0, 1) v = runif(1, -1, 1) t = u + v - 1 r[i] = u/v r (a) Explain why the while loop in the above ratio of uniform algorithm simulates (u, v) from Ch = {(u, v) 0