R quantile fonksiyonu ile tüm gün hayretler oldum.

Abonelik sözleşmeleri nasıl çalışır, sezgisel bir kavram ve istatistikleri M. S. bir ben var, ama aman aman, bu belgeleri bana kafa karıştırıcı.

Gelen dokümanlar:

Q[i](p) = (1 - gamma) x[j] gama x[j 1],

Bunu dinliyorum. Bir tür içinbenquantile, x[j] ve x [j 1] arasında bir ilişkilendirme, gizemli bir sabit dayanıyorgama

1 <= <= 9, (j-m)/n < p < (j-m 1)/ n, x[j] Google sırasıdır istatistik, örnek büyüklüğü n, m. sürekli örnek tarafından belirlenir quantile yazın. Burada gama bağlıdır g kesirli bölümü = m-j np.

Nasıl j hesaplamak? m?

Sürekli örnek quantile için türleri (4 üzerinden 9), örnek doğrusal abonelik sözleşmeleri ile elde edilebilir aradeğerleme hazırlanmasında sırası arasında ve(k) istatistiği p:

p(k) = (k - alpha) (n - alfa - beta / 1), burada A ve b sabitleri belirlenir Türüne göre. Ayrıca, = (1 m alfa p - alfa - beta ve gama = g.

Şimdi gerçekten kayboldum. önce sabit olan p,, şimdi anlaşılan bir fonksiyonudur.

Yazın için 7 abonelik sözleşmeleri, varsayılan...

Yazın 7

p(k) = (k - 1) / (n - 1). Bu durumda, p(k) = mod[F(x[k])]. Bu S. tarafından kullanılır

Kimse bana yardım etmek ister misin? Özellikle p gösterimde bir fonksiyon olması kafamı karıştırdı ve bir sabit, bana nemve şimdi bazı özel j hesaplamak içinp.

Cevapları buraya göre daha iyi burada neler olduğunu açıklayan bazı revize belgeler sunması umuyorum.

quantile.R source code ya da yazın: quantile.varsayılan

Şaşırman anlaşılabilir bir durum. Bu belge korkunç bir şey. Kağıt geri dön dayanan zorunda kaldım (Hyndman, R. J.; Fan, Y. (Kasım 1996). "İstatistik Paket örnek abonelik sözleşmeleri".Amerikalı İstatistikçi50 (4): 361-365. Anlamak içindoi:10.2307/2684934) . Hadi ilk sorun ile başlar.

1 <= <= 9, (j-m)/n < p < (j-m 1)/ n, x[j], Google sipariş istatistiği, n, örnek boyutu, ve m sabiti tarafından belirlenir örnek quantile yazın. Burada gama m-j g = np kesirli kısmını bağlıdır.

İlk bölümü düz kağıt geliyor, ama belgelere yazarlar atlanmış ne j = int(pn m). Bu Q[i](p) sadece iki sipariş istatistikleri ile (sıralı) p kısmını olmaya en yakın gözlem bağlıdır anlamına gelir. (Terim aşina olanlar, benim gibi, "" gözlemler bir dizi, sıralanmış bir dizi.) sipariş istatistikleri

Ayrıca, bu son cümle çok yanlış. Bunu şöyle olmalıdır

Burada gama np kesirli kısmını bağlıdır m, m-j g = np

m gelince her şey ortada. m seçilen bağlıdır. Bu yüzden sadece Q[i] gibi quantile fonksiyonu m kabul edilmelidir m[i]. Ve 2, m algoritmalar 1 için 3, m/2 -1, ve diğerleri için, bir sonraki bölümü olan 0,.

Sürekli örnek için quantile türleri (4 üzerinden 9), örnek abonelik sözleşmeleri hazırlanmasında, sipariş istatistiği ve p(k) arasında doğrusal aradeğerleme ile elde edilebilir:

p(k) = (k - alpha) A ve b / (n - alfa - beta 1), sabitler türü tarafından belirlenir. Ayrıca, m = alfa p(1 - alfa - beta ve gama = g.

Bu gerçekten kafa karıştırıcı. Belgeleri p(k) dediği önce p aynı şey değildir. p(k) plotting position. Gazetede, yazarlar p olarak yazkyardımcı olur., Özellikle m, p için ifade orjinali ** 17 m = alpha p * (1 - alpha - beta). Kavramsal olarak, algoritmalar için 4-9 puan (pk, x[k]) çözüm (, *p*23)almak için aralıklı vardır. Her algoritma sadece p algoritması farklıdırk.

Son olarak, sadece R Ler kullandığını belirten.

Orijinal kağıt "örnek bir quantile için istenen Özellikleri 1 tarafından tüm karşılayan." işlev ve #8 tercih Devletleri 6 bir listesini verir Hepsini karşılar, ama başka alanlarda da hoşlanmıyorlar #5 (ilkelerinden türetilmiş daha fenomenolojik). Non-stat abonelik sözleşmeleri düşünün ve wikipedia açıklanmıştır kendim gibi inekler nedir #2.

BTW, dreeves answer, * yanıt olarak önemli şeyleri farklı yapar. Eşleme anladım. * Kolay anlamak için ise, (a) daha kolay saçma parametreleri ile kendini ayağından vurmak ve (b) R algoritması #2 yapamaz. (Burada *dört parametre açısından diğer tüm algoritmaları #2, ama basit bir genelleme verir yapamaz Birleşik Devletleri* 30*.)