กลับไปหน้าบทความ

อ่าน 6 นาที

การถดถอยของเดมิง

ในทางสถิติการถดถอยแบบเดมิง (Deming regression ) ซึ่งตั้งชื่อตามดับเบิลยู. เอ็ดเวิร์ด เดมิงเป็นแบบจำลองความคลาดเคลื่อนในตัวแปรที่พยายามหาเส้นตรงที่เหมาะสมที่สุด สำหรับ

การถดถอยของเดมิง

การถดถอยแบบเดมิง เส้นสีแดงแสดงข้อผิดพลาดทั้งในxและyซึ่งแตกต่างจากวิธีการกำลังสองน้อยที่สุดแบบดั้งเดิมที่วัดข้อผิดพลาดขนานกับ แกน yกรณีที่แสดงนี้ซึ่งวัดค่าเบี่ยงเบนในแนวตั้งฉาก เกิดขึ้นเมื่อข้อผิดพลาดในxและyมีค่าความแปรปรวนเท่ากัน

ในทางสถิติการถดถอยแบบเดมิง (Deming regression ) ซึ่งตั้งชื่อตามดับเบิลยู. เอ็ดเวิร์ด เดมิงเป็นแบบจำลองความคลาดเคลื่อนในตัวแปรที่พยายามหาเส้นตรงที่เหมาะสมที่สุด สำหรับ ชุดข้อมูลสองมิติมันแตกต่างจากการถดถอยเชิงเส้นแบบง่ายตรงที่มันคำนึงถึงความคลาดเคลื่อนในการสังเกตทั้งบนแกนxและ แกน yมันเป็นกรณีพิเศษของกำลังสองน้อยที่สุดโดยรวม (Total Least Squares ) ซึ่งอนุญาตให้มีตัวแปรทำนายได้หลายตัวและโครงสร้างความคลาดเคลื่อนที่ซับซ้อนกว่า

การถดถอยแบบเดมิงเทียบเท่ากับ การประมาณ ค่าความน่าจะเป็นสูงสุดของแบบจำลองข้อผิดพลาดในตัวแปรโดยที่ข้อผิดพลาดสำหรับตัวแปรทั้งสองถือว่าเป็นอิสระและมีการกระจายแบบปกติและอัตราส่วนของความแปรปรวนของตัวแปรทั้งสอง ซึ่งแสดงด้วยδ นั้น เป็นที่ทราบ[ 1 ]ในทางปฏิบัติ อัตราส่วนนี้อาจได้รับการประมาณจากแหล่งข้อมูลที่เกี่ยวข้อง อย่างไรก็ตาม ขั้นตอนการถดถอยไม่ได้คำนึงถึงข้อผิดพลาดที่อาจเกิดขึ้นในการประมาณอัตราส่วนนี้

การคำนวณการถดถอยแบบเดมิงนั้นยากกว่าการคำนวณการถดถอยเชิงเส้นแบบง่าย เพียงเล็กน้อย โปรแกรมซอฟต์แวร์ทางสถิติส่วนใหญ่ที่ใช้ในเคมีคลินิกมีฟังก์ชันการคำนวณการถดถอยแบบเดมิงให้ใช้งาน

แบบจำลองนี้ได้รับการแนะนำครั้งแรกโดยAdcock (1878)ซึ่งพิจารณากรณีδ  =  1 และต่อมาโดยKummell (1879) โดยทั่วไป ด้วยδ ใดๆ อย่างไรก็ตาม แนวคิดของพวกเขายังคงไม่ได้รับความสนใจมากนักเป็นเวลากว่า 50 ปี จนกระทั่งได้รับการฟื้นฟูโดยKoopmans (1936)และต่อมาได้รับการเผยแพร่มากขึ้นโดยDeming (1943)หนังสือเล่มหลังนี้ได้รับความนิยมอย่างมากใน สาขา เคมีคลินิกและสาขาที่เกี่ยวข้อง จนกระทั่งวิธีการนี้ถูกเรียกว่าการถดถอยของ Demingในสาขาเหล่านั้น[ 2 ]

ข้อกำหนด

สมมติว่าข้อมูลที่มีอยู่ ( y , x ) เป็นค่าสังเกตที่วัดได้ของค่า "จริง" ( y * , x * ) ซึ่งอยู่บนเส้นถดถอย:

yฉัน=yฉัน*+εฉัน,xฉัน=xฉัน*+ηฉัน,{\displaystyle {\begin{aligned}y_{i}&=y_{i}^{*}+\varepsilon _{i},\\x_{i}&=x_{i}^{*}+\eta _{i},\end{aligned}}}

โดยที่ค่าความคลาดเคลื่อนεและηเป็นอิสระต่อกัน และถือว่าทราบอัตราส่วนของความแปรปรวนของค่าความคลาดเคลื่อนทั้งสองแล้ว:

δ=σε2ση2.{\displaystyle \delta ={\frac {\sigma _{\varepsilon }^{2}}{\sigma _{\eta }^{2}}}.}

ในทางปฏิบัติ ความแปรปรวนของx{\displaystyle x}และy{\displaystyle y}ค่าพารามิเตอร์มักไม่เป็นที่ทราบ ซึ่งทำให้การประมาณค่ามีความซับซ้อนมากขึ้นδ{\displaystyle \delta }โปรดทราบว่าเมื่อใช้วิธีการวัดสำหรับx{\displaystyle x}และy{\displaystyle y}ถ้าเหมือนกัน ความแปรปรวนเหล่านี้ก็มีแนวโน้มที่จะเท่ากัน ดังนั้นδ=1{\displaystyle \delta =1}สำหรับกรณีนี้

เราพยายามค้นหาเส้นแบ่งที่ "เหมาะสมที่สุด"

y*=เบต้า0+เบต้า1x*,{\displaystyle y^{*}=\beta _{0}+\beta _{1}x^{*},}

เพื่อให้ผลรวมถ่วงน้ำหนักของค่าตกค้างกำลังสองของแบบจำลองมีค่าน้อยที่สุด: [ 3 ]

เอสเอสอาร์=ฉัน=1n(εฉัน2σε2+ηฉัน2ση2)=1σϵ2ฉัน=1n((yฉันเบต้า0เบต้า1xฉัน*)2+δ(xฉันxฉัน*)2)  นาทีเบต้า0,เบต้า1,x1*,,xn*เอสเอสอาร์{\displaystyle SSR=\sum _{i=1}^{n}{\bigg (}{\frac {\varepsilon _{i}^{2}}{\sigma _{\varepsilon }^{2}}}+{\frac {\eta _{i}^{2}}{\sigma _{\eta }^{2}}}{\bigg )}={\frac {1}{\sigma _{\epsilon }^{2}}}\sum _{i=1}^{n}{\Big (}(y_{i}-\beta _{0}-\beta _{1}x_{i}^{*})^{2}+\delta (x_{i}-x_{i}^{*})^{2}{\Big )}\ \to \ \min _{\beta _{0},\beta _{1},x_{1}^{*},\ldots ,x_{n}^{*}}SSR}

ดูJensen (2007)สำหรับที่มาของสูตรทั้งหมด

สารละลาย

คำตอบสามารถแสดงได้ในรูปของโมเมนต์ตัวอย่างระดับที่สอง กล่าวคือ เราจะคำนวณปริมาณต่อไปนี้ก่อน (ผลรวมทั้งหมดมีค่าตั้งแต่i  =  1 ถึงn ):

x¯=1nxฉันy¯=1nyฉัน,xx=1n(xฉันx¯)2=x2¯x¯2,xy=1n(xฉันx¯)(yฉันy¯)=xy¯x¯y¯,yy=1n(yฉันy¯)2=y2¯y¯2.{\displaystyle {\begin{aligned}{\overline {x}}&={\tfrac {1}{n}}\sum x_{i}&{\overline {y}}&={\tfrac {1}{n}}\sum y_{i},\\s_{xx}&={\tfrac {1}{n}}\sum (x_{i}-{\overline {x}})^{2}&&={\overline {x^{2}}}-{\overline {x}}^{2},\\s_{xy}&={\tfrac {1}{n}}\sum (x_{i}-{\overline {x}})(y_{i}-{\overline {y}})&&={\overline {xy}}-{\overline {x}}\,{\overline {y}},\\s_{yy}&={\tfrac {1}{n}}\sum (y_{i}-{\overline {y}})^{2}&&={\overline {y^{2}}}-{\overline {y}}^{2}.\end{aligned}}\,}

สุดท้าย ค่าประมาณกำลังสองน้อยที่สุดของพารามิเตอร์ของแบบจำลองจะเป็น[ 4 ]

เบต้า^1=yyδxx+(yyδxx)2+4δxy22xy,เบต้า^0=y¯เบต้า^1x¯,x^ฉัน*=xฉัน+เบต้า^1เบต้า^12+δ(yฉันเบต้า^0เบต้า^1xฉัน).{\displaystyle {\begin{aligned}&{\hat {\beta }}_{1}={\frac {s_{yy}-\delta s_{xx}+{\sqrt {(s_{yy}-\delta s_{xx})^{2}+4\delta s_{xy}^{2}}}}{2s_{xy}}},\\&{\hat {\beta }}_{0}={\overline {y}}-{\hat {\beta }}_{1}{\overline {x}},\\&{\hat {x}}_{i}^{*}=x_{i}+{\frac {{\hat {\beta }}_{1}}{{\hat {\beta }}_{1}^{2}+\delta }}(y_{i}-{\hat {\beta }}_{0}-{\hat {\beta }}_{1}x_{i}).\end{aligned}}}

การถดถอยเชิงตั้งฉาก

ในกรณีที่ความแปรปรวนของข้อผิดพลาดเท่ากัน กล่าวคือ เมื่อδ=1{\displaystyle \delta =1}การถดถอยแบบเดมิง (Deming regression) จะกลายเป็นการถดถอยเชิงตั้งฉาก (orthogonal regression ) กล่าว คือ มันจะลดผลรวมของ กำลังสองของ ระยะทางตั้งฉากจากจุดข้อมูลไปยังเส้นถดถอยให้ เหลือน้อยที่สุด ในกรณีนี้ ให้แต่ละค่าสังเกตเป็นจุดzเจ=xเจ+ฉันyเจ{\displaystyle z_{j}=x_{j}+iy_{j}}ในระนาบเชิงซ้อน (เช่น จุด)(xเจ,yเจ){\displaystyle (x_{j},y_{j})}ที่ไหนฉัน{\displaystyle i}คือหน่วยจินตนาการ ) กำหนดให้เป็นเอส=(zเจz¯)2{\displaystyle S=\sum {(z_{j}-{\overline {z}})^{2}}}ผลรวมของกำลังสองของผลต่างของจุดข้อมูลจากจุดศูนย์กลางz¯=1nzเจ{\displaystyle {\overline {z}}={\tfrac {1}{n}}\sum z_{j}}(ระบุด้วยพิกัดเชิงซ้อน) ซึ่งเป็นจุดที่มีตำแหน่งแนวนอนและแนวตั้งเป็นค่าเฉลี่ยของจุดข้อมูล จากนั้น: [ 5 ]

  • ถ้าเอส=0{\displaystyle S=0}ดังนั้น ทุกเส้นที่ลากผ่านจุดศูนย์กลางจะเป็นเส้นตั้งฉากที่เหมาะสมที่สุด
  • ถ้าเอส0{\displaystyle S\neq 0}เส้นถดถอยเชิงตั้งฉากจะผ่านจุดศูนย์กลางและขนานกับเวกเตอร์จากจุดกำเนิดไปยังจุดนั้นเอส{\displaystyle {\sqrt {S}}}.

การ แสดง ตรีโกณมิติของเส้นถดถอยเชิงตั้งฉากได้รับการนำเสนอโดยคูลิดจ์ในปี พ.ศ. 2456 [ 6 ]ระยะทางยังสามารถคำนวณได้โดยใช้สมการเส้นตรงทั่วไป ซึ่งกำหนดไว้ดังนี้y=x+เค{\displaystyle y=mx+k}.

แอปพลิเคชัน

ในกรณีที่จุดสาม จุด ไม่เรียงตัวกันบนระนาบสามเหลี่ยมที่มีจุดเหล่านี้เป็นจุดยอดจะมีวงรีสไตเนอร์ที่ ไม่ซ้ำกัน ซึ่งสัมผัสกับด้านของสามเหลี่ยมที่จุดกึ่งกลางแกนหลักของวงรีนี้จะตกอยู่บนเส้นถดถอยตั้งฉากสำหรับจุดยอดทั้งสาม[ 7 ] การหาปริมาณของ สัญญาณรบกวนภายในเซลล์ชีวภาพสามารถหาปริมาณได้โดยการใช้การถดถอยเดมิงกับพฤติกรรมที่สังเกตได้ของวงจรชีวภาพสังเคราะห์ที่มีผู้ รายงานสองคน [ 8 ]

เมื่อมนุษย์ถูกขอให้วาดการถดถอยเชิงเส้นบนแผนภาพกระจายโดยการเดา คำตอบของพวกเขาจะใกล้เคียงกับการถดถอยเชิงตั้งฉากมากกว่าการถดถอยกำลังสองน้อยที่สุดแบบธรรมดา[ 9 ]

เวอร์ชันที่ไม่ขึ้นกับมาตราส่วน

ความสัมพันธ์เชิงฟังก์ชันค่าเฉลี่ยเรขาคณิตเป็นกรณีหนึ่งของการถดถอยแบบเดมิง โดยที่δ = σ² / σ² โดยถือว่าตัวแปรทั้งสองมีความน่าเชื่อถือ เท่ากัน กล่าวคือσ / σ = σ / σ ความชันที่ได้จะเป็นค่าเฉลี่ยเรขาคณิตของความชันกำลังสองน้อยที่สุดแบบธรรมดาและความชันกำลังสองน้อยที่สุดแบบผกผัน ซึ่งทำให้คำนวณได้ง่ายโดยใช้เครื่องมือที่มีอยู่ มีการแนะนำให้ใช้ในทางชีววิทยา[ 10 ]วิธีนี้ไม่เปลี่ยนแปลงภายใต้การปรับขนาด การแปล และการสลับตัวแปร เป็นวิธีเดียวที่ขึ้นอยู่กับโมเมนต์แรกและโมเมนต์ที่สองเท่านั้นที่มีคุณสมบัตินี้[ 11 ]

การถดถอยของยอร์ก

การถดถอยของ York ขยายการถดถอยของ Deming โดยอนุญาตให้มีข้อผิดพลาดที่สัมพันธ์กันใน x และ y [ 12 ]

ดูเพิ่มเติม

ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Deming_regression&oldid=1362308451#Orthogonal_regression "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ การถดถอยของเดมิง

ในทางสถิติการถดถอยแบบเดมิง (Deming regression ) ซึ่งตั้งชื่อตามดับเบิลยู. เอ็ดเวิร์ด เดมิงเป็นแบบจำลองความคลาดเคลื่อนในตัวแปรที่พยายามหาเส้นตรงที่เหมาะสมที่สุด สำหรับ

ข้อกำหนด

สมมติว่าข้อมูลที่มีอยู่ ( y , x ) เป็นค่าสังเกตที่วัดได้ของค่า "จริง" ( y * , x * ) ซึ่งอยู่บนเส้นถดถอย:

สารละลาย

คำตอบสามารถแสดงได้ในรูปของโมเมนต์ตัวอย่างระดับที่สอง กล่าวคือ เราจะคำนวณปริมาณต่อไปนี้ก่อน (ผลรวมทั้งหมดมีค่าตั้งแต่ i = 1 ถึง n ):

การถดถอยเชิงตั้งฉาก

ในกรณีที่ความแปรปรวนของข้อผิดพลาดเท่ากัน กล่าวคือ เมื่อ δ = 1 {\displaystyle \delta =1} การถดถอยแบบเดมิง (Deming regression) จะกลายเป็นการ ถดถอยเชิงตั้งฉาก (orthogonal regression ) กล่าว คือ มันจะลดผลรวมของ กำลังสองของ...