JW Gillard tarafından her iki Değişkende Hatalı Doğrusal Regresyonun Tarihsel Görünümünde açıklanan bir dizi olasılık vardır
Bir yöntemi başka bir yöntemle seçmekle ilgili ayrıntılar veya nedenlerle ilgilenmiyorsanız, çizgiyi centroid ile eğim , yani gözlemlenen standart sapmaların oranı (eğim işaretini ve kovaryansının işareti ile aynı yapmak ); muhtemelen çalışabileceğiniz gibi, bu 'nın eksenine kesişme sağlarβ = s y / s x x y y α = ˉ y - β ˉ X .(x¯,y¯)β^=sy/sxxyyα^=y¯−β^x¯.
Bu özel yaklaşımın esasları
- Bu karşılaştırma, aynı hat verir karşı olarak karşı ,y y xxyyx
- ölçek değişmez, bu yüzden birimler hakkında endişelenmenize gerek yok,
- iki sıradan lineer regresyon çizgisi arasında yer alır.
- gözlemlerin merkezinde, birbirlerini geçtikleri yerden geçer ve
- hesaplamak çok kolaydır.
Eğim, iki sıradan lineer regresyon eğiminin eğimlerinin geometrik ortalamasıdır. Ayrıca ve gözlemlerini standardize ederseniz , 45 ° 'de bir çizgi çizdiyseniz (veya negatif korelasyon varsa 135 °' de) ve daha sonra çizgiyi de standartlaştırırsanız elde edeceğiniz şey de budur . Ayrıca, iki hata kümesinin varyanslarının, iki gözlem kümesinin varyanslarıyla orantılı olduğu örtük bir varsayımın yapılmasına eşdeğer olarak görülebilir; anlayabildiğim kadarıyla bunun yanlış olduğunu bilmediğini iddia ediyorsun.yxy
Burada göstermek için bazı R kodudur: grafikte kırmızı çizgi OLS regresyonu ile , mavi çizgi OLS regresyonu ile ve yeşil çizgi, bu basit bir yöntemdir. Eğimin yaklaşık 5 olması gerektiğini unutmayın.X X YYXXY
X0 <- 1600:3600
Y0 <- 5*X0 + 700
X1 <- X0 + 400*rnorm(2001)
Y1 <- Y0 + 2000*rnorm(2001)
slopeOLSXY <- lm(Y1 ~ X1)$coefficients[2] #OLS slope of Y on X
slopeOLSYX <- 1/lm(X1 ~ Y1)$coefficients[2] #Inverse of OLS slope of X on Y
slopesimple <- sd(Y1)/sd(X1) *sign(cov(X1,Y1)) #Simple slope
c(slopeOLSXY, slopeOLSYX, slopesimple) #Show the three slopes
plot(Y1~X1)
abline(mean(Y1) - slopeOLSXY * mean(X1), slopeOLSXY, col="red")
abline(mean(Y1) - slopeOLSYX * mean(X1), slopeOLSYX, col="blue")
abline(mean(Y1) - slopesimple * mean(X1), slopesimple, col="green")