Kuantil regresyon için hangi tanı grafikleri var?


25

OLS ile ilgili sorumu takiben merak ediyorum: kuantil regresyon için hangi tanı grafikleri var? (ve bunların R uygulaması var mı?)

Hızlı bir google araması zaten (daha önce hiç duymadığım) solucan komplosuyla geldi ve hakkında bilmeniz gereken daha fazla yöntemi bilmekten mutluluk duyarım. (bunlardan herhangi biri kuantil regresyon için kullanılan OLS'tan biri mi?)


Sanırım, oyun kütüphanesinin solucan komplosunun bir uygulaması olduğunu buldunuz.
Peter Ellis,

Yanıtlar:


16

Kuantil regresyon, yanıt değişkeninin neredeyse sürekli olduğunu varsaymak dışında, dağılım varsayımları, yani artıklar hakkındaki varsayımlarda bulunmaz. Eğer bir fonksiyon tahmincisi X olarak tek bir nicelik tahmin etme problemini ele alıyorsanız, uzak ve uzakta yanlış giden başlıca şeyler, lineer öngörücü yanlış belirlenmesidir.XβAltta yatan, yani doğrusal olmayan efektleri (ortak bir problem) veya etkileşim etkilerini içerememek. Önerilen en az iki yaklaşım var. İlk olarak, örnek büyüklüğünüz büyükse, sadece daha esnek bir model kullanın. İyi bir uzlaşma, kısıtlı kübik splinelar (doğal splinelar) gibi regresyon spline'ları kullanarak tüm ana efektlerin lineer olmamasını sağlamaktır. O zaman etkileşimler dışında kontrol edilmesi gereken hiçbir şey yoktur. İkinci yaklaşım, modelin basit olduğunu (neden?) Ancak karmaşık olmasına izin vermek, sonra da karmaşık eklemelerin basit modele etkisini değerlendirmektir. Örneğin, doğrusal olmayan veya etkileşim terimlerinin birleşik katkılarını veya her ikisini de değerlendirebiliriz. Bir örnek, R rmsvequantregpaketler. Parametre sayısını sınırlamak için bir uzlaşma etkileşimi formu kullanılır. Etkileşimler iki kat lineer olmamakla sınırlıdır.

require(rms)
# Estimate 25th percentile of y as a function of x1 and x2
f <- Rq(y ~ rcs(x1, 4) + rcs(x2, 4) + rcs(x1, 4) %ia% rcs(x2, 4), tau=.25)
# rcs = restricted cubic spline, here with 4 default knots
# %ia% = restricted interaction
# To use general interactions (all cross product terms), use:
# f <- Rq(y ~ rcs(x1, 4)*rcs(x2, 4), tau=.25)
anova(f)   # get automatic combined 'chunk' tests: nonlinearity, interaction
# anova also provides the combined test of complexity (nonlin. + interact.)
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.