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

อ่าน 11 นาที

ไม่มีชื่อบทความ

ใน ทางสถิติ โดยเฉพาะอย่างยิ่งใน การทดสอบสมมติฐาน การแจกแจง Hotelling's T-squared ( T² ) ที่ เสนอ โดย Harold Hotelling [ 1 ] เป็นการ แจกแจง ความน่าจะเป็นแบบหลายตัวแปร...

การแจกแจง T -squared ของ Hotelling

ในทางสถิติโดยเฉพาะอย่างยิ่งในการทดสอบสมมติฐานการแจกแจง Hotelling's T-squared ( ) ที่เสนอโดย Harold Hotelling [ 1 ] เป็นการแจกแจงความน่าจะเป็นแบบหลายตัวแปรที่เกี่ยวข้องอย่างใกล้ชิดกับการแจกแจงFและโดดเด่นที่สุดเนื่องจากเกิดขึ้นจากการแจกแจงของชุดสถิติตัวอย่างที่เป็นการวางนัยทั่วไปตามธรรมชาติของสถิติที่อยู่เบื้องหลังการแจกแจงStudent's tสถิติ Hotelling's t-squared ( )เป็นการวางนัยทั่วไปของสถิติStudent's tที่ใช้ใน การ ทดสอบสมมติฐานแบบหลายตัวแปร[ 2 ]

แรงจูงใจ

การแจกแจงนี้เกิดขึ้นในสถิติหลายตัวแปรในการทดสอบความแตกต่างระหว่างค่าเฉลี่ย (หลายตัวแปร) ของประชากรที่แตกต่างกัน โดยที่การทดสอบสำหรับปัญหาตัวแปรเดียวจะใช้การทดสอบtการแจกแจงนี้ตั้งชื่อตามHarold Hotelling ซึ่งพัฒนาให้เป็นการวางนัยทั่วไปของ การแจกแจงtของ Student [ 1 ]

คำนิยาม

ถ้าเวกเตอร์ {\displaystyle d} เป็นการแจกแจงแบบเกาส์เซียนหลายตัวแปรโดยมีค่าเฉลี่ยเป็นศูนย์และเมทริกซ์ความแปรปรวน ร่วมเป็นหนึ่งเอ็น(0พี,ฉันพี,พี){\displaystyle N(\mathbf {0} _{p},\mathbf {I} _{p,p})}และเอ็ม{\displaystyle M}เป็นพี×พี{\displaystyle p\times p}เมทริกซ์สุ่มที่มีการแจกแจงแบบวิชาร์ต(ฉันพี,พี,){\displaystyle W(\mathbf {I} _{p,p},m)}โดยมี เมทริกซ์มาตราส่วนหน่วยและm องศาอิสระและdกับMเป็นอิสระต่อกันรูปแบบกำลังสอง จะเป็นดังนี้X{\displaystyle X}มีการกระจายแบบ Hotelling (พร้อมพารามิเตอร์)พี{\displaystyle p}และ{\displaystyle m}): [ 3 ]X=ทีเอ็ม1~ที2(พี,).{\displaystyle X=md^{T}M^{-1}d\sim T^{2}(p,m).}

สามารถแสดงได้ว่า ถ้าตัวแปรสุ่มXมีการแจกแจงแบบ Hotelling's T -squaredX~ทีพี,2{\displaystyle X\sim T_{p,m}^{2}}จากนั้น: [ 1 ]พี+1พีX~เอฟพี,พี+1{\displaystyle {\frac {m-p+1}{pm}}X\sim F_{p,m-p+1}} ที่ไหนเอฟพี,พี+1{\displaystyle F_{p,m-p+1}}คือการแจกแจงFที่มีพารามิเตอร์pและm p + 1    

สถิติt -squared ของ Hotelling

อนุญาตΣ^{\displaystyle {\hat {\mathbf {\Sigma } }}}ให้ค่าความแปรปรวนร่วมของตัวอย่าง เป็น :

Σ^=1n1ฉัน=1n(xฉันx¯)(xฉันx¯){\displaystyle {\hat {\mathbf {\Sigma } }}={\frac {1}{n-1}}\sum _{i=1}^{n}\left(\mathbf {x} _{i}-{\overline {\mathbf {x} }}\right)\left(\mathbf {x} _{i}-{\overline {\mathbf {x} }}\right)'}

โดยที่เราใช้เครื่องหมายอะพอสโทรฟีแทนการสลับแถวและ คอลัมน์ สามารถแสดงได้ว่าΣ^{\displaystyle {\hat {\mathbf {\Sigma } }}}เป็น เมทริกซ์ บวก (กึ่ง) กำหนดและ(n1)Σ^{\displaystyle (n-1){\hat {\mathbf {\Sigma } }}}เป็นไปตามการแจกแจง Wishart แบบ pตัวแปรที่มีองศาอิสระn − 1 [ 4 ] เมทริกซ์ความแปรปรวนร่วมของค่าเฉลี่ยของตัวอย่างมีดังนี้  Σ^x¯=Σ^/n{\displaystyle {\hat {\mathbf {\Sigma } }}_{\overline {\mathbf {x} }}={\hat {\mathbf {\Sigma } }}/n}[ 5 ]

สถิติt -squared ของ Hotellingถูกกำหนดดังนี้: [ 6 ]

ที2=(x¯μ)Σ^x¯1(x¯μ)=n(x¯μ)Σ^1(x¯μ),{\displaystyle t^{2}=({\overline {\mathbf {x} }}-{\boldsymbol {\mu }})'{\hat {\mathbf {\Sigma } }}_{\overline {\mathbf {x} }}^{-1}({\overline {\mathbf {x} }}-{\boldsymbol {\mathbf {\mu } }})=n({\overline {\mathbf {x} }}-{\boldsymbol {\mu }})'{\hat {\mathbf {\Sigma } }}^{-1}({\overline {\mathbf {x} }}-{\boldsymbol {\mathbf {\mu } }}),}

ซึ่งเป็นสัดส่วนกับระยะทางมาฮาลาโนบิสระหว่างค่าเฉลี่ยของตัวอย่างและμ{\displaystyle {\boldสัญลักษณ์ {\mu }}}ด้วยเหตุนี้ จึงควรคาดหวังว่าค่าสถิติจะมีค่าต่ำหากx¯μ{\displaystyle {\overline {\mathbf {x} }}\approx {\boldsymbol {\mu }}}และจะมีค่าสูงหากค่าเหล่านั้นแตกต่างกัน

จากการกระจายตัว

ที2~ทีพี,n12=พี(n1)nพีเอฟพี,nพี,{\displaystyle t^{2}\sim T_{p,n-1}^{2}={\frac {p(n-1)}{np}}F_{p,np},}

ที่ไหนเอฟพี,nพี{\displaystyle F_{p,np}}คือการแจกแจงFที่มีพารามิเตอร์pและn p 

เพื่อคำนวณค่าp (ซึ่งไม่เกี่ยวข้องกับ ตัวแปร pในที่นี้) โปรดสังเกตว่าการกระจายของที2{\displaystyle t^{2}}เทียบเท่ากับหมายความว่า

nพีพี(n1)ที2~เอฟพี,nพี.{\displaystyle {\frac {np}{p(n-1)}}t^{2}\sim F_{p,np}.}

จากนั้น ใช้ค่าทางด้านซ้ายมือเพื่อประเมิน ค่า pที่สอดคล้องกับตัวอย่าง ซึ่งได้มาจากค่าการแจกแจงFนอกจากนี้ยังสามารถกำหนดช่วงความเชื่อมั่น ได้โดยใช้ตรรกะที่คล้ายกัน

แรงจูงใจ

อนุญาตเอ็นพี(μ,Σ){\displaystyle {\mathcal {N}__{p}({\boldสัญลักษณ์ {\mu }},{\mathbf {\Sigma } })}แสดงถึงการแจกแจงปกติp ตัวแปร ที่มีตำแหน่งμ{\displaystyle {\boldสัญลักษณ์ {\mu }}}และค่าความแปรปรวน ร่วมที่ทราบΣ{\displaystyle {\mathbf {\Sigma } }}. อนุญาต

x1,,xn~เอ็นพี(μ,Σ){\displaystyle {\mathbf {x} __{1},\dots ,{\mathbf {x} __{n}\sim {\mathcal {N}__{p}({\boldสัญลักษณ์ {\mu }},{\mathbf {\Sigma } })}

เป็นตัวแปรสุ่มอิสระที่มีการแจกแจงเหมือนกัน (iid) จำนวนn ตัว ซึ่งสามารถแสดงได้ดังนี้พี×1{\displaystyle p\times 1}เวกเตอร์คอลัมน์ของจำนวนจริง กำหนด

x¯=x1++xnn{\displaystyle {\overline {\mathbf {x} }}={\frac {\mathbf {x} _{1}+\cdots +\mathbf {x} _{n}}{n}}}

เพื่อเป็นค่าเฉลี่ยตัวอย่างที่มีความแปรปรวนร่วมΣx¯=Σ/n{\displaystyle {\mathbf {\Sigma } }_{\overline {\mathbf {x} }}={\mathbf {\Sigma } }/n}สามารถแสดงให้เห็นได้ว่า

(x¯μ)Σx¯1(x¯μ)~χพี2,{\displaystyle ({\overline {\mathbf {x} }}-{\boldsymbol {\mu }})'{\mathbf {\Sigma } }_{\overline {\mathbf {x} }}^{-1}({\overline {\mathbf {x} }}-{\boldsymbol {\mathbf {\mu } }})\sim \chi _{p}^{2},}

ที่ไหนχพี2{\displaystyle \chi _{p}^{2}}คือการแจกแจงไคกำลังสองที่มีองศาอิสระp [ 7 ]

อีกทางเลือกหนึ่งคือ สามารถโต้แย้งโดยใช้ฟังก์ชันความหนาแน่นและฟังก์ชันลักษณะเฉพาะได้ ดังต่อไปนี้

สถิติสองกลุ่มตัวอย่าง

ถ้าx1,,xnx~เอ็นพี(μ,Σ){\displaystyle {\mathbf {x} }_{1},\dots ,{\mathbf {x} }_{n_{x}}\sim N_{p}({\boldsymbol {\mu }},{\mathbf {\Sigma } })}และy1,,yny~เอ็นพี(μ,Σ){\displaystyle {\mathbf {y} }_{1},\dots ,{\mathbf {y} }_{n_{y}}\sim N_{p}({\boldsymbol {\mu }},{\mathbf {\Sigma } })}โดยที่ตัวอย่าง ถูกสุ่ม อย่างอิสระจากสองการแจกแจงปกติหลายตัวแปรอิสระที่ มีค่าเฉลี่ยและความแปรปรวนร่วมเดียวกัน และเรากำหนด

x¯=1nxฉัน=1nxxฉันy¯=1nyฉัน=1nyyฉัน{\displaystyle {\overline {\mathbf {x} }}={\frac {1}{n_{x}}}\sum _{i=1}^{n_{x}}\mathbf {x} _{i}\qquad {\overline {\mathbf {y} }}={\frac {1}{n_{y}}}\sum _{i=1}^{n_{y}}\mathbf {y} _{i}}

ตามค่าเฉลี่ยของตัวอย่าง และ

Σ^x=1nx1ฉัน=1nx(xฉันx¯)(xฉันx¯)Σ^y=1ny1ฉัน=1ny(yฉันy¯)(yฉันy¯){\displaystyle {\begin{aligned}{\hat {\mathbf {\Sigma } }}_{\mathbf {x} }&={\frac {1}{n_{x}-1}}\sum _{i=1}^{n_{x}}\left(\mathbf {x} _{i}-{\overline {\mathbf {x} }}\right)\left(\mathbf {x} _{i}-{\overline {\mathbf {x} }}\right)'\\{\hat {\mathbf {\Sigma } }}_{\mathbf {y} }&={\frac {1}{n_{y}-1}}\sum _{i=1}^{n_{y}}\left(\mathbf {y} _{i}-{\overline {\mathbf {y} }}\right)\left(\mathbf {y} _{i}-{\overline {\mathbf {y} }}\right)'\end{aligned}}}

โดยใช้เมทริกซ์ความแปรปรวนร่วมของตัวอย่างตามลำดับ จากนั้น

Σ^=(nx1)Σ^x+(ny1)Σ^ynx+ny2{\displaystyle {\hat {\mathbf {\Sigma } }}={\frac {(n_{x}-1){\hat {\mathbf {\Sigma } }}_{\mathbf {x} }+(n_{y}-1){\hat {\mathbf {\Sigma } }}_{\mathbf {y} }}{n_{x}+n_{y}-2}}}

คือ ค่าประมาณ เมทริกซ์ความแปรปรวนร่วมแบบรวม ที่ไม่เอนเอียง (ซึ่งเป็นส่วนขยายของความแปรปรวนแบบรวม )

สุดท้ายนี้ ค่าสถิติ t -squared สองกลุ่มตัวอย่างของ Hotellingคือ

ที2=nxnynx+ny(x¯y¯)Σ^1(x¯y¯)~ที2(พี,nx+ny2){\displaystyle t^{2}={\frac {n_{x}n_{y}}{n_{x}+n_{y}}}({\overline {\mathbf {x} }}-{\overline {\mathbf {y} }})'{\hat {\mathbf {\Sigma } }}^{-1}({\overline {\mathbf {x} }}-{\overline {\mathbf {y} }})\sim T^{2}(p,n_{x}+n_{y}-2)}

สามารถเชื่อมโยงกับการแจกแจง F ได้โดย[ 4 ​​]

nx+nyพี1(nx+ny2)พีที2~เอฟ(พี,nx+ny1พี).{\displaystyle {\frac {n_{x}+n_{y}-p-1}{(n_{x}+n_{y}-2)p}}t^{2}\sim F(p,n_{x}+n_{y}-1-p).}

การแจกแจงที่ไม่เป็นศูนย์ของสถิตินี้คือการแจกแจง F แบบไม่ศูนย์กลาง (อัตราส่วนของตัวแปรสุ่มไคกำลังสองแบบไม่ศูนย์กลาง และตัวแปรสุ่ม ไคกำลังสองแบบ ศูนย์กลางที่เป็นอิสระ ) nx+nyพี1(nx+ny2)พีที2~เอฟ(พี,nx+ny1พี;δ),{\displaystyle {\frac {n_{x}+n_{y}-p-1}{(n_{x}+n_{y}-2)p}}t^{2}\sim F(p,n_{x}+n_{y}-1-p;\delta ),} กับ δ=nxnynx+nyΣ1,{\displaystyle \delta ={\frac {n_{x}n_{y}}{n_{x}+n_{y}}}{\boldsymbol {d}}'\mathbf {\Sigma } ^{-1}{\boldsymbol {d}},} ที่ไหน=x¯y¯{\displaystyle {\boldsymbol {d}}=\mathbf {{\overline {x}}-{\overline {y}}} }คือเวกเตอร์ผลต่างระหว่างค่าเฉลี่ยของประชากร

ในกรณีที่มีตัวแปรสองตัว สูตรจะง่ายขึ้นอย่างเห็นได้ชัด ทำให้เข้าใจถึงความสัมพันธ์ได้ง่ายขึ้นρ{\displaystyle \rho }ระหว่างตัวแปรต่างๆ นั้นมีผลกระทบที2{\displaystyle t^{2}}ถ้าเรากำหนด 1=x¯1y¯1,2=x¯2y¯2{\displaystyle d_{1}={\overline {x}}_{1}-{\overline {y}}_{1},\qquad d_{2}={\overline {x}}_{2}-{\overline {y}}_{2}} และ 1=Σ112=Σ22ρ=Σ12/(12)=Σ21/(12){\displaystyle s_{1}={\sqrt {\Sigma _{11}}}\qquad s_{2}={\sqrt {\Sigma _{22}}}\qquad \rho =\Sigma _{12}/(s_{1}s_{2})=\Sigma _{21}/(s_{1}s_{2})} แล้ว ที2=nxny(nx+ny)(1ρ2)[(11)2+(22)22ρ(11)(22)]{\displaystyle t^{2}={\frac {n_{x}n_{y}}{(n_{x}+n_{y})(1-\rho ^{2})}}\left[\left({\frac {d_{1}}{s_{1}}}\right)^{2}+\left({\frac {d_{2}}{s_{2}}}\right)^{2}-2\rho \left({\frac {d_{1}}{s_{1}}}\right)\left({\frac {d_{2}}{s_{2}}}\right)\right]} ดังนั้น หากความแตกต่างในสองแถวของเวกเตอร์=x¯y¯{\displaystyle \mathbf {d} ={\overline {\mathbf {x} }}-{\overline {\mathbf {y} }}}โดยทั่วไปแล้วจะมีราศีเดียวกันที2{\displaystyle t^{2}}เล็กลงเรื่อยๆρ{\displaystyle \rho }จะกลายเป็นค่าบวกมากขึ้น หากความแตกต่างมีเครื่องหมายตรงข้ามกันที2{\displaystyle t^{2}}มีขนาดใหญ่ขึ้นเมื่อρ{\displaystyle \rho }กลายเป็นเชิงบวกมากขึ้น

กรณีพิเศษของตัวแปรเดี่ยวสามารถพบได้ในการทดสอบ t ของเวลช์ (Welch's t-test )

มีการเสนอการทดสอบที่แข็งแกร่งและมีประสิทธิภาพมากกว่าการทดสอบสองตัวอย่างของ Hotelling ในเอกสารทางวิชาการ ดูตัวอย่างเช่น การทดสอบตามระยะห่างระหว่างจุด ซึ่งสามารถนำไปใช้ได้แม้ในกรณีที่จำนวนตัวแปรเทียบเท่าหรือมากกว่าจำนวนผู้ถูกทดสอบ[ 9 ] [ 10 ]

ดูเพิ่มเติม

สรุปเนื้อหา

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

ข้อมูลสำคัญเกี่ยวกับ ไม่มีชื่อบทความ

ใน ทางสถิติ โดยเฉพาะอย่างยิ่งใน การทดสอบสมมติฐาน การแจกแจง Hotelling's T-squared ( T² ) ที่ เสนอ โดย Harold Hotelling [ 1 ] เป็นการ แจกแจง ความน่าจะเป็นแบบหลายตัวแปร...

แรงจูงใจ

การแจกแจงนี้เกิดขึ้นใน สถิติหลายตัวแปร ในการ ทดสอบ ความแตกต่างระหว่างค่าเฉลี่ย (หลายตัวแปร) ของประชากรที่แตกต่างกัน โดยที่การทดสอบสำหรับปัญหาตัวแปรเดียวจะใช้การ ทดสอบ t การแจกแจงนี้ตั้งชื่อตาม Harold Hotelling ซึ่งพัฒนาให้เป็นการวางนัยทั่วไปของ การแจกแจง t...

คำนิยาม

ถ้าเวกเตอร์ ง {\displaystyle d} เป็นการ แจกแจงแบบเกาส์เซียนหลายตัวแปร โดยมีค่าเฉลี่ยเป็นศูนย์และ เมทริกซ์ความแปรปรวน ร่วมเป็นหนึ่ง เอ็น ( 0 พี , ฉัน พี , พี ) {\displaystyle N(\mathbf {0} _{p},\mathbf {I} _{p,p})} และ เอ็ม {\displaystyle M} เป็น พี × พี...

สถิติ t -squared ของ Hotelling

อนุญาต Σ ^ {\displaystyle {\hat {\mathbf {\Sigma } }}} ให้ค่า ความแปรปรวนร่วมของตัวอย่าง เป็น :