กระบวนการร้านอาหารจีน
ในทฤษฎีความน่าจะเป็น กระบวนการ ร้านอาหารจีนเป็นกระบวนการสุ่มแบบเวลาไม่ต่อเนื่อง คล้ายกับการจัดที่นั่งให้ลูกค้าที่โต๊ะในร้านอาหาร ลองนึกภาพร้านอาหารที่มีโต๊ะกลมจำนวนอนันต์ แต่ละโต๊ะมีความจุไม่จำกัด ลูกค้าคนที่ 1 นั่งที่โต๊ะแรก ลูกค้าคนถัดไปอาจนั่งที่โต๊ะเดียวกับลูกค้าคนที่ 1 หรือโต๊ะถัดไป กระบวนการนี้ดำเนินต่อไป โดยลูกค้าแต่ละคนเลือกที่จะนั่งที่โต๊ะที่มีคนนั่งอยู่แล้ว ด้วยความน่าจะเป็นที่แปรผันตามจำนวนลูกค้าที่นั่งอยู่แล้ว (กล่าวคือ พวกเขามีโอกาสนั่งที่โต๊ะที่มีคนเยอะมากกว่าโต๊ะที่มีคนน้อย) หรือโต๊ะว่าง ณ เวลาnลูกค้าnคนจะถูกแบ่งออกเป็นm ≤ nโต๊ะ (หรือกลุ่มของส่วนแบ่ง) ผลลัพธ์ของกระบวนการนี้สามารถสลับเปลี่ยนได้หมายความว่าลำดับที่ลูกค้านั่งไม่มีผลต่อความน่าจะเป็นของการกระจาย สุดท้าย คุณสมบัตินี้ช่วยลดความซับซ้อนของปัญหาหลายอย่างในพันธุศาสตร์ประชากรการวิเคราะห์ทางภาษาศาสตร์และการจดจำภาพได้ อย่างมาก
การเปรียบเทียบกับร้านอาหารปรากฏครั้งแรกในงานเขียนของDavid Aldous ในปี 1985 [ 1 ] ซึ่งอ้างอิงถึงJim Pitman (ซึ่งให้เครดิตLester Dubins เพิ่มเติมด้วย ) [ 2 ]
กระบวนการแบ่งส่วนที่เทียบเท่ากันได้รับการตีพิมพ์เมื่อปีก่อนโดยFred Hoppe [ 3 ]โดยใช้ "แผนผังโถ" ที่คล้ายกับโถของ Pólya เมื่อเปรียบเทียบกับ แบบจำลองโถของ Hoppe กระบวนการร้านอาหารจีนมีข้อได้เปรียบตรงที่สามารถอธิบายการเรียงสับเปลี่ยนแบบสุ่มผ่านโครงสร้างวงจรได้อย่างเป็นธรรมชาติ นอกเหนือจากการอธิบายการแบ่งส่วนแบบสุ่ม
คำจำกัดความอย่างเป็นทางการ
สำหรับจำนวนเต็มบวกใดๆ, อนุญาตแทนเซตของการแบ่งย่อยทั้งหมดของเซตกระบวนการของร้านอาหารจีนนั้นใช้ค่าในผลคูณคาร์ทีเซียนอนันต์.
คุณค่าของกระบวนการ ณ เวลาเป็นพาร์ติชันของชุดซึ่งการแจกแจงความน่าจะเป็นถูกกำหนดดังต่อไปนี้ ณ เวลาการแบ่งส่วนแบบไม่สำคัญได้รับ (ด้วยความน่าจะเป็นหนึ่ง) ณ เวลาองค์ประกอบ "" คืออย่างใดอย่างหนึ่ง:
- เพิ่มเข้าไปในบล็อกหนึ่งของพาร์ติชันโดยแต่ละบล็อกจะถูกเลือกด้วยความน่าจะเป็นที่ไหนคือขนาดของบล็อก (เช่น จำนวนองค์ประกอบ) หรือ
- เพิ่มเข้าไปในพาร์ติชั่นในฐานะบล็อกเดี่ยวใหม่ ด้วยความน่าจะเป็น.
การแบ่งกลุ่มแบบสุ่มที่สร้างขึ้นนั้นมีคุณสมบัติพิเศษบางประการ กล่าวคือสามารถสลับเปลี่ยนได้ในแง่ของการติดป้ายกำกับใหม่ไม่เปลี่ยนแปลงการกระจายของส่วนแบ่ง และมีความสอดคล้องกันในแง่ที่ว่ากฎของการแบ่งส่วนของได้มาจากการลบองค์ประกอบออกจากการแบ่งแบบสุ่มก็เหมือนกับกฎของการแบ่งแบบสุ่มนั่นเอง.
ความน่าจะเป็นที่กำหนดให้กับที่นั่งแต่ละที่ (โดยไม่คำนึงถึงลำดับการนั่งของลูกค้าแต่ละโต๊ะ) คือ
ที่ไหนเป็นบล็อกในพาร์ติชันและมีขนาดเท่ากับ.
สามารถขยายความนิยามได้โดยการเพิ่มพารามิเตอร์เข้าไปซึ่งจะปรับเปลี่ยนความน่าจะเป็นที่ลูกค้ารายใหม่จะนั่งที่โต๊ะใหม่เป็นและปรับเปลี่ยนความน่าจะเป็นที่พวกเขาจะนั่งที่โต๊ะขนาดดังกล่าวตามไปด้วยถึงกระบวนการวานิลลาที่กล่าวถึงข้างต้นสามารถกู้คืนได้โดยการตั้งค่าโดยสัญชาตญาณแล้วสามารถตีความได้ว่าเป็นจำนวนลูกค้าที่นั่งอยู่ที่โต๊ะว่างโต๊ะแรกอย่างแท้จริง
คำจำกัดความทางเลือก
อีกวิธีหนึ่งที่เทียบเท่ากัน แต่แตกต่างกันเล็กน้อยในการกำหนดกระบวนการร้านอาหารจีน คือการให้ลูกค้าใหม่เลือกเพื่อนร่วมโต๊ะแทนที่จะเลือกโต๊ะ[ 4 ]ลูกค้าเลือกที่จะนั่งที่โต๊ะเดียวกันกับคนใดคนหนึ่งในกลุ่มนั้นลูกค้าที่นั่งด้วยความน่าจะเป็นหรือเลือกที่จะนั่งที่โต๊ะใหม่ที่ว่างอยู่ด้วยความน่าจะเป็นโปรดสังเกตว่าในรูปแบบนี้ ลูกค้าเลือกโต๊ะโดยไม่ต้องนับจำนวนโต๊ะที่ว่างอยู่—เราไม่จำเป็นต้อง.
การกระจายจำนวนโต๊ะ
การกระจายโต๊ะร้านอาหารจีน ( CRT ) คือการกระจายความน่าจะเป็นของจำนวนโต๊ะในกระบวนการร้านอาหารจีน[ 5 ]สามารถเข้าใจได้ว่าเป็นผลรวมของตัวแปร สุ่มเบอร์นูลลีอิสระแต่ละตัวมีพารามิเตอร์ต่างกัน:
ฟังก์ชันความน่าจะเป็นมวลของกำหนดโดย[ 6 ]
ที่ไหนหมายถึง จำนวนสเตอร์ลิง ชนิดแรก
การสรุปทั่วไปแบบสองพารามิเตอร์
โครงสร้างนี้สามารถนำไปปรับใช้กับแบบจำลองที่มีพารามิเตอร์สองตัวได้&[ 2 ] [ 7 ]โดยทั่วไปเรียกว่า พารามิเตอร์ ความแข็งแกร่ง (หรือความเข้มข้น)และส่วนลดตามลำดับ ณ เวลาลูกค้ารายต่อไปที่มาถึงพบว่า...โต๊ะที่มีคนนั่งอยู่แล้ว และตัดสินใจนั่งที่โต๊ะว่างด้วยความน่าจะเป็น
หรือที่โต๊ะที่มีคนนั่งอยู่ขนาดด้วยความน่าจะเป็น
เพื่อให้โครงสร้างดังกล่าวสามารถกำหนดมาตรวัดความน่าจะเป็น ที่ถูกต้องได้ จำเป็นต้องสมมติว่าข้อใดข้อหนึ่งต่อไปนี้เป็น จริงและสำหรับบางคนหรือที่และ.
ภายใต้แบบจำลองนี้ ความน่าจะเป็นที่กำหนดให้กับพาร์ติชันใดพาร์ติชันหนึ่งโดยเฉพาะของสามารถแสดงได้ในกรณีทั่วไป (สำหรับค่าใดๆ ของที่ตรงตามข้อจำกัดที่กล่าวไว้ข้างต้น) ในแง่ของสัญลักษณ์ k ของ Pochhammerดังนี้
โดยที่สัญลักษณ์ k ของ Pochhammer ถูกกำหนดไว้ดังนี้: ตามธรรมเนียมปฏิบัติและสำหรับ
ที่ไหนคือแฟกทอเรียลที่เพิ่มขึ้นและคือแฟกทอเรียลแบบลดลง เป็นที่น่าสังเกตว่าสำหรับการตั้งค่าพารามิเตอร์ที่และ, แล้วซึ่งจะมีค่าเป็นศูนย์เมื่อใดก็ตามที่ดังนั้นเป็นค่าขอบเขตบนของจำนวนบล็อกในพาร์ติชันโปรดดูรายละเอียดเพิ่มเติม ในหัวข้อย่อยเกี่ยวกับ แบบจำลอง Dirichlet-categorical ด้านล่าง
สำหรับกรณีที่และความน่าจะเป็นของการแบ่งกลุ่มสามารถเขียนใหม่ได้ในรูปของฟังก์ชันแกมมาดังนี้
ในกรณีที่มีพารามิเตอร์เดียว โดยที่มีค่าเป็นศูนย์ และสิ่งนี้ทำให้ง่ายขึ้นเป็น
หรือเมื่อมีค่าเป็นศูนย์ และ
เช่นเดียวกับที่ผ่านมา ความน่าจะเป็นที่กำหนดให้กับพาร์ติชันใดๆ ขึ้นอยู่กับขนาดของบล็อกเท่านั้น ดังนั้นพาร์ติชันแบบสุ่มจึงสามารถสลับเปลี่ยนได้ในความหมายที่อธิบายไว้ข้างต้น คุณสมบัติความสอดคล้องยังคงใช้ได้เช่นเดิม โดยอาศัยโครงสร้าง
ถ้าการแจกแจงความน่าจะเป็นของการแบ่งส่วนแบบสุ่มของจำนวนเต็มดังนั้นจึงเกิดการแจกแจงแบบ Ewensที่มีพารามิเตอร์ใช้ในพันธุศาสตร์ประชากรและทฤษฎีความเป็นกลางแบบเอกภาพของความหลากหลายทางชีวภาพ
อนุพันธ์
นี่คือวิธีหนึ่งในการหาความน่าจะเป็นของการแบ่งกลุ่มนี้ ให้เป็นบล็อกสุ่มที่หมายเลขนั้นเข้าไปมีการเพิ่มสำหรับ. แล้ว
ความน่าจะเป็นที่คือการแบ่งส่วนเฉพาะใดๆ ของเซตผลคูณของความน่าจะเป็นเหล่านี้คือวิ่งจากถึงทีนี้ลองพิจารณาขนาดของบล็อกดู: มันจะเพิ่มขึ้นทีละหนึ่งทุกครั้งที่เราเพิ่มองค์ประกอบเข้าไป เมื่อองค์ประกอบสุดท้ายในบล็อกจะต้องเพิ่มเข้าไป โดยขนาดบล็อกคือตัวอย่างเช่น ลองพิจารณาลำดับของตัวเลือกต่อไปนี้: (สร้างบล็อกใหม่))(เข้าร่วม)(เข้าร่วม)(เข้าร่วม). ในที่สุด บล็อกมี 4 องค์ประกอบ และผลคูณของตัวเศษในสมการข้างต้นคือเมื่อใช้ตรรกะนี้ เราจะได้ตามที่กล่าวมาข้างต้น
จำนวนโต๊ะที่คาดหวัง
สำหรับกรณีที่มีพารามิเตอร์เดียวและจำนวนโต๊ะจะถูกจัดสรรตามรูปแบบการจัดสรรโต๊ะในร้านอาหารจีนค่าที่คาดหวังของตัวแปรสุ่มนี้ โดยกำหนดให้มีจำนวนโต๊ะเท่ากับจำนวนโต๊ะทั้งหมดลูกค้าที่นั่งคือ[ 9 ]
ที่ไหนคือฟังก์ชันไดแกมมาสำหรับกรณีสองพารามิเตอร์ สำหรับจำนวนโต๊ะที่ถูกจองโดยประมาณคือ[ 7 ]
ที่ไหนคือแฟกทอเรียลที่เพิ่มขึ้น (ตามที่นิยามไว้ข้างต้น)
แบบจำลอง Dirichlet-categorical
สำหรับการเลือกพารามิเตอร์และ, ที่ไหนกระบวนการร้านอาหารจีนแบบสองพารามิเตอร์นั้นเทียบเท่ากับแบบจำลอง Dirichlet-categoricalซึ่งเป็นแบบจำลองลำดับชั้นที่สามารถกำหนดได้ดังนี้ โปรดสังเกตว่าสำหรับการตั้งค่าพารามิเตอร์นี้ ความน่าจะเป็นของการเข้าใช้โต๊ะใหม่เมื่อมีโต๊ะว่างอยู่แล้วจำนวนโต๊ะที่มีคนนั่งอยู่เป็นศูนย์ ดังนั้นจำนวนโต๊ะที่มีคนนั่งอยู่จึงมีค่าสูงสุดไม่เกินหากเราเลือกที่จะระบุตารางด้วยป้ายกำกับที่รับค่าในจากนั้นจึงสร้างการแบ่งแบบสุ่มของชุดข้อมูลแบบจำลองลำดับชั้นจะสร้างการกระจายป้ายกำกับตามหมวดหมู่ ขึ้นมาก่อนจากการกระจายแบบ Dirichlet สมมาตร โดยมีพารามิเตอร์ความเข้มข้นจากนั้น ดำเนินการแยกกันสำหรับแต่ละส่วนลูกค้า ป้ายกำกับตารางถูกดึงมาจากหมวดหมู่เนื่องจากการแจกแจงแบบ Dirichlet เป็นการแจกแจงแบบคู่ควบกับการแจกแจงแบบจัดกลุ่ม ดังนั้นตัวแปรที่ซ่อนอยู่สามารถแยกส่วนย่อยออกเพื่อให้ได้การแจกแจงการทำนายภายหลังสำหรับสถานะป้ายกำกับถัดไป ได้, ที่ให้ไว้ป้ายกำกับก่อนหน้า
ที่ไหนคือจำนวนลูกค้าที่นั่งอยู่ที่โต๊ะแล้ว. กับและซึ่งสอดคล้องกับสูตรทั่วไปข้างต้นสำหรับความน่าจะเป็นของการนั่งที่โต๊ะที่มีคนนั่งอยู่ เมื่อความน่าจะเป็นที่จะนั่งที่ใดที่หนึ่งโต๊ะที่ว่างอยู่ก็สอดคล้องกับสูตรทั่วไปและกำหนดโดย
ความน่าจะเป็นแบบมาร์จินัลสำหรับป้ายกำกับต่างๆ กำหนดโดย
ที่ไหนและคือแฟกทอเรียลที่เพิ่มขึ้นโดยทั่วไปแล้วจะมีสถานะป้ายกำกับหลายสถานะที่สอดคล้องกับ พาร์ติชัน เดียวกันสำหรับพาร์ติชันที่กำหนดซึ่งมีบล็อก จำนวนสถานะป้ายกำกับที่สอดคล้องกับพาร์ติชันนี้ทั้งหมด กำหนดโดย แฟกทอเรี ยลแบบลดลงเมื่อพิจารณาถึงเรื่องนี้แล้ว ความน่าจะเป็นสำหรับการแบ่งกลุ่มคือ
ซึ่งสามารถตรวจสอบได้ว่าสอดคล้องกับความน่าจะเป็นของการแบ่งส่วนทั่วไปที่ระบุไว้ข้างต้นในรูปของสัญลักษณ์ k ของ Pochhammer โปรดสังเกตอีกครั้งว่า ถ้าอยู่นอกเหนือการสนับสนุน เช่นแฟกทอเรียลที่ลดลงมีค่าเป็นศูนย์ตามที่ควรจะเป็น (การใช้งานจริงที่ประเมินความน่าจะเป็นของลอการิทึมสำหรับพาร์ติชันผ่าน จะกลับมาเมื่อใดก็ตามที่(ตามความจำเป็น)
ความสัมพันธ์ระหว่าง CRP แบบ Dirichlet-categorical และ CRP แบบพารามิเตอร์เดียว
ในอีกด้านหนึ่ง ลองพิจารณาขั้นตอนการดำเนินงานร้านอาหารจีนที่มีพารามิเตอร์เดียวและซึ่งเราใช้สัญลักษณ์แทนและในอีกด้านหนึ่ง แบบจำลอง Dirichlet-categorical ที่มีจำนวนเต็มบวกและตำแหน่งที่เราเลือกซึ่งดังที่แสดงไว้ข้างต้น เทียบเท่ากับสิ่งนี้แสดงให้เห็นว่าแบบจำลอง Dirichlet-categorical สามารถทำให้ใกล้เคียงกับค่าที่กำหนดได้มากเท่าใดก็ได้โดยการสร้างใหญ่.
กระบวนการหักแท่ง
กระบวนการร้านอาหารจีนสองพารามิเตอร์สามารถกำหนดได้อย่างเทียบเท่าในแง่ของกระบวนการหักแท่ง [ 10 ] สำหรับกรณีที่และกระบวนการหักแท่งสามารถอธิบายได้ว่าเป็นแบบจำลองลำดับชั้น คล้ายกับแบบจำลอง Dirichlet-categorical ข้างต้น ยกเว้นว่ามีสถานะป้ายกำกับจำนวนอนันต์ ป้ายกำกับตารางถูกสุ่มอย่างอิสระจากการกระจายเชิงหมวดหมู่ที่ไม่มีที่สิ้นสุดโดยส่วนประกอบต่างๆ จะถูกสุ่มตัวอย่างโดยใช้การหักแท่งไม้ : เริ่มต้นด้วยแท่งไม้ที่มีความยาว 1 และหักมันออกเป็นสองส่วนโดยสุ่ม ความยาวของครึ่งซ้ายคือและครึ่งขวาถูกแบ่งออกอีกครั้งแบบวนซ้ำเพื่อให้ได้กล่าวให้แม่นยำยิ่งขึ้นคือ เศษส่วนทางซ้ายของ-th break ถูกสุ่มมาจากการแจกแจงแบบเบต้า :
ความน่าจะเป็นเชิงหมวดหมู่มีดังนี้:
สำหรับการตั้งค่าพารามิเตอร์และ, ที่ไหนเป็นจำนวนเต็มบวก และในกรณีที่หมวดหมู่มีค่าจำกัด:เราสามารถสุ่มตัวอย่างได้จากการแจกแจงแบบ Dirchlet ทั่วไปตามที่อธิบายไว้ข้างต้นแต่ยังสามารถสุ่มตัวอย่างได้โดยใช้ สูตรการหักแท่ง แบบตัดทอนโดยที่สูตรสำหรับการสุ่มเศษส่วนจะถูกแก้ไขเป็น:
และ.
กระบวนการบุฟเฟต์แบบอินเดีย
เป็นไปได้ที่จะปรับโมเดลเพื่อให้แต่ละจุดข้อมูลไม่เกี่ยวข้องกับคลาสใดคลาสหนึ่งโดยเฉพาะอีกต่อไป (กล่าวคือ เราไม่ได้สร้างพาร์ติชันอีกต่อไป) แต่สามารถเกี่ยวข้องกับคลาสใดๆ ก็ได้ ซึ่งทำให้การเปรียบเทียบโต๊ะในร้านอาหารไม่ชัดเจน จึงเปรียบเทียบกับกระบวนการที่ผู้รับประทานอาหารหลายคนเลือกตัวอย่างจากชุดย่อยของอาหารที่มีให้เลือกมากมายไม่จำกัดในบุฟเฟต์ ความน่าจะเป็นที่ผู้รับประทานอาหารคนใดคนหนึ่งจะเลือกอาหารจานใดจานหนึ่งนั้นเป็นสัดส่วนกับความนิยมของอาหารจานนั้นในหมู่ผู้รับประทานอาหารที่ผ่านมา และนอกจากนี้ผู้รับประทานอาหารอาจเลือกจากอาหารที่ยังไม่ได้ลองด้วย กระบวนการนี้ได้รับการตั้งชื่อว่ากระบวนการบุฟเฟต์แบบอินเดีย และสามารถใช้เพื่ออนุมานคุณลักษณะแฝงในข้อมูลได้[ 11 ]
แอปพลิเคชัน
กระบวนการร้านอาหารจีนมีความเชื่อมโยงอย่างใกล้ชิดกับกระบวนการ Dirichletและแผนการโถของ Pólyaดังนั้นจึงมีประโยชน์ในการประยุกต์ใช้สถิติแบบเบย์เซียนรวมถึงวิธีการแบบเบย์เซียนที่ไม่ใช่พาราเมตริก กระบวนการร้านอาหารจีนแบบทั่วไปมีความเกี่ยวข้องอย่างใกล้ชิดกับกระบวนการ Pitman–Yorกระบวนการเหล่านี้ถูกนำไปใช้ในแอปพลิเคชันมากมาย รวมถึงการสร้างแบบจำลองข้อความ การจัดกลุ่มข้อมูลไมโครอาร์เรย์ ทางชีวภาพ [ 12 ]การสร้างแบบจำลองความหลากหลายทางชีวภาพและการสร้างภาพขึ้นใหม่[ 13 ] [ 14 ]
ดูเพิ่มเติม
ลิงก์ภายนอก
- หนังสือแนะนำเกี่ยวกับการกระจายแบบ Dirichlet และกระบวนการที่เกี่ยวข้อง โดย Frigyik, Kapila และ Gupta
- การบรรยายโดยไมเคิล ไอ. จอร์แดนเกี่ยวกับโครงการ CRP:
- http://videolectures.net/icml05_jordan_dpcrp/