วิธีตัวอย่างจากการกระจายคันทอร์?

อะไรจะเป็นวิธีที่ดีที่สุดในการสุ่มตัวอย่างจากการกระจายคันทอร์ ? มันมี cdf เท่านั้นและเราไม่สามารถกลับด้านได้

ง่าย ๆ : ตัวอย่างจากการแจกแจงเครื่องแบบและบันทึกจากไบนารีเป็นสามส่วนตีความแต่ละ "1" เป็น "2" (นี่คือวิธีการแปลงความน่าจะเป็นแบบผกผัน: มันจะกลับ CDF!)$(0,1)$

นี่คือRการนำไปใช้งานซึ่งเขียนขึ้นในลักษณะที่ควรปรับพอร์ตให้พร้อมกับสภาพแวดล้อมการประมวลผลเกือบทุกชนิด

binary.to.ternary <- function(x) {
  y <- 0
  x <- round(2^52 * x)
  for (i in 1:52) {
    y <- y + 2*(x %% 2)
    y <- y/3
    x <- floor(x/2)
  }
  y
}

n <- 1000
x <- runif(n)
y <- binary.to.ternary(x)
plot(ecdf(y), pch=".")