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

อ่าน 7 นาที

การแจกแจงแบบ Dirichlet ทั่วไป

ผสานการแจกแจงก่อนหน้า/การกระจายอย่างต่อเนื่อง/การกระจายตัวของครอบครัวแบบเอ็กซ์โปเนนเชียล/การแจกแจงต่อเนื่องหลายตัวแปร

ในทางสถิติการแจกแจง Dirichlet แบบทั่วไป ( GD ) เป็นการขยายความของการแจกแจง Dirichletด้วยโครงสร้างความแปรปรวนร่วมที่ทั่วไปกว่าและมีจำนวนพารามิเตอร์เกือบสองเท่า...

การแจกแจงแบบ Dirichlet ทั่วไป

ในทางสถิติการแจกแจง Dirichlet แบบทั่วไป ( GD ) เป็นการขยายความของการแจกแจง Dirichletด้วยโครงสร้างความแปรปรวนร่วมที่ทั่วไปกว่าและมีจำนวนพารามิเตอร์เกือบสองเท่า เวกเตอร์สุ่มที่มีการแจกแจง GD เป็นกลาง โดย สมบูรณ์[ 1 ]

ฟังก์ชันความหนาแน่นของพี1,,พีเค1{\displaystyle p_{1},\ldots ,p_{k-1}}เป็น

[ฉัน=1เค1บี(เอฉัน,ฉัน)]1พีเคเค11ฉัน=1เค1[พีฉันเอฉัน1(เจ=ฉันเคพีเจ)ฉัน1(เอฉัน+ฉัน)]{\displaystyle \left[\prod _{i=1}^{k-1}B(a_{i},b_{i})\right]^{-1}p_{k}^{b_{k-1}-1}\prod _{i=1}^{k-1}\left[p_{i}^{a_{i}-1}\left(\sum _{j=i}^{k}p_{j}\right)^{b_{i-1}-(a_{i}+b_{i})}\right]}

โดยที่เรากำหนดพีเค=1ฉัน=1เค1พีฉัน{\textstyle p_{k}=1-\sum _{i=1}^{k-1}p_{i}}. ที่นี่บี(x,y){\displaystyle B(x,y)}หมายถึงฟังก์ชันเบตาซึ่งจะลดรูปเป็นฟังก์ชันการแจกแจงแบบดิริชเลต์มาตรฐานหากฉัน1=เอฉัน+ฉัน{\displaystyle b_{i-1}=a_{i}+b_{i}}สำหรับ2ฉันเค1{\displaystyle 2\leqslant i\leqslant k-1}(0{\displaystyle b_{0}}(เป็นการกำหนดโดยพลการ)

ตัวอย่างเช่น ถ้าk=4แล้วฟังก์ชันความหนาแน่นของพี1,พี2,พี3{\displaystyle p_{1},p_{2},p_{3}}เป็น

[ฉัน=13บี(เอฉัน,ฉัน)]1พี1เอ11พี2เอ21พี3เอ31พี431(พี2+พี3+พี4)1(เอ2+2)(พี3+พี4)2(เอ3+3){\displaystyle \left[\prod _{i=1}^{3}B(a_{i},b_{i})\right]^{-1}p_{1}^{a_{1}-1}p_{2}^{a_{2}-1}p_{3}^{a_{3}-1}p_{4}^{b_{3}-1}\left(p_{2}+p_{3}+p_{4}\right)^{b_{1}-\left(a_{2}+b_{2}\right)}\left(p_{3}+p_{4}\right)^{b_{2}-\left(a_{3}+b_{3}\right)}}

ที่ไหนพี1+พี2+พี3<1{\displaystyle p_{1}+p_{2}+p_{3}<1}และพี4=1พี1พี2พี3{\displaystyle p_{4}=1-p_{1}-p_{2}-p_{3}}.

คอนเนอร์และโมซิมานน์กำหนดนิยามของ PDF ในลักษณะดังกล่าวด้วยเหตุผลดังต่อไปนี้ กำหนดนิยามของตัวแปรสุ่มz1,,zเค1{\displaystyle z_{1},\ldots ,z_{k-1}}กับz1=พี1,z2=พี2/(1พี1),z3=พี3/(1(พี1+พี2)),,zฉัน=พีฉัน/(1(พี1++พีฉัน1)){\displaystyle z_{1}=p_{1},z_{2}=p_{2}/\left(1-p_{1}\right),z_{3}=p_{3}/\left(1-(p_{1}+p_{2})\right),\ldots ,z_{i}=p_{i}/\left(1-\left(p_{1}+\cdots +p_{i-1}\right)\right)}. แล้วพี1,,พีเค{\displaystyle p_{1},\ldots ,p_{k}}มีการกระจายแบบ Dirichlet ทั่วไปตามพารามิเตอร์ที่กำหนดข้างต้น หากzฉัน{\displaystyle z_{i}}เป็นเบต้า อิสระ ที่มีพารามิเตอร์เอฉัน,ฉัน{\displaystyle a_{i},b_{i}},ฉัน=1,,เค1{\displaystyle i=1,\ldots ,k-1}.

รูปแบบทางเลือกที่เสนอโดยหว่อง

Wong [ 2 ]ให้รูปแบบที่กระชับกว่าเล็กน้อยสำหรับx1++xเค1{\displaystyle x_{1}+\cdots +x_{k}\leq 1}

ฉัน=1เคxฉันαฉัน1(1x1xฉัน)γฉันบี(αฉัน,เบต้าฉัน){\displaystyle \prod _{i=1}^{k}{\frac {x_{i}^{\alpha _{i}-1}\left(1-x_{1}-\cdots -x_{i}\right)^{\gamma _{i}}}{B(\alpha _{i},\beta _{i})}}}

ที่ไหนγเจ=เบต้าเจαเจ+1เบต้าเจ+1{\displaystyle \gamma _{j}=\beta _{j}-\alpha _{j+1}-\beta _{j+1}}สำหรับ1เจเค1{\displaystyle 1\leq j\leq k-1}และγเค=เบต้าเค1{\displaystyle \gamma _{k}=\beta _{k}-1}โปรดทราบว่า หว่องกำหนดการกระจายตัวเหนือ aเค{\displaystyle k}พื้นที่มิติ (โดยนัยกำหนด)xเค+1=1ฉัน=1เคxฉัน{\textstyle x_{k+1}=1-\sum _{i=1}^{k}x_{i}}) ในขณะที่คอนเนอร์และโมซิแมนใช้เค1{\displaystyle k-1}พื้นที่มิติที่มีxเค=1ฉัน=1เค1xฉัน{\textstyle x_{k}=1-\sum _{i=1}^{k-1}x_{i}}.

ฟังก์ชันโมเมนต์ทั่วไป

ถ้าX=(X1,,Xเค)~จีดีเค(α1,,αเค;เบต้า1,,เบต้าเค){\displaystyle X=\left(X_{1},\ldots ,X_{k}\right)\sim GD_{k}\left(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k}\right)}, แล้ว

อี[X11X22Xเคเค]=เจ=1เคΓ(αเจ+เบต้าเจ)Γ(αเจ+เจ)Γ(เบต้าเจ+δเจ)Γ(αเจ)Γ(เบต้าเจ)Γ(αเจ+เบต้าเจ+เจ+δเจ){\displaystyle E\left[X_{1}^{r_{1}}X_{2}^{r_{2}}\cdots X_{k}^{r_{k}}\right]=\prod _{j=1}^{k}{\frac {\Gamma \left(\alpha _{j}+\beta _{j}\right)\Gamma \left(\alpha _{j}+r_{j}\right)\Gamma \left(\beta _{j}+\delta _{j}\right)}{\Gamma \left(\alpha _{j}\right)\Gamma \left(\beta _{j}\right)\Gamma \left(\alpha _{j}+\beta _{j}+r_{j}+\delta _{j}\right)}}}

ที่ไหนδเจ=เจ+1+เจ+2++เค{\displaystyle \delta _{j}=r_{j+1}+r_{j+2}+\cdots +r_{k}}สำหรับเจ=1,2,,เค1{\displaystyle j=1,2,\cdots ,k-1}และδเค=0{\displaystyle \delta _{k}=0}. ดังนั้น

อี(Xเจ)=αเจαเจ+เบต้าเจ=1เจ1เบต้าα+เบต้า.{\displaystyle E\left(X_{j}\right)={\frac {\alpha _{j}}{\alpha _{j}+\beta _{j}}}\prod _{m=1}^{j-1}{\frac {\beta _{m}}{\alpha _{m}+\beta _{m}}}.}

ลดรูปเป็นการกระจายแบบ Dirichlet มาตรฐาน

ดังที่กล่าวไว้ข้างต้น หากฉัน1=เอฉัน+ฉัน{\displaystyle b_{i-1}=a_{i}+b_{i}}สำหรับ2ฉันเค{\displaystyle 2\leq i\leq k}จากนั้นการแจกแจงจะลดลงเหลือการแจกแจงแบบ Dirichlet มาตรฐาน เงื่อนไขนี้แตกต่างจากกรณีปกติ ซึ่งการตั้งค่าพารามิเตอร์เพิ่มเติมของการแจกแจงแบบทั่วไปให้เป็นศูนย์จะทำให้ได้การแจกแจงแบบเดิม อย่างไรก็ตาม ในกรณีของ GDD นั้น จะส่งผลให้ได้ฟังก์ชันความหนาแน่นที่ซับซ้อนมาก

การวิเคราะห์แบบเบย์เซียน

สมมติX=(X1,,Xเค)~จีดีเค(α1,,αเค;เบต้า1,,เบต้าเค){\displaystyle X=\left(X_{1},\ldots ,X_{k}\right)\sim GD_{k}\left(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k}\right)}เป็น Dirichlet แบบทั่วไป และว่าวายX{\displaystyle Y\mid X}เป็นพหุนามที่มีn{\displaystyle n}การทดลอง (ที่นี่)วาย=(วาย1,,วายเค){\displaystyle Y=\left(Y_{1},\ldots ,Y_{k}\right)}). การเขียนวายเจ=yเจ{\displaystyle Y_{j}=y_{j}}สำหรับ1เจเค{\displaystyle 1\leq j\leq k}และyเค+1=nฉัน=1เคyฉัน{\textstyle y_{k+1}=n-\sum _{i=1}^{k}y_{i}}ข้อต่อด้านหลังของX|วาย{\displaystyle X|Y}เป็นการแจกแจงแบบ Dirichlet ทั่วไปที่มี

Xวาย~จีดีเค(α1,,αเค;เบต้า1,,เบต้าเค){\displaystyle X\mid Y\sim GD_{k}\left({\alpha '}_{1},\ldots ,{\alpha '}_{k};{\beta '}_{1},\ldots ,{\beta '}_{k}\right)}

ที่ไหนαเจ=αเจ+yเจ{\displaystyle {\alpha '}_{j}=\alpha _{j}+y_{j}}และเบต้าเจ=เบต้าเจ+ฉัน=เจ+1เค+1yฉัน{\displaystyle {\beta '}_{j}=\beta _{j}+\sum _{i=j+1}^{k+1}y_{i}}สำหรับ1เจเค.{\displaystyle 1\leqslant j\leqslant k.}

การทดลองสุ่มตัวอย่าง

หว่องได้ยกตัวอย่างระบบต่อไปนี้เพื่อแสดงความแตกต่างระหว่างการแจกแจงแบบ Dirichlet และการแจกแจงแบบ Dirichlet ทั่วไป เขาตั้งสมมติฐานว่าโถขนาดใหญ่บรรจุลูกบอลจำนวนหนึ่งเค+1{\displaystyle k+1}สีต่าง ๆ กัน สัดส่วนของแต่ละสีไม่ทราบแน่ชัด จงเขียนX=(X1,,Xเค){\displaystyle X=(X_{1},\ldots ,X_{k})}สำหรับสัดส่วนของลูกบอลที่มีสีเจ{\displaystyle j}ในโกศ

การทดลองที่ 1นักวิเคราะห์คนที่ 1 เชื่อว่าX~ดี(α1,,αเค,αเค+1){\displaystyle X\sim D(\alpha _{1},\ldots ,\alpha _{k},\alpha _{k+1})}(เช่น,X{\displaystyle X}เป็นสมการ Dirichlet ที่มีพารามิเตอร์αฉัน{\displaystyle \alpha _{i}}จากนั้นนักวิเคราะห์ก็ทำการเค+1{\displaystyle k+1}กล่องแก้วและวางαฉัน{\displaystyle \alpha _{i}}ลูกแก้วหลากสีฉัน{\displaystyle i}ในกล่องฉัน{\displaystyle i}(สันนิษฐานว่าαฉัน{\displaystyle \alpha _{i}}เป็นจำนวนเต็ม1{\displaystyle \geq 1}จากนั้นนักวิเคราะห์คนที่ 1 หยิบลูกบอลจากโถ แล้วสังเกตสีของลูกบอล (เช่น สี)เจ{\displaystyle j}) แล้วใส่ลงในกล่องเจ{\displaystyle j}เขาสามารถระบุกล่องที่ถูกต้องได้เพราะกล่องเหล่านั้นโปร่งใสและมองเห็นสีของลูกแก้วที่อยู่ข้างในได้ กระบวนการนี้ดำเนินต่อไปจนกระทั่งn{\displaystyle n}ลูกบอลถูกสุ่มเลือกแล้ว การแจกแจงความน่าจะเป็นภายหลังจะเป็นแบบ Dirichlet โดยมีพารามิเตอร์เป็นจำนวนลูกแก้วในแต่ละกล่อง

การทดลองที่ 2นักวิเคราะห์คนที่ 2 เชื่อว่าX{\displaystyle X}เป็นไปตามการแจกแจงแบบ Dirichlet ทั่วไป:X~จีดี(α1,,αเค;เบต้า1,,เบต้าเค){\displaystyle X\sim GD(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k})}พารามิเตอร์ทั้งหมดถือว่าเป็นจำนวนเต็มบวกอีกครั้ง นักวิเคราะห์ทำการวิเคราะห์เค+1{\displaystyle k+1}กล่องไม้ กล่องเหล่านี้มีสองส่วน: ส่วนหนึ่งสำหรับลูกบอล และอีกส่วนสำหรับลูกแก้ว ลูกบอลมีสี แต่ลูกแก้วไม่มีสี จากนั้นสำหรับเจ=1,,เค{\displaystyle j=1,\ldots ,k}เขาวางαเจ{\displaystyle \alpha _{j}}ลูกบอลหลากสีเจ{\displaystyle j}, และเบต้าเจ{\displaystyle \beta _{j}}ลูกแก้ว ใส่ลงในกล่องเจ{\displaystyle j}จากนั้นเขาก็ใส่ลูกบอลสีลงไปเค+1{\displaystyle k+1}ในกล่องเค+1{\displaystyle k+1}จากนั้น นักวิเคราะห์จะหยิบลูกบอลจากโถ เนื่องจากกล่องทำจากไม้ นักวิเคราะห์จึงไม่สามารถบอกได้ว่าควรใส่ลูกบอลลงในกล่องใด (เหมือนกับในทดลองที่ 1 ข้างต้น) นอกจากนี้ เขายังมีความจำไม่ดีและจำไม่ได้ว่ากล่องใดมีลูกบอลสีอะไร เขาต้องค้นหาว่ากล่องใดเป็นกล่องที่ถูกต้องที่จะใส่ลูกบอลลงไป เขาทำเช่นนั้นโดยการเปิดกล่องที่ 1 และเปรียบเทียบลูกบอลในกล่องนั้นกับลูกบอลที่หยิบได้ หากสีแตกต่างกัน แสดงว่ากล่องนั้นผิด นักวิเคราะห์จะใส่ลูกบอลลงในกล่องที่ 1 แล้วไปยังกล่องที่ 2 เขาทำซ้ำกระบวนการนี้จนกว่าลูกบอลในกล่องจะตรงกับลูกบอลที่หยิบได้ ณ จุดนั้น เขาจะใส่ลูกบอลลงในกล่องที่มีลูกบอลสีเดียวกัน จากนั้น นักวิเคราะห์จะหยิบลูกบอลอีกลูกจากโถและทำซ้ำจนกระทั่งn{\displaystyle n}มีการวาดลูกบอล จากนั้นจึงใช้ Dirichlet แบบทั่วไปที่มีพารามิเตอร์α{\displaystyle \alpha }โดยเป็นจำนวนลูกบอล และเบต้า{\displaystyle \beta }จำนวนลูกแก้วในแต่ละกล่อง

โปรดสังเกตว่าในการทดลองที่ 2 การเปลี่ยนลำดับของกล่องมีผลกระทบที่ไม่เล็กน้อย ซึ่งแตกต่างจากการทดลองที่ 1

ดูเพิ่มเติม

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

สรุปเนื้อหา

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

ข้อมูลสำคัญเกี่ยวกับ การแจกแจงแบบ Dirichlet ทั่วไป

ในทางสถิติการแจกแจง Dirichlet แบบทั่วไป ( GD ) เป็นการขยายความของการแจกแจง Dirichletด้วยโครงสร้างความแปรปรวนร่วมที่ทั่วไปกว่าและมีจำนวนพารามิเตอร์เกือบสองเท่า...

รูปแบบทางเลือกที่เสนอโดยหว่อง

Wong [ 2 ] ให้รูปแบบที่กระชับกว่าเล็กน้อยสำหรับ x 1 + ⋯ + x เค ≤ 1 {\displaystyle x_{1}+\cdots +x_{k}\leq 1}

ฟังก์ชันโมเมนต์ทั่วไป

ถ้า X = ( X 1 , … , X เค ) ~ จี ดี เค ( α 1 , … , α เค ; เบต้า 1 , … , เบต้า เค ) {\displaystyle X=\left(X_{1},\ldots ,X_{k}\right)\sim GD_{k}\left(\alpha _{1},\ldots ,\alpha _{k};\beta _{1},\ldots ,\beta _{k}\right)} , แล้ว

ลดรูปเป็นการกระจายแบบ Dirichlet มาตรฐาน

ดังที่กล่าวไว้ข้างต้น หาก ข ฉัน − 1 = เอ ฉัน + ข ฉัน {\displaystyle b_{i-1}=a_{i}+b_{i}} สำหรับ 2 ≤ ฉัน ≤ เค {\displaystyle 2\leq i\leq k} จากนั้นการแจกแจงจะลดลงเหลือการแจกแจงแบบ Dirichlet มาตรฐาน เงื่อนไขนี้แตกต่างจากกรณีปกติ...