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

อ่าน 14 นาที

กระบวนการร้านอาหารจีน

ใน ทฤษฎีความน่าจะเป็น กระบวนการ ร้าน อาหารจีน เป็น กระบวนการสุ่มแบบ เวลาไม่ต่อเนื่อง คล้ายกับการจัดที่นั่งให้ลูกค้าที่โต๊ะในร้านอาหาร ลองนึกภาพร้านอาหารที่มีโต๊ะกลมจำนวนอนันต์...

กระบวนการร้านอาหารจีน

ในทฤษฎีความน่าจะเป็น กระบวนการ ร้านอาหารจีนเป็นกระบวนการสุ่มแบบเวลาไม่ต่อเนื่อง คล้ายกับการจัดที่นั่งให้ลูกค้าที่โต๊ะในร้านอาหาร ลองนึกภาพร้านอาหารที่มีโต๊ะกลมจำนวนอนันต์ แต่ละโต๊ะมีความจุไม่จำกัด ลูกค้าคนที่ 1 นั่งที่โต๊ะแรก ลูกค้าคนถัดไปอาจนั่งที่โต๊ะเดียวกับลูกค้าคนที่ 1 หรือโต๊ะถัดไป กระบวนการนี้ดำเนินต่อไป โดยลูกค้าแต่ละคนเลือกที่จะนั่งที่โต๊ะที่มีคนนั่งอยู่แล้ว ด้วยความน่าจะเป็นที่แปรผันตามจำนวนลูกค้าที่นั่งอยู่แล้ว (กล่าวคือ พวกเขามีโอกาสนั่งที่โต๊ะที่มีคนเยอะมากกว่าโต๊ะที่มีคนน้อย) หรือโต๊ะว่าง ณ เวลาnลูกค้าnคนจะถูกแบ่งออกเป็นm nโต๊ะ (หรือกลุ่มของส่วนแบ่ง) ผลลัพธ์ของกระบวนการนี้สามารถสลับเปลี่ยนได้หมายความว่าลำดับที่ลูกค้านั่งไม่มีผลต่อความน่าจะเป็นของการกระจาย สุดท้าย คุณสมบัตินี้ช่วยลดความซับซ้อนของปัญหาหลายอย่างในพันธุศาสตร์ประชากรการวิเคราะห์ทางภาษาศาสตร์และการจดจำภาพได้ อย่างมาก  

การเปรียบเทียบกับร้านอาหารปรากฏครั้งแรกในงานเขียนของDavid Aldous ในปี 1985 [ 1 ] ซึ่งอ้างอิงถึงJim Pitman (ซึ่งให้เครดิตLester Dubins เพิ่มเติมด้วย ) [ 2 ]

กระบวนการแบ่งส่วนที่เทียบเท่ากันได้รับการตีพิมพ์เมื่อปีก่อนโดยFred Hoppe [ 3 ]โดยใช้ "แผนผังโถ" ที่คล้ายกับโถของ Pólya เมื่อเปรียบเทียบกับ แบบจำลองโถของ Hoppe กระบวนการร้านอาหารจีนมีข้อได้เปรียบตรงที่สามารถอธิบายการเรียงสับเปลี่ยนแบบสุ่มผ่านโครงสร้างวงจรได้อย่างเป็นธรรมชาติ นอกเหนือจากการอธิบายการแบ่งส่วนแบบสุ่ม

คำจำกัดความอย่างเป็นทางการ

สำหรับจำนวนเต็มบวกใดๆn{\displaystyle n}, อนุญาตพีn{\displaystyle {\mathcal {P}}_{n}}แทนเซตของการแบ่งย่อยทั้งหมดของเซต{1,2,3,...,n}[n]{\displaystyle \{1,2,3,...,n\}\triangleq [n]}กระบวนการของร้านอาหารจีนนั้นใช้ค่าในผลคูณคาร์ทีเซียนอนันต์n1พีn{\displaystyle \prod _{n\geq 1}{\คณิตศาสตร์ {P}__{n}}.

คุณค่าของกระบวนการ ณ เวลาn{\displaystyle n}เป็นพาร์ติชันบีn{\displaystyle B_{n}}ของชุด[n]{\displaystyle [n]}ซึ่งการแจกแจงความน่าจะเป็นถูกกำหนดดังต่อไปนี้ ณ เวลาn=1{\displaystyle n=1}การแบ่งส่วนแบบไม่สำคัญบี1={{1}}{\displaystyle B_{1}=\{\{1\}\}}ได้รับ (ด้วยความน่าจะเป็นหนึ่ง) ณ เวลาn+1{\displaystyle n+1}องค์ประกอบ "n+1{\displaystyle n+1}" คืออย่างใดอย่างหนึ่ง:

  1. เพิ่มเข้าไปในบล็อกหนึ่งของพาร์ติชันบีn{\displaystyle B_{n}}โดยแต่ละบล็อกจะถูกเลือกด้วยความน่าจะเป็น||/(n+1){\displaystyle |b|/(n+1)}ที่ไหน||{\displaystyle |b|}คือขนาดของบล็อก (เช่น จำนวนองค์ประกอบ) หรือ
  2. เพิ่มเข้าไปในพาร์ติชั่นบีn{\displaystyle B_{n}}ในฐานะบล็อกเดี่ยวใหม่ ด้วยความน่าจะเป็น1/(n+1){\displaystyle 1/(n+1)}.

การแบ่งกลุ่มแบบสุ่มที่สร้างขึ้นนั้นมีคุณสมบัติพิเศษบางประการ กล่าวคือสามารถสลับเปลี่ยนได้ในแง่ของการติดป้ายกำกับใหม่{1,...,n}{\displaystyle \{1,...,n\}}ไม่เปลี่ยนแปลงการกระจายของส่วนแบ่ง และมีความสอดคล้องกันในแง่ที่ว่ากฎของการแบ่งส่วนของ[n1]{\displaystyle [n-1]}ได้มาจากการลบองค์ประกอบออกn{\displaystyle n}จากการแบ่งแบบสุ่มบีn{\displaystyle B_{n}}ก็เหมือนกับกฎของการแบ่งแบบสุ่มนั่นเองบีn1{\displaystyle B_{n-1}}.

ความน่าจะเป็นที่กำหนดให้กับที่นั่งแต่ละที่ (โดยไม่คำนึงถึงลำดับการนั่งของลูกค้าแต่ละโต๊ะ) คือ

ปร.(บีn=บี)=บี(||1)!n!,บีพีn{\displaystyle \Pr(B_{n}=B)={\frac {\prod _{b\in B}(|b|-1)!}{n!}},\qquad B\in {\mathcal {P}}_{n}}

ที่ไหน{\displaystyle b}เป็นบล็อกในพาร์ติชันบี{\displaystyle B}และ||{\displaystyle |b|}มีขนาดเท่ากับ{\displaystyle b}.

สามารถขยายความนิยามได้โดยการเพิ่มพารามิเตอร์เข้าไปθ>0{\displaystyle \theta >0}ซึ่งจะปรับเปลี่ยนความน่าจะเป็นที่ลูกค้ารายใหม่จะนั่งที่โต๊ะใหม่เป็นθn+θ{\displaystyle {\frac {\theta }{n+\theta }}}และปรับเปลี่ยนความน่าจะเป็นที่พวกเขาจะนั่งที่โต๊ะขนาดดังกล่าวตามไปด้วย||{\displaystyle |b|}ถึง||n+θ{\displaystyle {\frac {|b|}{n+\theta }}}กระบวนการวานิลลาที่กล่าวถึงข้างต้นสามารถกู้คืนได้โดยการตั้งค่าθ=1{\displaystyle \theta =1}โดยสัญชาตญาณแล้วθ{\displaystyle \theta }สามารถตีความได้ว่าเป็นจำนวนลูกค้าที่นั่งอยู่ที่โต๊ะว่างโต๊ะแรกอย่างแท้จริง

คำจำกัดความทางเลือก

อีกวิธีหนึ่งที่เทียบเท่ากัน แต่แตกต่างกันเล็กน้อยในการกำหนดกระบวนการร้านอาหารจีน คือการให้ลูกค้าใหม่เลือกเพื่อนร่วมโต๊ะแทนที่จะเลือกโต๊ะ[ 4 ]ลูกค้าn+1{\displaystyle n+1}เลือกที่จะนั่งที่โต๊ะเดียวกันกับคนใดคนหนึ่งในกลุ่มนั้นn{\displaystyle n}ลูกค้าที่นั่งด้วยความน่าจะเป็น1n+θ{\displaystyle {\frac {1}{n+\theta }}}หรือเลือกที่จะนั่งที่โต๊ะใหม่ที่ว่างอยู่ด้วยความน่าจะเป็นθn+θ{\displaystyle {\frac {\theta }{n+\theta }}}โปรดสังเกตว่าในรูปแบบนี้ ลูกค้าเลือกโต๊ะโดยไม่ต้องนับจำนวนโต๊ะที่ว่างอยู่—เราไม่จำเป็นต้อง||{\displaystyle |b|}.

การกระจายจำนวนโต๊ะ

การกระจายโต๊ะร้านอาหารจีน ( CRT ) คือการกระจายความน่าจะเป็นของจำนวนโต๊ะในกระบวนการร้านอาหารจีน[ 5 ]สามารถเข้าใจได้ว่าเป็นผลรวมของn{\displaystyle n}ตัวแปร สุ่มเบอร์นูลลีอิสระแต่ละตัวมีพารามิเตอร์ต่างกัน:

เค=ฉัน=1nฉันฉัน~เบอร์นูลลี(θฉัน1+θ){\displaystyle {\begin{aligned}K&=\sum _{i=1}^{n}b_{i}\\[4pt]b_{i}&\sim \operatorname {Bernoulli} \left({\frac {\theta }{i-1+\theta }}\right)\end{aligned}}}

ฟังก์ชันความน่าจะเป็นมวลของเค{\displaystyle K}กำหนดโดย[ 6 ]

เอฟ(เค)=Γ(θ)Γ(n+θ)|(n,เค)|θเค,เค=1,,n,{\displaystyle f(k)={\frac {\Gamma (\theta )}{\Gamma (n+\theta )}}|s(n,k)|\theta ^{k},\quad k=1,\dots ,n,}

ที่ไหน{\displaystyle s}หมายถึง จำนวนสเตอร์ลิง ชนิดแรก

การสรุปทั่วไปแบบสองพารามิเตอร์

โครงสร้างนี้สามารถนำไปปรับใช้กับแบบจำลองที่มีพารามิเตอร์สองตัวได้θ{\displaystyle \theta }&α{\displaystyle \alpha }[ 2 ] [ 7 ]โดยทั่วไปเรียกว่า พารามิเตอร์ ความแข็งแกร่ง (หรือความเข้มข้น)และส่วนลดตามลำดับ ณ เวลาn+1{\displaystyle n+1}ลูกค้ารายต่อไปที่มาถึงพบว่า...|บี|{\displaystyle |B|}โต๊ะที่มีคนนั่งอยู่แล้ว และตัดสินใจนั่งที่โต๊ะว่างด้วยความน่าจะเป็น

θ+|บี|αn+θ,{\displaystyle {\frac {\theta +|B|\alpha }{n+\theta }},}

หรือที่โต๊ะที่มีคนนั่งอยู่{\displaystyle b}ขนาด||{\displaystyle |b|}ด้วยความน่าจะเป็น

||αn+θ.{\displaystyle {\frac {|b|-\alpha }{n+\theta }}.}

เพื่อให้โครงสร้างดังกล่าวสามารถกำหนดมาตรวัดความน่าจะเป็น ที่ถูกต้องได้ จำเป็นต้องสมมติว่าข้อใดข้อหนึ่งต่อไปนี้เป็น จริงα<0{\displaystyle \alpha <0}และθ=แอลα{\displaystyle \theta =-L\alpha }สำหรับบางคนแอล{1,2,,...}{\displaystyle L\in \{1,2,,...\}}หรือที่0α<1{\displaystyle 0\leq \alpha <1}และθ>α{\displaystyle \theta >-\alpha }.

ภายใต้แบบจำลองนี้ ความน่าจะเป็นที่กำหนดให้กับพาร์ติชันใดพาร์ติชันหนึ่งโดยเฉพาะบี{\displaystyle B}ของ[n]{\displaystyle [n]}สามารถแสดงได้ในกรณีทั่วไป (สำหรับค่าใดๆ ของθ,α{\displaystyle \theta ,\alpha }ที่ตรงตามข้อจำกัดที่กล่าวไว้ข้างต้น) ในแง่ของสัญลักษณ์ k ของ Pochhammerดังนี้

ปร.(บีn=บีθ,α)=(θ+α)|บี|1,α(θ+1)n1,1บี(1α)||1,1{\displaystyle \Pr(B_{n}=B\mid \theta ,\alpha )={\frac {(\theta +\alpha )_{|B|-1,\alpha }}{(\theta +1)_{n-1,1}}}\prod _{b\in B}(1-\alpha )_{|b|-1,1}}

โดยที่สัญลักษณ์ k ของ Pochhammer ถูกกำหนดไว้ดังนี้: ตามธรรมเนียมปฏิบัติ(เอ)0,เค=1{\displaystyle (a)_{0,k}=1}และสำหรับ>0{\displaystyle m>0}

(เอ),เค=ฉัน=01(เอ+ฉันเค)={เอถ้า เค=0,เค(เอเค)¯ถ้า เค>0,|เค|(เอ|เค|)_ถ้า เค<0{\displaystyle (a)_{m,k}=\prod _{i=0}^{m-1}(a+ik)={\begin{cases}a^{m}&{\text{if }}k=0,\\\\k^{m}\,({\frac {a}{k}})^{\overline {m}}&{\text{if }}k>0,\\\\\left|k\right|^{m}\,({\frac {a}{\left|k\right|}})^{\underline {m}}&{\text{if }}k<0\end{cases}}}

ที่ไหนx¯=ฉัน=01(x+ฉัน){\displaystyle x^{\overline {m}}=\prod _{i=0}^{m-1}(x+i)}คือแฟกทอเรียลที่เพิ่มขึ้นและx_=ฉัน=01(xฉัน){\displaystyle x^{\underline {m}}=\prod _{i=0}^{m-1}(x-i)}คือแฟกทอเรียลแบบลดลง เป็นที่น่าสังเกตว่าสำหรับการตั้งค่าพารามิเตอร์ที่α<0{\displaystyle \alpha <0}และθ=แอลα{\displaystyle \theta =-L\alpha }, แล้ว(θ+α)|บี|1,α=(|α|(แอล1))|บี|1,α{\displaystyle (\theta +\alpha )_{|B|-1,\alpha }=(|\alpha |(L-1))_{|B|-1,\alpha }}ซึ่งจะมีค่าเป็นศูนย์เมื่อใดก็ตามที่|บี|>แอล{\displaystyle |B|>L}ดังนั้นแอล{\displaystyle L}เป็นค่าขอบเขตบนของจำนวนบล็อกในพาร์ติชันโปรดดูรายละเอียดเพิ่มเติม ในหัวข้อย่อยเกี่ยวกับ แบบจำลอง Dirichlet-categorical ด้านล่าง

สำหรับกรณีที่θ>0{\displaystyle \theta >0}และ0<α<1{\displaystyle 0<\alpha <1}ความน่าจะเป็นของการแบ่งกลุ่มสามารถเขียนใหม่ได้ในรูปของฟังก์ชันแกมมาดังนี้

ปร.(บีn=บีθ,α)=Γ(θ)Γ(θ+n)α|บี|Γ(θ/α+|บี|)Γ(θ/α)บีΓ(||α)Γ(1α).{\displaystyle \Pr(B_{n}=B\mid \theta ,\alpha )={\frac {\Gamma (\theta )}{\Gamma (\theta +n)}}{\dfrac {\alpha ^{|B|}\,\Gamma (\theta /\alpha +|B|)}{\Gamma (\theta /\alpha )}}\prod _{b\in B}{\dfrac {\Gamma (|b|-\alpha )}{\Gamma (1-\alpha )}}.}

ในกรณีที่มีพารามิเตอร์เดียว โดยที่α{\displaystyle \alpha }มีค่าเป็นศูนย์ และθ>0{\displaystyle \theta >0}สิ่งนี้ทำให้ง่ายขึ้นเป็น

ปร.(บีn=บีθ)=Γ(θ)θ|บี|Γ(θ+n)บีΓ(||).{\displaystyle \Pr(B_{n}=B\mid \theta )={\frac {\Gamma (\theta )\,\theta ^{|B|}}{\Gamma (\theta +n)}}\prod _{b\in B}\Gamma (|b|).}

หรือเมื่อθ{\displaystyle \theta }มีค่าเป็นศูนย์ และ0<α<1{\displaystyle 0<\alpha <1}

ปร.(บีn=บีα)=α|บี|1Γ(|บี|)Γ(n)บีΓ(||α)Γ(1α).{\displaystyle \Pr(B_{n}=B\mid \alpha )={\frac {\alpha ^{|B|-1}\,\Gamma (|B|)}{\Gamma (n)}}\prod _{b\in B}{\frac {\Gamma (|b|-\alpha )}{\Gamma (1-\alpha )}}.}

เช่นเดียวกับที่ผ่านมา ความน่าจะเป็นที่กำหนดให้กับพาร์ติชันใดๆ ขึ้นอยู่กับขนาดของบล็อกเท่านั้น ดังนั้นพาร์ติชันแบบสุ่มจึงสามารถสลับเปลี่ยนได้ในความหมายที่อธิบายไว้ข้างต้น คุณสมบัติความสอดคล้องยังคงใช้ได้เช่นเดิม โดยอาศัยโครงสร้าง

ถ้าα=0{\displaystyle \alpha =0}การแจกแจงความน่าจะเป็นของการแบ่งส่วนแบบสุ่มของจำนวนเต็มn{\displaystyle n}ดังนั้นจึงเกิดการแจกแจงแบบ Ewensที่มีพารามิเตอร์θ{\displaystyle \theta }ใช้ในพันธุศาสตร์ประชากรและทฤษฎีความเป็นกลางแบบเอกภาพของความหลากหลายทางชีวภาพ

ภาพเคลื่อนไหวแสดงขั้นตอนการทำงานในร้านอาหารจีน พร้อมพารามิเตอร์การปรับขนาดθ=0.5, α=0{\displaystyle \theta =0.5,\ \alpha =0}โต๊ะจะถูกซ่อนเมื่อไม่สามารถแสดงลูกค้าของโต๊ะได้อีกต่อไป อย่างไรก็ตาม โต๊ะทุกโต๊ะมีที่นั่งไม่จำกัดจำนวน (บันทึกภาพเคลื่อนไหวแบบโต้ตอบ[ 8 ] )

อนุพันธ์

นี่คือวิธีหนึ่งในการหาความน่าจะเป็นของการแบ่งกลุ่มนี้ ให้ซีฉัน{\displaystyle C_{i}}เป็นบล็อกสุ่มที่หมายเลขนั้นเข้าไปฉัน{\displaystyle i}มีการเพิ่มสำหรับฉัน=1,2,3,...{\displaystyle i=1,2,3,...}. แล้ว

ปร.(ซีฉัน=ซี1,,ซีฉัน1)={θ+|บี|αθ+ฉัน1ถ้า บล็อกใหม่,||αθ+ฉัน1ถ้า ;{\displaystyle \Pr(C_{i}=c\mid C_{1},\ldots ,C_{i-1})={\begin{cases}{\dfrac {\theta +|B|\alpha }{\theta +i-1}}&{\text{if }}c\in {\text{new block}},\\\\{\dfrac {|b|-\alpha }{\theta +i-1}}&{\text{if }}c\in b;\end{cases}}}

ความน่าจะเป็นที่บีn{\displaystyle B_{n}}คือการแบ่งส่วนเฉพาะใดๆ ของเซต{1,...,n}{\displaystyle \{1,...,n\}}ผลคูณของความน่าจะเป็นเหล่านี้คือฉัน{\displaystyle i}วิ่งจาก1{\displaystyle 1}ถึงn{\displaystyle n}ทีนี้ลองพิจารณาขนาดของบล็อกดู{\displaystyle b}: มันจะเพิ่มขึ้นทีละหนึ่งทุกครั้งที่เราเพิ่มองค์ประกอบเข้าไป เมื่อองค์ประกอบสุดท้ายในบล็อก{\displaystyle b}จะต้องเพิ่มเข้าไป โดยขนาดบล็อกคือ||1{\displaystyle |b|-1}ตัวอย่างเช่น ลองพิจารณาลำดับของตัวเลือกต่อไปนี้: (สร้างบล็อกใหม่){\displaystyle b})(เข้าร่วม{\displaystyle b})(เข้าร่วม{\displaystyle b})(เข้าร่วม{\displaystyle b}). ในที่สุด บล็อก{\displaystyle b}มี 4 องค์ประกอบ และผลคูณของตัวเศษในสมการข้างต้นคือθ123{\displaystyle \theta \cdot 1\cdot 2\cdot 3}เมื่อใช้ตรรกะนี้ เราจะได้ปร.(บีn=บี){\displaystyle \Pr(B_{n}=B)}ตามที่กล่าวมาข้างต้น

จำนวนโต๊ะที่คาดหวัง

สำหรับกรณีที่มีพารามิเตอร์เดียวα=0{\displaystyle \alpha =0}และ0<θ<{\displaystyle 0<\theta <\infty }จำนวนโต๊ะจะถูกจัดสรรตามรูปแบบการจัดสรรโต๊ะในร้านอาหารจีนค่าที่คาดหวังของตัวแปรสุ่มนี้ โดยกำหนดให้มีจำนวนโต๊ะเท่ากับจำนวนโต๊ะทั้งหมดn{\displaystyle n}ลูกค้าที่นั่งคือ[ 9 ]

เค=1nθθ+เค1=θ(Ψ(θ+n)Ψ(θ)){\displaystyle {\begin{aligned}\sum _{k=1}^{n}{\frac {\theta }{\theta +k-1}}=\theta \cdot (\Psi (\theta +n)-\Psi (\theta ))\end{aligned}}}

ที่ไหนΨ(θ){\displaystyle \Psi (\theta )}คือฟังก์ชันไดแกมมาสำหรับกรณีสองพารามิเตอร์ สำหรับα0{\displaystyle \alpha \neq 0}จำนวนโต๊ะที่ถูกจองโดยประมาณคือ[ 7 ]

(θ+α)n¯α(θ+1)n1¯θα,{\displaystyle {\begin{aligned}{\frac {(\theta +\alpha )^{\overline {n}}}{\alpha (\theta +1)^{\overline {n-1}}}}-{\frac {\theta }{\alpha }},\end{aligned}}}

ที่ไหนx¯{\displaystyle x^{\overline {m}}}คือแฟกทอเรียลที่เพิ่มขึ้น (ตามที่นิยามไว้ข้างต้น)

แบบจำลอง Dirichlet-categorical

สำหรับการเลือกพารามิเตอร์α<0{\displaystyle \alpha <0}และθ=แอลα{\displaystyle \theta =-L\alpha }, ที่ไหนแอล{1,2,3,}{\displaystyle L\in \{1,2,3,\ldots \}}กระบวนการร้านอาหารจีนแบบสองพารามิเตอร์นั้นเทียบเท่ากับแบบจำลอง Dirichlet-categoricalซึ่งเป็นแบบจำลองลำดับชั้นที่สามารถกำหนดได้ดังนี้ โปรดสังเกตว่าสำหรับการตั้งค่าพารามิเตอร์นี้ ความน่าจะเป็นของการเข้าใช้โต๊ะใหม่เมื่อมีโต๊ะว่างอยู่แล้วแอล{\displaystyle L}จำนวนโต๊ะที่มีคนนั่งอยู่เป็นศูนย์ ดังนั้นจำนวนโต๊ะที่มีคนนั่งอยู่จึงมีค่าสูงสุดไม่เกินแอล{\displaystyle L}หากเราเลือกที่จะระบุตารางด้วยป้ายกำกับที่รับค่าใน{1,2,,แอล}{\displaystyle \{1,2,\ldots ,L\}}จากนั้นจึงสร้างการแบ่งแบบสุ่มของชุดข้อมูล[n]={1,2,,n}{\displaystyle [n]=\{1,2,\ldots ,n\}}แบบจำลองลำดับชั้นจะสร้างการกระจายป้ายกำกับตามหมวดหมู่ ขึ้นมาก่อนพี=(พี1,พี2,,พีแอล){\displaystyle \mathbf {p} =(p_{1},p_{2},\ldots ,p_{L})}จากการกระจายแบบ Dirichlet สมมาตร โดยมีพารามิเตอร์ความเข้มข้นγ=α>0{\displaystyle \gamma =-\alpha >0}จากนั้น ดำเนินการแยกกันสำหรับแต่ละส่วนn{\displaystyle n}ลูกค้า ป้ายกำกับตารางถูกดึงมาจากหมวดหมู่พี{\displaystyle \mathbf {p} }เนื่องจากการแจกแจงแบบ Dirichlet เป็นการแจกแจงแบบคู่ควบกับการแจกแจงแบบจัดกลุ่ม ดังนั้นตัวแปรที่ซ่อนอยู่พี{\displaystyle \mathbf {p} }สามารถแยกส่วนย่อยออกเพื่อให้ได้การแจกแจงการทำนายภายหลังสำหรับสถานะป้ายกำกับถัดไป ได้n+1{\displaystyle \ell _{n+1}}, ที่ให้ไว้n{\displaystyle n}ป้ายกำกับก่อนหน้า

พี(n+1=ฉัน1,,n)=γ+|ฉัน|แอลγ+n{\displaystyle P(\ell _{n+1}=i\mid \ell _{1},\ldots ,\ell _{n})={\frac {\gamma +\left|{b_{i}}\right|}{L\gamma +n}}}

ที่ไหน|ฉัน|0{\displaystyle \left|{b_{i}}\right|\geq 0}คือจำนวนลูกค้าที่นั่งอยู่ที่โต๊ะแล้วฉัน{\displaystyle i}. กับα=γ{\displaystyle \alpha =-\gamma }และθ=แอลγ{\displaystyle \theta =L\gamma }ซึ่งสอดคล้องกับสูตรทั่วไปข้างต้น|ฉัน|αn+θ{\displaystyle {\frac {|b_{i}|-\alpha }{n+\theta }}}สำหรับความน่าจะเป็นของการนั่งที่โต๊ะที่มีคนนั่งอยู่ เมื่อ|ฉัน|1{\displaystyle |b_{i}|\geq 1}ความน่าจะเป็นที่จะนั่งที่ใดที่หนึ่งแอล|บี|{\displaystyle L-|B|}โต๊ะที่ว่างอยู่ก็สอดคล้องกับสูตรทั่วไปและกำหนดโดย

ฉัน:|ฉัน|=0พี(n+1=ฉัน1,,n)=(แอล|บี|)γn+แอลγ=θ+|บี|αn+θ{\displaystyle \sum _{i:|b_{i}|=0}P(\ell _{n+1}=i\mid \ell _{1},\ldots ,\ell _{n})={\frac {(L-|B|)\gamma }{n+L\gamma }}={\frac {\theta +|B|\alpha }{n+\theta }}}

ความน่าจะเป็นแบบมาร์จินัลสำหรับป้ายกำกับต่างๆ กำหนดโดย

พี(1,,n)=พี(1)ที=1n1พี(ที+11,,ที)=ฉัน=1แอลγ|ฉัน|¯(แอลγ)n¯{\displaystyle P(\ell _{1},\ldots ,\ell _{n})=P(\ell _{1})\prod _{t=1}^{n-1}P(\ell _{t+1}\mid \ell _{1},\ldots ,\ell _{t})={\frac {\prod _{i=1}^{L}\gamma ^{\overline {\left|{b_{i}}\right|}}}{(L\gamma )^{\overline {n}}}}}

ที่ไหนพี(1)=1แอล{\displaystyle P(\ell _{1})={\frac {1}{L}}}และx¯=ฉัน=01(x+ฉัน){\displaystyle x^{\overline {m}}=\prod _{i=0}^{m-1}(x+i)}คือแฟกทอเรียลที่เพิ่มขึ้นโดยทั่วไปแล้วจะมีสถานะป้ายกำกับหลายสถานะที่สอดคล้องกับ พาร์ติชัน เดียวกันสำหรับพาร์ติชันที่กำหนดบี{\displaystyle B}ซึ่งมี|บี|แอล{\displaystyle \left|B\right|\leq L}บล็อก จำนวนสถานะป้ายกำกับที่สอดคล้องกับพาร์ติชันนี้ทั้งหมด กำหนดโดย แฟกทอเรี ยลแบบลดลงแอล|บี|_=ฉัน=0|บี|1(แอลฉัน){\displaystyle L^{\underline {\left|B\right|}}=\prod _{i=0}^{\left|B\right|-1}(L-i)}เมื่อพิจารณาถึงเรื่องนี้แล้ว ความน่าจะเป็นสำหรับการแบ่งกลุ่มคือ

ปร.(บีn=บีγ,แอล)=แอล|บี|_ฉัน=1แอลγ|ฉัน|¯(แอลγ)n¯{\displaystyle {\text{Pr}}(B_{n}=B\mid \gamma ,L)=L^{\underline {\left|B\right|}}\,{\frac {\prod _{i=1}^{L}\gamma ^{\overline {\left|{b_{i}}\right|}}}{(L\gamma )^{\overline {n}}}}}

ซึ่งสามารถตรวจสอบได้ว่าสอดคล้องกับความน่าจะเป็นของการแบ่งส่วนทั่วไปที่ระบุไว้ข้างต้นในรูปของสัญลักษณ์ k ของ Pochhammer โปรดสังเกตอีกครั้งว่า ถ้าบี{\displaystyle B}อยู่นอกเหนือการสนับสนุน เช่น|บี|>แอล{\displaystyle |B|>L}แฟกทอเรียลที่ลดลงแอล|บี|_{\displaystyle L^{\underline {|B|}}}มีค่าเป็นศูนย์ตามที่ควรจะเป็น (การใช้งานจริงที่ประเมินความน่าจะเป็นของลอการิทึมสำหรับพาร์ติชันผ่าน บันทึกแอล|บี|_=บันทึก|Γ(แอล+1)|บันทึก|Γ(แอล+1|บี|)|{\displaystyle \log L^{\underline {|B|}}=\log \left|\Gamma (L+1)\right|-\log \left|\Gamma (L+1-|B|)\right|}จะกลับมา{\displaystyle -\infty }เมื่อใดก็ตามที่|บี|>แอล{\displaystyle |B|>L}(ตามความจำเป็น)

ความสัมพันธ์ระหว่าง CRP แบบ Dirichlet-categorical และ CRP แบบพารามิเตอร์เดียว

ในอีกด้านหนึ่ง ลองพิจารณาขั้นตอนการดำเนินงานร้านอาหารจีนที่มีพารามิเตอร์เดียวα=0{\displaystyle \alpha =0}และθ>0{\displaystyle \theta >0}ซึ่งเราใช้สัญลักษณ์แทนซีอาร์พี(α=0,θ){\displaystyle {\text{CRP}}(\alpha =0,\theta )}และในอีกด้านหนึ่ง แบบจำลอง Dirichlet-categorical ที่มีแอล{\displaystyle L}จำนวนเต็มบวกและตำแหน่งที่เราเลือกγ=θแอล{\displaystyle \gamma ={\frac {\theta }{L}}}ซึ่งดังที่แสดงไว้ข้างต้น เทียบเท่ากับซีอาร์พี(α=θแอล,θ){\displaystyle {\text{CRP}}(\alpha =-{\frac {\theta }{L}},\theta )}สิ่งนี้แสดงให้เห็นว่าแบบจำลอง Dirichlet-categorical สามารถทำให้ใกล้เคียงกับค่าที่กำหนดได้มากเท่าใดก็ได้ซีอาร์พี(0,θ){\displaystyle {\text{CRP}}(0,\theta )}โดยการสร้างแอล{\displaystyle L}ใหญ่.

กระบวนการหักแท่ง

กระบวนการร้านอาหารจีนสองพารามิเตอร์สามารถกำหนดได้อย่างเทียบเท่าในแง่ของกระบวนการหักแท่ง [ 10 ] สำหรับกรณีที่0α<1{\displaystyle 0\leq \alpha <1}และθ>α{\displaystyle \theta >-\alpha }กระบวนการหักแท่งสามารถอธิบายได้ว่าเป็นแบบจำลองลำดับชั้น คล้ายกับแบบจำลอง Dirichlet-categorical ข้างต้น ยกเว้นว่ามีสถานะป้ายกำกับจำนวนอนันต์ ป้ายกำกับตารางถูกสุ่มอย่างอิสระจากการกระจายเชิงหมวดหมู่ที่ไม่มีที่สิ้นสุดพี=(พี1,พี2,){\displaystyle \mathbf {p} =(p_{1},p_{2},\ldots )}โดยส่วนประกอบต่างๆ จะถูกสุ่มตัวอย่างโดยใช้การหักแท่งไม้ : เริ่มต้นด้วยแท่งไม้ที่มีความยาว 1 และหักมันออกเป็นสองส่วนโดยสุ่ม ความยาวของครึ่งซ้ายคือพี1{\displaystyle p_{1}}และครึ่งขวาถูกแบ่งออกอีกครั้งแบบวนซ้ำเพื่อให้ได้พี2,พี3,{\displaystyle p_{2},p_{3},\ldots }กล่าวให้แม่นยำยิ่งขึ้นคือ เศษส่วนทางซ้ายเอฟเค{\displaystyle f_{k}}ของเค{\displaystyle k}-th break ถูกสุ่มมาจากการแจกแจงแบบเบต้า :

เอฟเค~บี(1α,θ+เคα),สำหรับ เค1 และ 0α<1{\displaystyle f_{k}\sim B(1-\alpha ,\theta +k\alpha ),\;{\text{for }}k\geq 1{\text{ and }}0\leq \alpha <1}

ความน่าจะเป็นเชิงหมวดหมู่มีดังนี้:

พีเค=เอฟเคฉัน=1เค1(1เอฟเค),โดยที่ผลคูณว่างเปล่าจะมีค่าเท่ากับหนึ่ง{\displaystyle p_{k}=f_{k}\prod _{i=1}^{k-1}(1-f_{k}),\;{\text{where the empty product evaluates to one.}}}

สำหรับการตั้งค่าพารามิเตอร์α<0{\displaystyle \alpha <0}และθ=αแอล{\displaystyle \theta =-\alpha L}, ที่ไหนแอล{\displaystyle L}เป็นจำนวนเต็มบวก และในกรณีที่หมวดหมู่มีค่าจำกัด:พี=(พี1,,พีแอล){\displaystyle \mathbf {p} =(p_{1},\ldots ,p_{L})}เราสามารถสุ่มตัวอย่างได้พี{\displaystyle \mathbf {p} }จากการแจกแจงแบบ Dirchlet ทั่วไปตามที่อธิบายไว้ข้างต้นแต่ยังสามารถสุ่มตัวอย่างได้โดยใช้ สูตรการหักแท่ง แบบตัดทอนโดยที่สูตรสำหรับการสุ่มเศษส่วนจะถูกแก้ไขเป็น:

เอฟเค~บี(α,θ+เคα),สำหรับ 1เคแอล1 และ α<0{\displaystyle f_{k}\sim B(-\alpha ,\theta +k\alpha ),\;{\text{for }}1\leq k\leq L-1{\text{ and }}\alpha <0}

และเอฟแอล=1{\displaystyle f_{L}=1}.

กระบวนการบุฟเฟต์แบบอินเดีย

เป็นไปได้ที่จะปรับโมเดลเพื่อให้แต่ละจุดข้อมูลไม่เกี่ยวข้องกับคลาสใดคลาสหนึ่งโดยเฉพาะอีกต่อไป (กล่าวคือ เราไม่ได้สร้างพาร์ติชันอีกต่อไป) แต่สามารถเกี่ยวข้องกับคลาสใดๆ ก็ได้ ซึ่งทำให้การเปรียบเทียบโต๊ะในร้านอาหารไม่ชัดเจน จึงเปรียบเทียบกับกระบวนการที่ผู้รับประทานอาหารหลายคนเลือกตัวอย่างจากชุดย่อยของอาหารที่มีให้เลือกมากมายไม่จำกัดในบุฟเฟต์ ความน่าจะเป็นที่ผู้รับประทานอาหารคนใดคนหนึ่งจะเลือกอาหารจานใดจานหนึ่งนั้นเป็นสัดส่วนกับความนิยมของอาหารจานนั้นในหมู่ผู้รับประทานอาหารที่ผ่านมา และนอกจากนี้ผู้รับประทานอาหารอาจเลือกจากอาหารที่ยังไม่ได้ลองด้วย กระบวนการนี้ได้รับการตั้งชื่อว่ากระบวนการบุฟเฟต์แบบอินเดีย และสามารถใช้เพื่ออนุมานคุณลักษณะแฝงในข้อมูลได้[ 11 ]

แอปพลิเคชัน

กระบวนการร้านอาหารจีนมีความเชื่อมโยงอย่างใกล้ชิดกับกระบวนการ Dirichletและแผนการโถของ Pólyaดังนั้นจึงมีประโยชน์ในการประยุกต์ใช้สถิติแบบเบย์เซียนรวมถึงวิธีการแบบเบย์เซียนที่ไม่ใช่พาราเมตริก กระบวนการร้านอาหารจีนแบบทั่วไปมีความเกี่ยวข้องอย่างใกล้ชิดกับกระบวนการ Pitman–Yorกระบวนการเหล่านี้ถูกนำไปใช้ในแอปพลิเคชันมากมาย รวมถึงการสร้างแบบจำลองข้อความ การจัดกลุ่มข้อมูลไมโครอาร์เรย์ ทางชีวภาพ [ 12 ]การสร้างแบบจำลองความหลากหลายทางชีวภาพและการสร้างภาพขึ้นใหม่[ 13 ] [ 14 ]

ดูเพิ่มเติม

  • หนังสือแนะนำเกี่ยวกับการกระจายแบบ Dirichlet และกระบวนการที่เกี่ยวข้อง โดย Frigyik, Kapila และ Gupta
  • การบรรยายโดยไมเคิล ไอ. จอร์แดนเกี่ยวกับโครงการ CRP:
    • http://videolectures.net/icml05_jordan_dpcrp/

สรุปเนื้อหา

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

ข้อมูลสำคัญเกี่ยวกับ กระบวนการร้านอาหารจีน

ใน ทฤษฎีความน่าจะเป็น กระบวนการ ร้าน อาหารจีน เป็น กระบวนการสุ่มแบบ เวลาไม่ต่อเนื่อง คล้ายกับการจัดที่นั่งให้ลูกค้าที่โต๊ะในร้านอาหาร ลองนึกภาพร้านอาหารที่มีโต๊ะกลมจำนวนอนันต์...

คำจำกัดความอย่างเป็นทางการ

สำหรับจำนวนเต็มบวกใดๆ n {\displaystyle n} , อนุญาต พี n {\displaystyle {\mathcal {P}}_{n}} แทนเซตของการแบ่งย่อยทั้งหมดของเซต { 1 , 2 , 3 , . . . , n } ≜ [ n ] {\displaystyle \{1,2,3,...

คำจำกัดความทางเลือก

อีกวิธีหนึ่งที่เทียบเท่ากัน แต่แตกต่างกันเล็กน้อยในการกำหนดกระบวนการร้านอาหารจีน คือการให้ลูกค้าใหม่เลือกเพื่อนร่วมโต๊ะแทนที่จะเลือกโต๊ะ [ 4 ] ลูกค้า n + 1 {\displaystyle n+1} เลือกที่จะนั่งที่โต๊ะเดียวกันกับคนใดคนหนึ่งในกลุ่มนั้น n {\displaystyle n}...

การกระจายจำนวนโต๊ะ

การ กระจายโต๊ะร้านอาหารจีน ( CRT ) คือ การกระจายความน่าจะ เป็นของจำนวนโต๊ะในกระบวนการร้านอาหารจีน [ 5 ] สามารถเข้าใจได้ว่าเป็นผลรวมของ n {\displaystyle n} ตัวแปร สุ่ม เบอร์นูลลีอิสระ แต่ละตัวมีพารามิเตอร์ต่างกัน: