JW Gillard द्वारा वर्णित कई प्रकार की संभावनाएँ हैं जो दोनों वैरिएबल्स में त्रुटियों के साथ रैखिक प्रतिगमन के एक ऐतिहासिक अवलोकन में वर्णित हैं
यदि आप किसी अन्य पर एक विधि को चुनने के लिए विवरण या कारणों में कोई दिलचस्पी नहीं कर रहे हैं, सिर्फ आसान, केन्द्रक के माध्यम से लाइन आकर्षित करने के लिए है, जिसके साथ जाना ढलान के साथ , अर्थात मनाया गया मानक विचलन का अनुपात (ढलान के संकेत को और के सहसंयोजक के संकेत के समान ); जैसा कि आप शायद बाहर काम कर सकते हैं, यह के -axis पर एक अवरोधन देता है(x¯,y¯)β^=sy/sxxyyα^=y¯−β^x¯.
इस विशेष दृष्टिकोण के गुण हैं
- यह एक ही पंक्ति की तुलना देता है के खिलाफ के रूप में के खिलाफ ,xyyx
- यह पैमाने पर अपरिवर्तनीय है, इसलिए आपको इकाइयों के बारे में चिंता करने की आवश्यकता नहीं है,
- यह दो साधारण रेखीय प्रतिगमन रेखाओं के बीच स्थित है
- यह उन्हें पार करता है जहां वे टिप्पणियों के केंद्रक पर एक दूसरे को पार करते हैं, और
- इसकी गणना करना बहुत आसान है।
ढलान दो साधारण रैखिक प्रतिगमन ढलानों के ढलान का ज्यामितीय माध्य है। यह भी है कि यदि आप और अवलोकनों का मानकीकरण करते हैं, तो आपको 45 ° (या यदि नकारात्मक सहसंबंध है तो 135 °) पर एक रेखा खींची जाए और फिर लाइन को डी-मानकीकृत कर दिया जाए। इसे एक निहित धारणा बनाने के बराबर के रूप में भी देखा जा सकता है कि त्रुटियों के दो सेटों के भिन्न रूप, टिप्पणियों के दो सेटों के संस्करण के समानुपाती हैं; जहां तक मैं बता सकता हूं, आप दावा करते हैं कि यह नहीं पता कि यह कौन सा तरीका गलत है।xy
यह बताने के लिए कुछ R कोड है: चार्ट में लाल रेखा पर का OLS प्रतिगमन है, नीली रेखा पर का OLS प्रतिगमन है , और हरी रेखा यह सरल विधि है। ध्यान दें कि ढलान लगभग 5 होना चाहिए।YXXY
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")