Qqline () tarafından R'de üretilen hattın kullanımı nedir?


20

qqnorm()R fonksiyonu normal QQ-arsa üreten ve qqline()birinci ve üçüncü çeyrek içinden geçen bir çizgi ekler. Bu hattın kökeni nedir? Normalliği kontrol etmek yardımcı olur mu? Bu klasik çizgi değildir ( muhtemelen doğrusal ölçeklemeden sonra diyagonal ).y=x

İşte bir örnek. İlk önce ampirik dağılım fonksiyonunu teorik dağılım fonksiyonu ile karşılaştırıyorum : Şimdi qq-plot'u ; bu grafik kabaca bir önceki grafiğin (doğrusal olmayan) ölçeklenmesine karşılık gelir: Ama burada R qqline ile qq-plot: Bu son grafik ilk grafikte olduğu gibi kalkış göstermiyor.y = μ +N-(μ^,σ^2)birikimli dağılım fonksiyonlarının karşılaştırılmasıy=μ^+σ^xqqnorm "iyi" çizgi ile birlikteqqnorm ve qqline

Yanıtlar:


9

Resimde görebileceğiniz gibi,resim açıklamasını buraya girin

tarafından edinilmiş

> y <- rnorm(2000)*4-4
> qqnorm(y); qqline(y, col = 2,lwd=2,lty=2)

köşegen anlamlı olmaz çünkü ilk eksen bir dağılımının teorik nicelikleri açısından ölçeklendirilir . Çizgiyi ayarlamak için birinci ve üçüncü çeyrekleri kullanmanın, ampirik ortalama ve varyans kullanımı ile karşılaştırıldığında, normal dağılımın parametrelerini tahmin etmek için sağlam bir yaklaşım verdiğini düşünüyorum. Hattan sapmalar (kuyruklar hariç) normallik eksikliğinin bir göstergesidir.N(0,1)


"Doğrusal ölçeklemeden sonra" köşegen burada abline (ortalama (y), sd (y)) ile elde edilir. Burada normal verileri simüle edersiniz, bu nedenle bu iki çizgi yakındır. Ancak bazen veriler normal dağılıma yakın değildir ancak qqplot qqline'a yakındır, ancak "ölçeklemeden sonra" diyagonaline yakın değildir.
Stéphane Laurent

... soruma bir örnek ekleyeceğim
Stéphane Laurent

4
Sanırım bu, çeyrekleri kullanmanın ampirik ortalama ve varyans kullanmaktan daha sağlam olduğunu belirtmekti.
Xi'an

1
Tamam çok teşekkür ederim. Şimdi bu açık görünüyor. Qqline tercih edilebilir çünkü bazen pratikte kuyruklardaki normallik kabul edilebilir. Ancak qqline'ı çizmek için gerçek bir ihtiyaç yoktur: görsel bir kontrol yeterlidir - ihtiyacımız olan tek şey QQ-grafiğini anlamaktır :)
Stéphane Laurent

1
Tamam - etiketliyorum, ancak cevabın kendisi tatmin edici değildi: tartışmamızla birlikte cevap; ama bu benim hatam: örneği eklemeden önce sorum net değildi. Bu arada sorum KS testi ile biraz ilgili: ks.test (x, "pnorm", mu.hat, sigma.hat yazdığımızda ve tahminlerinin seçimi ne olacak ? )? σμ^σ^
Stéphane Laurent
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.