Cantor dağıtımından nasıl örnek alınır?


19

Cantor dağıtımından örneklemenin en iyi yolu nedir ? Sadece cdf var ve tersine çeviremiyoruz.


4
Aslında, birisi matematik üzerine sordu: math.stackexchange.com/questions/1115907/…
RUser4512

İşte bazı ilginç takip soruları: standart sapma nedir? Moment üreten fonksiyon nedir? Tekdüze dağılımı için olan muadilleriyle nasıl kıyaslanıyorlar ? :-)(0,1)
whuber

5
Buraya bağlanan math.stackexchange yayınına başvurarak oluşturduğunuz sonsuz döngüyü seviyorum: p
Tasos Papastylianou

Yanıtlar:


23

Kolay: Tekdüze dağılımdan numune alın ve her "1" i "2" olarak yorumlayan ikili dosyadan üçlü dosyaya yeniden kodlayın. (Bu ters olasılık dönüşüm yaklaşımıdır: gerçekten CDF'yi tersine çevirir!)(0,1)

şekil

İşte Rhemen hemen her bilgi işlem ortamına kolayca taşınacak şekilde yazılmış bir uygulama.

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=".")

3
Bu yılın başlarında ben de biraz dolgun uygulamasını başlatan github.com/Henrygb/CantorDist.R fonksiyonları ile rCantor(), qCantor(), pCantor()ve daha az anlamlıdCantor()
Henry

1
@Henry Ne dcantoruygular? Tim'in belirttiği gibi, bu dağılımın yoğunluğu yoktur. Ayrık atomları da yoktur. Kesintisiz ama kesinlikle sürekli olmayan bir dağılımın arketip örneği. ( qcantor
BTW'nin uygulanmasını seviyorum

1
Gerçek dağılıma sadece sınırlı bir yaklaşımla uğraştığımızı unutmamalıyız. Diyelim ki 10 üçlü basamak hassas sayısı vardı (pratikte daha uzun olacaklar) ve basamakları daha da genişleyen bir değişkeni "temsil etmek" için 0.0222020002 ürettik. Aynı yorum, sürekli bir rv'ye sahip herhangi bir gerçek değerli rv için geçerli olmakla birlikte, tüm "temsil edilen" değerler, sonlu uzunluk yaklaşımının dayanabileceği "kümede" dir. Gerçek Cantor dağılımında, bu on haneli dizinin neredeyse tüm "devamları" kümede değildir. ...
ctd

1
@whuber Rastgele sayılar üretmenin her yönteminin ikinci cümlemde son derece hassas olduğunu açıkça kabul ettim. Tekrarlamayı seçtiğin ve verdiğiniz vurgu, orada gerçek noktamı kaçırdığınızı gösteriyor; sonlu hassasiyetten sürekli bir değişkeni temsil ettiğimde, böyle bir sonlu yaklaşımın temsil edebileceği gerçek değerler, üretmek istediğimiz "sette" dir. Bunun gibi bir değişkeni sonlu hassasiyete temsil ettiğimde, sonlu bir yaklaşımın temsil edebileceği gerçek değerlerin neredeyse hepsi kümede değildir. Oldukça farklı bir durum. ...
ctd

1
ctd ... yayınınıza hiçbir eleştiri gelmedi; okuyucuların göz ardı edebileceği ve özellikle Cantor setinin özelliklerini simüle ederek çıkarım yapmaya çalıştıkları takdirde düşünmek isteyebilecekleri bir noktaydı.
Glen_b -Monica
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.