পয়সোঁ বিন্যাস

পরিসংখ্যান ও সম্ভাবনা তত্ত্বে পয়সোঁ বিন্যাস একটি বিচ্ছিন্ন সম্ভাবনা বিন্যাস। ফরাসি গণিতবিদ সিমেওঁ দ্যনি পোয়াসোঁ এর নাম থেকে বিন্যাসটির নাম নেওয়া হয়েছে। বিন্যাসটি নির্দিষ্ট পরিমাণ সময় বা স্থানের ব্যাপ্তিতে ঘটা ঘটনার সংখ্যার সম্ভাবনা প্রকাশ করে, যেখানে ঘটনাগুলো একটি জানা নির্দিষ্ট হারে ঘটে এবং সর্বশেষ ঘটনার পরের সময়ের ওপর অনির্ভরশীল হয় দূরত্ব, ক্ষেত্রফল বা আয়তন বা এ ধরনের অন্য নির্দিষ্ট ব্যাপ্তির ক্ষেত্রেও পয়সোঁ বিন্যাস ব্যবহার করা যেতে পারে।

উদাহরণস্বরূপ, কেউ যদি প্রতিদিন পাওয়া চিঠির পরিমাণের হিসাব রাখেন, তাহলে হয়ত দেখা যাবে প্রতি দিন গড়ে ৪টি চিঠি আসছে। যদি নির্দিষ্ট কোনো চিঠি ভবিষ্যতের কোনো চিঠি আসার সময়কে প্রভাবিত না করে, অর্থাৎ যদি চিঠিগুলো অনেকগুলো আলাদা আলাদা উৎস থেকে স্বাধীনভাবে আসে, তাহলে প্রতি দিন পাওয়া চিঠির সংখ্যা পয়সোঁ বিন্যাস মেনে চলবে ধরে নেওয়া একটি যুক্তিসঙ্গত অনুমান হবে। পয়সোঁ বিন্যাসের অন্যান্য উদাহরণের মধ্যে রয়েছে কোনো কল সেন্টারে প্রতি ঘণ্টায় আসা ফোন কলের সংখ্যা ও কোনো তেজস্ক্রিয় উৎস থেকে প্রতি সেকেন্ডে ক্ষয়কৃত কণার সংখ্যা।

মৌলিক ধারণা

নির্দিষ্ট পরিমাণ সময় বা স্থানের ব্যাপ্তিতে ঘটা ঘটনার সংখ্যার মডেল তৈরিতে পয়সোঁ বিন্যাস খুব জনপ্রিয়।

উদাহরণ

নিচের ঘটনার ক্ষেত্রে পয়সোঁ বিন্যাসের সাহায্যে মডেল তৈরি করা যেতে পারে-

• প্রতি বছর এক মিটারের বেশি ব্যাসের যে পরিমাণ উল্কাপিণ্ড পৃথিবীর বুকে আঘাত হানে
• সকাল ১০টা থেকে রাত ১১টার মধ্যে যে পরিমাণ রোগী ইমারজেন্সি কক্ষে আসেন
• নির্দিষ্ট পরিমাণ সময়ের মধ্যে যে পরিমাণ ফোটন একটি ডিটেক্টরে ধরা পড়ে

অনুমান ও বৈধতা

নিচের অনুমানগুলো সত্য হলে পয়সোঁ বিন্যাসকে উপযুক্ত মডেল হিসেবে বিবেচনা করা যাবে।

• কোনো একটি নির্দিষ্ট ব্যাপ্তিতে ঘটনার সংখ্যাকে টেমপ্লেট:Mvar দিয়ে প্রকাশ করলে যদি টেমপ্লেট:Mvar এর মান হতে পারে ০, ১, ২, ...।
• কোনো ঘটনা পরবর্তী ঘটনা ঘটার সম্ভাবনাকে প্রভাবিত করবে না। অর্থাৎ, ঘটনাগুলো হবে স্বাধীন।
• ঘটনা ঘটার গড় হার হবে ধ্রুব।
• দুটি ঘটনা ঠিক একই সময়ে ঘটবে না। বরং, অতি ক্ষুদ্র ব্যাপ্তিতে কোনো ঘটনা হয় ঘটবে নয়ত ঘটবে না।

অথবা

• প্রকৃত সম্ভাবনা বিন্যাস হবে দ্বিপদী বিন্যাস এবং চেষ্টার (trial) সংখ্যা সংশ্লিষ্ট সফলতার সংখ্যার চেয়ে যথেষ্ট বড় হবে।

এই শর্তগুলো সত্য হলে টেমপ্লেট:Mvar হবে একটি পয়সোঁ দৈব চলক আর টেমপ্লেট:Mvar এর বিন্যাস হবে একটি পয়সোঁ বিন্যাস।

পয়সোঁ বিন্যাসে ঘটনার সম্ভাবনা

কোনো ব্যাপ্তিতে একটি ঘটনা ০, ১, ২, ... বার ঘটতে পারে। ব্যাপ্তির ঘটনার গড় সংখ্যাকে $\lambda$ (ল্যামডা) দ্বারা প্রকাশ করা হয়। $\lambda$ হলো ঘটনার হার, যাকে হার পরামিতিও বলা হয়। কোনো ব্যাপ্তিতে টেমপ্লেট:Mvar সংখ্যক ঘটনা পর্যবেক্ষণ করার সম্ভাবনা

$P(k\text{ events in interval})=e^{-\lambda }\frac{{\lambda }^k}{k!}$

সমীকরণটি দিয়ে প্রকাশ করা যাবে, যেখানে

• $\lambda$ হলো প্রতি ব্যাপ্তিতে ঘটনার গড় সংখ্যা
• e হলো একটি সংখ্যা, যার মান 2.71828... (অয়লার সংখ্যা ও প্রাকৃতিক লগারিদমের ভিত্তি)
• টেমপ্লেট:Mvar এর মান হতে পারে ০, ১, ২, ...
• টেমপ্লেট:Mvar! = টেমপ্লেট:Mvar × (টেমপ্লেট:Mvar − 1) × (টেমপ্লেট:Mvar − 2) × … × 2 × 1 হলো টেমপ্লেট:Mvar এর ফ্যাক্টরিয়াল।

এই সমীকরণটি হলো পয়সোঁ বিন্যাসের সম্ভাবনা ভর ফাংশন। লক্ষ্যনীয় যে, গড় ঘটনা $\lambda$ এর বদলে ঘটনা ঘটার সময়ের হার $r$ দেওয়া থাকলেও সমীকরণটির পরিবর্তিত রূপ ব্যবহার করা যাবে। সেক্ষেত্রে $\lambda =rt$ হবে (যেখানে $r$ এর একক হলো ১/সময়) এবং

$P(k\text{ events in interval }t)=e^{-rt}\frac{(rt{)}^k}{k!}$

পয়সোঁ বিন্যাসে সম্ভাবনার উদাহরণ

কোনো একটি নির্দিষ্ট নদীতে গড়ে প্রতি একশত বছরে একবার অতিপ্রবাহের কারণে বন্যা হয়। পয়সোঁ মডেলকে উপযুক্ত ধরে নিয়ে ১০০ বছরের ব্যাপ্তিতে এমন টেমপ্লেট:Mvar = ০, ১, ২, ৩, ৪, ৫ বা ৬টি বন্যা হবে তার সম্ভাবনা বের করা সম্ভব। এখানে গড় ঘটনার হার হলো প্রতি ১০০ বছরে একটি অতিপ্রবাহ। অর্থাৎ, λ = 1

$P(k\text{ overflow floods in 100 years})=\frac{{\lambda }^k{e}^{-\lambda }}{k!}=\frac{1^k{e}^{-1}}{k!}$

$P(k=0\text{ overflow floods in 100 years})=\frac{1^0{e}^{-1}}{0!}=\frac{e^{-1}}{1}\approx 0.368$

$P(k=1\text{ overflow flood in 100 years})=\frac{1^1{e}^{-1}}{1!}=\frac{e^{-1}}{1}\approx 0.368$

$P(k=2\text{ overflow floods in 100 years})=\frac{1^2{e}^{-1}}{2!}=\frac{e^{-1}}{2}\approx 0.184$

১০০ বছর সময়কালের মধ্যে ০ থেকে ৬টি অতিপ্রবাহের সম্ভাবনা নিচের সারণিতে দেওয়া আছে।

টেমপ্লেট:Mvar	P(টেমপ্লেট:Mvar overflow floods in 100 years)
0	0.368
1	0.368
2	0.184
3	0.061
4	0.015
5	0.003
6	0.0005

উগারতে ও তার সহকর্মীরা জানিয়েছেন, ফুটবল বিশ্বকাপের একটি ম্যাচে গড় গোলের সংখ্যা প্রায় ২.৫ এবং পয়সোঁ মডেলের ব্যবহার যথাযথ। যেহেতু প্রতি ম্যাচে গোলের গড় সংখ্যা ২.৫, অতএব λ = 2.5।

$P(k\text{ goals in a match})=\frac{2.5^k{e}^{-2.5}}{k!}$

$P(k=0\text{ goals in a match})=\frac{2.5^0{e}^{-2.5}}{0!}=\frac{e^{-2.5}}{1}\approx 0.082$

$P(k=1\text{ goal in a match})=\frac{2.5^1{e}^{-2.5}}{1!}=\frac{2.5e^{-2.5}}{1}\approx 0.205$

$P(k=2\text{ goals in a match})=\frac{2.5^2{e}^{-2.5}}{2!}=\frac{6.25e^{-2.5}}{2}\approx 0.257$

নিচের সারণিতে কোনো ম্যাচে ০ থেকে ৭টি গোল হবার সম্ভাবনা দেওয়া আছে।

টেমপ্লেট:Mvar	P(টেমপ্লেট:Mvar goals in a World Cup soccer match)
0	0.082
1	0.205
2	0.257
3	0.213
4	0.133
5	0.067
6	0.028
7	0.010

প্রতি ব্যাপ্তিতে একবার ঘটা ঘটনা: λ = 1 ও k = 0 এর বিশেষ অবস্থা

ধরা যাক, জ্যোতির্বিদগণ হিসেব করে পেলেন যে বড় বড় উল্কাপিণ্ড (নির্দিষ্ট আকারের চেয়ে বড়) প্রতি ১০০ বছরে একবার পৃথিবীতে আঘাত হানে (প্রতি ১০০ বছরে λ = 1টি ঘটনা)। আরও দেখলেন যে পৃথিবীকে আঘাত করা এই আকারে উল্কাপিণ্ডের সংখ্যা পয়সোঁ বিন্যাস মেনে চলে। তাহলে পরবর্তী ১০০ বছরে টেমপ্লেট:Mvar = 0টি উল্কাপিন্ড আঘাত হানবে তার সম্ভবনা কত?

$P(k=\text{0 meteorites hit in next 100 years})=\frac{1^0{e}^{-1}}{0!}=\frac{1}{e}\approx 0.37$

এ অনুমানগুলো মেনে নিলে দেখা যায়, পরবর্তী ১০০ বছরে বড় কোনো উল্কাপিণ্ড পৃথিবীতে আঘাত হানবে না এমন সম্ভাবনা প্রায় ০.৩৭। বাকি ১ − ০.৩৭ = ০.৬৩ হলো পরবর্তী ১০০ বছরে ১, ২, ৩ বা আরও বেশি সংখ্যক বড় উল্কাপিণ্ড আঘাত হানার সম্ভাবনা। ওপরের একটি উদাহরণে অতিপ্রবাহজনিত বন্যা প্রতি ১০০ বছরে ১ বার ঘটেছিল (λ = 1)। একই হিসাব অনুসারেই ১০০ বছরে অতিপ্রবাহজনিত কোনো বন্যা না হবার সম্ভাবনা ছিল ০.৩৭। সাধারণভাবে প্রতি ব্যাপ্তিতে কোনো ঘটনা গড়ে একবার ঘটলে (λ = 1) এবং ঘটনাগুলো পয়সোঁ বিন্যাস মেনে চললে টেমপ্লেট:Nowrap। এছাড়া, P(পরবর্তী ব্যাপ্তিতে শুধু একটি ঘটনা) = ০.৩৭, যেটা অতিপ্রবাহজনিত বন্যার সারণিতে দেখানো হয়েছে।

পয়সোঁ অনুমান মেনে চলে না এমন উদাহরণ

কোনো ছাত্র পরিষদে প্রতি মিনিটে উপস্থিত হওয়া ছাত্রদের সংখ্যা পয়সোঁ বিন্যাস নাও মেনে চলতে পারে, কারণ এখানে হার ধ্রুবক নয় (ক্লাস চলাকালে হার কম এবং ক্লাসের ফাঁকে হার বেশি)। আবার ছাত্রদের আসার ঘটনা স্বাধীনও নয় (ছাত্ররা সাধারণত দল বেঁধে আসে)। একটি বড় ভূমিকম্পের কারণে সমমাত্রার আফটারশকের সম্ভাবনা বেড়ে গেলে কোনো দেশে প্রতি বছর সংঘটিত ৫ মাত্রার ভূমিকম্পের সংখ্যা পয়সোঁ বিন্যাস নাও মেনে চলতে পারে। কোনো হাসপাতালের নিবিড় পরিচর্যা কেন্দ্রে ভর্তি রোগীদের ক্ষেত্রে অবস্থানের দিনের সংখ্যা পয়সোঁ বিন্যাস মেনে চলবে না, কারণ দিনের সংখ্যা শূন্য হওয়া সম্ভব নয়। এই বিন্যাসকে শূন্য-বিহীন পয়সোঁ বিন্যাসের সাহায্যে মডেল করা যেতে পারে। যে সকল গণনা বিন্যাসে শূন্যটি ঘটনার ব্যাপ্তির সংখ্যা পয়সোঁ মডেলের অনুমানের চেয়ে বেশি সেক্ষেত্রে শূন্য-স্ফীত মডেল ব্যবহার করা যেতে পারে।

পয়সোঁ নির্ভরণ ও ঋণাত্মক দ্বিপদী নির্ভরণ

অধীন চলক গণনাবাচক হলে অর্থাৎ কোনো ব্যাপ্তিতে ঘটনার সংখ্যা (০, ১, ২, ...) হলে পয়সোঁ নির্ভরণ ও ঋণাত্মক দ্বিপদী নির্ভরণ খুব ভালোভাবে কাজে লাগানো যায়।

ইতিহাস

সিমেওঁ দ্যনি পোয়াসোঁ (১৭৮১-১৮৪০) তাঁর ১৮৩১ সালের Recherches sur la probabilité des jugements en matière criminelle et en matière civile (ফৌজদারী ও বেসামরিক বিষয়াদি সম্পর্কিত রায়ের সম্ভাবনা বিষয়ক গবেষণা) কাজে সম্ভাবনা তত্ত্বের সাথে বিন্যাসটি সর্বপ্রথম প্রবর্তন ও প্রকাশ করেন। কাজটিতে একটি দৈব চলক টেমপ্লেট:Mvar ব্যবহারের মাধ্যমে কোনো নির্দিষ্ট দেশের অন্যায় রায়ের সংখ্যা সম্পর্কে তত্ত্ব দেওয়া হয়। এই টেমপ্লেট:Mvar অন্যান্য জিনিসের মধ্যে গণনা করে যে একটি নির্দিষ্ট সময় ব্যাপ্তিতে কতটি বিচ্ছিন্ন ঘটনা ঘটেছে। এই ফলাফল এর আগে আব্রাআম দ্য মোয়াভ্র্‌ও (১৭১১) ফিলোসোফিকেল ট্রাঞ্জেকশন অব রয়েল সোসাইটি জার্নালে De Mensura Sortis seu; de Probabilitate Eventuum in Ludis a Casu Fortuito Pendentibus নামে প্রকাশিত নিবন্ধে দেখিয়েছিলেন। এ কারণে এটি স্টিগলারের নিয়মেরও একটি উদাহরণ। আর এ কারণে অনেকে বলেছেন পয়সোঁ বিন্যাসে দ্য মোয়াভ্ররের নাম থাকা উচিত। ১৮৯৮ সালে বিন্যাসটির একটি বাস্তব প্রয়োগ দেখিয়েছেন লাদিসলাউস বরতিকিউইসজ। তাঁকে প্রুশিয়ান সেনাবাহিনীতে ঘোড়ার পদাঘাতে দূর্ঘটনাক্রমে মারা যাওয়া সৈন্যের সংখ্যা নিয়ে তদন্ত করতে বলা হলে সে কাজের অংশ হিসেবে তিনি পয়সোঁ বিন্যাসের প্রয়োগ ঘটনা। এই পরীক্ষণের মাধ্যমে পয়সোঁ বিন্যাস নির্ভরযোগ্যতা প্রকৌশল শাস্ত্রে অন্তর্ভূক্ত হয়।

সংজ্ঞা

একটি বিচ্ছিন্ন দৈব চলক Xটেমপ্লেট:Space, λ > 0 পরামিতির একটি পয়সোঁ বিন্যাস মেনে চলবে যদি k = 0, 1, 2, ..., এর জন্য টেমপ্লেট:Mvar এর সম্ভাবনা ভর ফাংশন এ রকম হয়:

$f(k;\lambda )=\mathrm{Pr}(X=k)=\frac{{\lambda }^k{e}^{-\lambda }}{k!},$

যেখানে

• e হলো অয়লার সংখ্যা (e = 2.71828...)।
• k! হলো k এর ফ্যাক্টোরিয়াল।

ধনাত্মক বাস্তব সংখ্যা λ হবে X এর প্রত্যাশিত মান ও ভেদাঙ্কের সমান।

$\lambda =\mathrm{E}(X)=\mathrm{Var}(X).$

অনেক বেশি সংখ্যক সম্ভাব্য ঘটনার প্রতিটি দুর্লভ ঘটনা হলে পয়সোঁ বিন্যাস ব্যবহার করা যেতে পারে। একটি নির্দিষ্ট সময় ব্যাপ্তিতে ঐ রকম কতগুলো ঘটনা ঘটবে সেটা উপযুক্ত পরিস্থিতিতে পয়সোঁ বিন্যাসের একটি দৈব সংখ্যা হবে। পয়সোঁ বিন্যাসের প্রচলিত সংজ্ঞায় এমন দুটি এমন দুটি পদ আছে যেগুলোর কারণে কম্পিউটার দিয়ে বিন্যাসটির কাজ করা সহজেই অসম্ভব হয়ে পড়ে। এগুলো হলো λ^k ও k!। এছাড়াও λ^k কে k! দ্বারা ভাগ দিলে যে আসন্নীকরণ ত্রুটি ঘটে তাও e^−λ এর তুলনায় অনেক বড়। এ কারণে ভুল ফলাফল পাওয়া যায়। সংখ্যাভিত্তিক স্থিতিশীলতার জন্য পয়সোঁ সম্ভাবনা ভর ফাংশন এভাবে বের করা উচিত:

$f(k;\lambda )=\exp \left\{k\ln \lambda -\lambda -\ln \mathrm{\Gamma }(k+1)\right\},$

যা গাণিতিকভাবে সমতুল্য কিন্তু সংখ্যাভিত্তিকভাবে স্থিতিশীল। C প্রোগ্রামিং ল্যাংগুয়েজের আদর্শ লাইব্রেরি (C99) থেকে lgamma ফাংশন, R প্রোগ্রামিং ল্যাংগুয়েজ, ম্যাটল্যাবের gammaln ফাংশন, SciPy বা ফোরট্রানের ২০০৮ বা তার পরবর্তী সংস্করণের log_gamma ফাংশন ব্যবহার করে গামা ফাংশনের স্বাভাবিক লগারিদম বের করা যায়।

বৈশিষ্ট্য

বর্ণনামূলক পরিসংখ্যান

• পয়সোঁ বিন্যাস মেনে চলা দৈব চলকের প্রত্যাশিত মান ও ভেদাঙ্ক দুটিই λ।
• বিভেদাঙ্ক হলো ${\lambda }^{-1/2}$ এবং বিস্তার সূচকের মান ১।
• গড় পরম ব্যবধান হলো:

$\mathrm{E}|X-\lambda |=2\exp (-\lambda )\frac{{\lambda }^{\lfloor \lambda \rfloor +1}}{\lfloor \lambda \rfloor !}.$

• λ অপূর্ণ সংখ্যা হলে পয়সোঁ বিন্যাসের দৈব চলকের প্রচুরক হবে $\lfloor \lambda \rfloor$, যার অর্থ হলো λ এর সমান বা ছোট পূর্ণ সংখ্যা। একে floor(λ) আকারেও লেখা হয়। λ ধনাত্মক পূর্ণসংখ্যা হলে প্রচুরক হবে λ ও λ-1।
• পয়সোঁ বিন্যাসের সবগুলো ক্রমযোজিত মান প্রত্যাশিত মান λ এর সমান হবে। পয়সোঁ বিন্যাসের nতম ফ্যাক্টোরিয়াল পরিঘাত হলো λⁿ।
• পয়সোঁ প্রক্রিয়ার প্রত্যাশিত মানকে অনেক সময় তীব্রতা ও এক্সপোজারের (যাকে আরও সাধারণভাবে বলা যায় সময় বা স্থানের সাপেক্ষে তীব্রতা ফাংশনের ইন্টিগ্রাল) গুণফল আকারে লেখা হয়।

মধ্যমা

বিন্যাসটির মধ্যমা (ν) জানা আছে এবং এটি খুব তীক্ষ্ণ:

$\lambda -\ln 2\le \nu <\lambda +\frac{1}{3}.$

উচ্চতর পরিঘাত

• মূলের সাপেক্ষে পয়সোঁ বিন্যাসের উচ্চতর পরিঘাত m_kগুলো হলো λ এর তাচার্দ বহুপদী:

$m_k=\sum _{i=0}^k{\lambda }^i\left\{\begin{array}{c}k\\ i\end{array}\right\},$

যেখানে দ্বিতীয় বন্ধনীর ভেতরের পদটি হলো দ্বিতীয় প্রকারের স্টারলিং সংখ্যা। বহুপদীর সহগগুলো আসে সমাবেশের সূত্র থেকে। বস্তুত, পয়সোঁ বিন্যাসের প্রত্যাশিত মান 1 হলে দোবিন্সকির সূত্র অনুসারে nতম পরিঘাত হবে n আকারের একটি সেটের বিভাজনসংখ্যার সমান।

পয়সোঁ দৈব চলকের সমষ্টি

যদি $i=1,\dots ,n$ এর জন্য $X_i\sim \mathrm{Pois}({\lambda }_i)$ স্বাধীন হয় এবং $\lambda =\sum _{i=1}^n{\lambda }_i$ হয়, তাহলে $Y=\left(\sum _{i=1}^n{X}_i\right)\sim \mathrm{Pois}(\lambda )$। বিপরীতভাবে বিবেচনা করলে পাওয়া যাবে রাইকভের উপপাদ্য, যেটি অনুসারে দুটি স্বাধীন দৈব চলকের সমষ্টি পয়সোঁ বিন্যাস মেনে চললে স্বাধীন দৈব চলকদুটিও পয়সোঁ বিন্যাস মেনে চলবে।

অন্যান্য বৈশিষ্ট্য

• পয়সোঁ বিন্যাসগুলো হলো অসীমতক বিভাজ্য সম্ভাবনা বিন্যাস।টেমপ্লেট:R
• Pois(λ) থেকে Pois(λ₀) এর নির্দেশিত কুলব্যাক-লেইবলার অপসরণ হলো

$D_{\text{KL}}(\lambda \mid {\lambda }_0)={\lambda }_0-\lambda +\lambda \log \frac{\lambda }{{\lambda }_0}$।

• পয়সোঁ দৈব চলক $X\sim \mathrm{Pois}(\lambda )$ এর প্রান্তীয় সম্ভাবনার সীমা চেরনোফ সীমা থেকে নির্ণয় করা যায়।

$P(X\ge x)\le \frac{e^{-\lambda }(e\lambda {)}^x}{x^x},\text{ for }x>\lambda$,

$P(X\le x)\le \frac{e^{-\lambda }(e\lambda {)}^x}{x^x},\text{ for }x<\lambda .$

• যে অসমতাগুলো পয়সোঁ বিন্যাসের দৈব চলক $X\sim \mathrm{Pois}(\lambda )$ এর সম্ভাবনা ফাংশনকে আদর্শ পরিমিত বিন্যাস ফাংশনের সাথে সম্পর্কিত করে সেগুলো হলো:

$\mathrm{\Phi }(sign(k-\lambda )\sqrt{2D_{\text{KL}}(k\mid \lambda )})<P(X\le k)<\mathrm{\Phi }(sign(k-\lambda +1)\sqrt{2D_{\text{KL}}(k+1\mid \lambda )}),\text{ for }k>0.$

যেখানে $D_{\text{KL}}(k\mid \lambda )$ হলো ওপরে উল্লিখিত নির্দেশিত কুলব্যাক-লেইবলার।

পয়সোঁ রেস

ধরা যাক, $X\sim \mathrm{Pois}(\lambda )$ ও $Y\sim \mathrm{Pois}(\mu )$ স্বাধীন দৈব চলক, যেখানে $\lambda <\mu$, তাহলে

$\frac{e^{-(\sqrt{\mu }-\sqrt{\lambda }{)}^2}}{(\lambda +\mu {)}^2}-\frac{e^{-(\lambda +\mu )}}{2\sqrt{\lambda \mu }}-\frac{e^{-(\lambda +\mu )}}{4\lambda \mu }\le P(X-Y\ge 0)\le e^{-(\sqrt{\mu }-\sqrt{\lambda }{)}^2}$

এখানে ঊর্ধ্ব-সীমা প্রমাণ করা হয় আদর্শ চেরনোফ সীমা দিয়ে।

আবার $P(X-Y\ge 0\mid X+Y=i)$ হলো $Z\ge \frac{i}{2}$ এর সম্ভাবনা, যেখানে $Z\sim \mathrm{Bin}(i,\frac{\lambda }{\lambda +\mu })$, যার নিম্ন-সীমা হলো $\frac{1}{(i+1{)}^2}{e}^{(-iD(0.5\Vert \frac{\lambda }{\lambda +\mu }))}$, যেখানে $D$ হলো আপেক্ষিক এনট্রপি। এভাবে নিম্ন-সীমাও প্রমাণ করা যায়। এছাড়াও $X+Y\sim \mathrm{Pois}(\lambda +\mu )$ সম্পর্ক দিয়ে শর্তহীন সম্ভাবনার নিম্ন-সীমা বের করলে ফলাফলটি পাওয়া যায়।

সম্পর্কিত বিন্যাস

• $X_1\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}({\lambda }_1)$ ও $X_2\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}({\lambda }_2)$ স্বাধীন হলে $Y=X_1-X_2$ পার্থক্যটি স্কেলাম বিন্যাস মেনে চলবে।
• $X_1\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}({\lambda }_1)$ ও $X_2\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}({\lambda }_2)$ স্বাধীন হলে $X_1+X_2$ শর্তে $X_1$ একটি দ্বিপদী বিন্যাস মেনে চলবে।

নির্দিষ্ট করে উল্লেখ করলে, $X_1+X_2=k$ হলে $X_1\sim \mathrm{B}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}(k,{\lambda }_1/({\lambda }_1+{\lambda }_2))$

আরও সাধারণভাবে, X₁, X₂,..., X_n λ₁, λ₂,..., λ_n এর স্বাধীন পয়সোঁ দৈব চলক হলে এবং

$\sum _{j=1}^n{X}_j=k$ হলে $X_i\sim \mathrm{B}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}(k,\frac{{\lambda }_i}{\sum _{j=1}^n{\lambda }_j})$। বস্তুত, $\{X_i\}\sim \mathrm{M}\mathrm{u}\mathrm{l}\mathrm{t}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}(k,\left\{\frac{{\lambda }_i}{\sum _{j=1}^n{\lambda }_j}\right\})$।

• যদি $X\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}(\lambda )$ এবং X = k শর্তে $Y$ একটি দ্বিপদী বিন্যাস $Y\mid (X=k)\sim \mathrm{B}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}(k,p)$ মেনে চলে তাহলে Y একটি পয়সোঁ বিন্যাস $Y\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}(\lambda \cdot p)$ মেনে চলবে। বস্তুত, X = k শর্তে $\{Y_i\}$ একটি বহুপদী বিন্যাস মেনে চলবে। অর্থাৎ, $\{Y_i\}\mid (X=k)\sim \mathrm{M}\mathrm{u}\mathrm{l}\mathrm{t}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}(k,p_i)$ হলে প্রতিটি $Y_i$ একটি স্বাধীন পয়সোঁ বিন্যাস $Y_i\sim \mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}(\lambda \cdot p_i),\rho (Y_i,Y_j)=0$ মেনে চলবে।
• দ্বিপদী বিন্যাসের চেষ্টাসংখ্যা অসীম হলে এবং সফলতার সংখ্যার প্রত্যাশিত মান ধ্রুব থাকলে বিন্যাসটির একটি সীমাস্থ অবস্থা হিসেবে পয়সোঁ বিন্যাস নির্ণয় করা যায়। অতএব, দ্বিপদী বিন্যাসের

n যথেষ্ট বড় হলে এবং pp যথেষ্ট ছোট হলে পয়সোঁ বিন্যাস দিয়ে দ্বিপদী বিন্যাসের খুব কাছাকাছি মান বের করা সম্ভব। একটি সাধারণ নীতি হলো, পয়সোঁ বিন্যাস দিয়ে দ্বিপদী বিন্যাসের মান ভালোভাবে বের করা যাবে যদি n অন্তত ২০ এবং pp ০.০৫ এর সমান বা ছোট হয়। আর n ≥ 100 ও np ≤ 10 হলে আসন্ন মান হবে খুবই নিখুঁত।

$F_{\mathrm{B}\mathrm{i}\mathrm{n}\mathrm{o}\mathrm{m}\mathrm{i}\mathrm{a}\mathrm{l}}(k;n,p)\approx F_{\mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}\mathrm{s}\mathrm{o}\mathrm{n}}(k;\lambda =np)$

• পয়সোঁ বিন্যাস বিচ্ছিন্ন যৌগিক পয়সোঁ বিন্যাসের একটিমাত্র পরমাতিযুক্ত একটি বিশেষ অবস্থা। একচলকযুক্ত বহুপদী বিন্যাসের সীমাস্থ বিন্যাস থেকে বিচ্ছিন্ন যৌগিক পয়সোঁ বিন্যাসে পৌঁছা যায়। এটি আবার যৌগিক পয়সোঁ বিন্যাসের একটি বিশেষ অবস্থা।
• λ এর বড় মানের জন্য (ধরা যাক λ>1000) λ গড় ও λ ভেদাঙ্কের (অর্থাৎ, পরিমিতি ব্যবধান $\sqrt{\lambda }$) পরিমিত বিন্যাস দিয়ে পয়সোঁ বিন্যাসের আসন্ন মান খুব ভালোভাবে বের করা যাবে। λ ১০ এর বড় হলে পরিমিত বিন্যাস দিয়ে ভালো আসন্ন মান পাওয়া যাবে যদি অবিচ্ছিন্নতা সংশোধন সঠিকভাবে করা হয়। অর্থাৎ, P(X ≤ x)কে P(X ≤ x + 0.5) দ্বারা প্রতিস্থাপিত করলে, যেখানে X একটি অঋণাত্মক পূর্ণসংখ্যা

$F_{\mathrm{P}\mathrm{o}\mathrm{i}\mathrm{s}\mathrm{s}\mathrm{o}\mathrm{n}}(x;\lambda )\approx F_{\mathrm{n}\mathrm{o}\mathrm{r}\mathrm{m}\mathrm{a}\mathrm{l}}(x;\mu =\lambda ,{\sigma }^2=\lambda )$

• ভেদাঙ্ক স্থিতিশীলকরণ রুপান্তর: কোনো চলক পয়সোঁ বিন্যাস মেনে চললে এর বর্গমূলের বিন্যাস পরিমিত বিন্যাসের কাছাকাছি হবে, যেখানে প্রত্যাশিত মান হবে $\sqrt{\lambda }$ এবং ভেদাঙ্ক হবে ১/৪।টেমপ্লেট:R এই রূপান্তরের ফলে অরূপান্তরিত চলকের চেয়ে দ্রুত হারে পরিমিত বিন্যাসের বৈশিষ্ট্য (λ বাড়ার সাথে সাথে) অর্জন করা যায়। ভেদাঙ্ক স্থিতিশীলকরণ রুপান্তরের আরও কিছু ও কিছুটা জটিল পদ্ধতিও আছে। এর মধ্যে অন্যতম হলো অ্যান্সকম্ব রূপান্তর।
• যদি প্রতিটি t > 0 এর জন্য সময় ব্যাপ্তি [0, t]তে ঘটনার সংখ্যা λt গড়ের পয়সোঁ বিন্যাস মেনে চলে তাহলে ঘটনা ঘটনার মধ্যবর্তী সময়ের অনুক্রম স্বাধীন হবে এবং সবাই 1/λ গড়ের সূচকীয় বিন্যাস মেনে চলবে।
• পয়সোঁ ও কাই-বর্গ বিন্যাসের ক্রমযোজিত বিন্যাস ফাংশনের সম্পর্ক এমন:

$F_{\text{Poisson}}(k;\lambda )=1-F_{{\chi }^2}(2\lambda ;2(k+1))\text{ integer }k,$

ও

$\mathrm{Pr}(X=k)=F_{{\chi }^2}(2\lambda ;2(k+1))-F_{{\chi }^2}(2\lambda ;2k).$

সংঘটন

গণনা বিষয়ক অসংখ্য শাখায় পয়সোঁ বিন্যাসের প্রয়োগ দেখা যায়:

• টেলিযোগাযোগে উদাহরণ: কোনো সিস্টেমে আসা কলের সংখ্যা
• জ্যোতির্বিদ্যায় উদাহরণ: টেলিস্কোপে আসা ফোটন কণার সংখ্যা
• রসায়নে উদাহরণ: জীবন্ত পলিমারকরণের মোলার ভর বিন্যাস
• জীববিদ্যায় উদাহরণ: ডিএনএ-এর প্রতি একক দৈর্ঘ্যের সুতায় পরিব্যক্তির সংখ্যা।
• ব্যবস্থাপনায় উদাহরণ: কোনো কাউন্টার বা কল সেন্টারে আসা গ্রাহকের সংখ্যা।
• অর্থসংস্থান ও বিমা উদাহরণ: একটি নির্দিষ্ট সময়ে ক্ষতি বা ক্ষতিপূরণের সংখ্যা।
• ভূমিকম্পবিদ্যায় উদাহরণ: বড় ভূমিকম্পের ক্ষেত্রে ঝুঁকি বিষয়ক অসীমতক পয়সোঁ মডেল।
• তেজস্ক্রিয়তায় উদাহরণ: নির্দিষ্ট সময়ে একটি তেজস্ক্রিয় নমুনার ক্ষয়ের সংখ্যা।

পয়সোঁ বিন্যাস মূলত পয়সোঁ প্রক্রিয়া থেকে আসে। বিচ্ছিন্ন বৈশিষ্ট্যযুক্ত (অর্থাৎ, যে ঘটনাগুলো কোনো নির্দিষ্ট সময় বা স্থানে ০, ১, ২, ৩, ... ইত্যাদি বার ঘটতে পারে) বিভিন্ন ঘটনার ক্ষেত্রে এটি ব্যবহার করা যায় যদি ঐ ঘটনা ঘটার সম্ভাবনা সময় বা স্থানে ধ্রুব হয়। পয়সোঁ মডেল দিয়ে ব্যাখ্যা করা যাবে এমন কিছু উদাহরণ হলো:

• প্রুশিয়ান ঘোড়সওয়ার বাহিনীর প্রতিটি বিভাগে প্রতি বছর ঘোড়ার পদাঘাতে মৃত সৈন্যের সংখ্যা। লাদিসলাউস বরতিকিউইসজ (১৮৬৮-১৯৩১) একটি বইয়ে এই উদাহরণ ব্যবহার করেছেন।
• গিনেস বিয়ার চোলাইয়ের সময় ব্যবহৃত ঈস্ট কোষের সংখ্যা। এই উদাহরণ ব্যবহার করেছিলেন উইলিয়াম সিলি গসেট (১৮৭৬-১৯৩৭)।
• এক মিনিটের মধ্যে কোনো কল সেন্টারে আসা কলের সংখ্যা। উদাহরণটির বিবরণ দিয়েছিলেন এ. কে. এরল্যাং।
• ইন্টারনেট ট্র্যাফিক।
• দুটি দলের কোনো খেলায় গোলের সংখ্যা।

• কোনো নির্দিষ্ট বয়সসীমার মানুষের মধ্যে প্রতি বছর মৃত্যুর সংখ্যা।
• কোনো নির্দিষ্ট সময়ের মধ্যে শেয়ার বাজারের ঊর্ধ্বগতির সংখ্যা।
• সমধর্মীতা অনুমান সঠিক ধরে নিলে প্রতি মিনিটে একটি ওয়েব সার্ভারে যতবার প্রবেশ করা হয়।
• নির্দিষ্ট পরিমাণ বিকিরণের পরে একটি নির্দিষ্ট পরিমাণ ডিএনএ-তে ঘটা পরিব্যক্তির সংখ্যা।
• সংক্রমণ ঘটানো জীবাণু ও সংক্রমণযোগ্য কোষের হার ধ্রুব থাকলে সংক্রমিত কোষের অনুপাত।
• নির্দিষ্ট পরিমাণ তরলে ব্যাকটেরিয়ার সংখ্যা।
• একটি পিক্সেল বর্তনীতে নির্দিষ্ট পরিমাণ আলো ফেললে নির্দিষ্ট সময়ে পৌঁছা ফোটনের সংখ্যা।
• দ্বিতীয় বিশ্বযুদ্ধের সময় লন্ডনের আকাশে ভি-১-ফ্লায়িং বোম্বিং এর পরিমাণ। ১৯৪৬ সালে এ নিয়ে কাজ করেন আর. ডি. ক্লার্ক।

১৯৭৬ সালে গ্যালাঘার দেখান যে ছোট ব্যাপ্তিতে মৌলিক সংখ্যার পরিমাণ পয়সোঁ বিন্যাস মেনে চলে, যদি হার্ডি ও লিটলউডের একটি অপ্রমাণিত অমুমানের নির্দিষ্ট সংস্করণ সত্য হয়।

দুর্লভ ঘটনার বিধি

পয়সোঁ (কাল রেখা) ও বিভিন্ন দ্বিপদী বিন্যাস: n = 10 (লাল বৃত্ত), n = 20 (নীল বৃত্ত), n = 1000 (সবুজ বৃত্ত)। সকল বিন্যাসের গড় ৫। অনুভূমিক অক্ষে ঘটনার সংখ্যা k দেখানো হয়েছে। n যত বড় হয়, একই গড়ের জন্য পয়সোঁ বিন্যাস দ্বিপদী বিন্যাসকে তত ভাল করে আসন্নীকৃত করতে পারে।

কোনো ঘটনা ঘটার হারের সাথে কোনো একটি ছোট উপব্যাপ্তিতে (স্থান, সময় বা অন্য কিছুর) ঐ ঘটনা ঘটার সম্ভাবনার সাথে সম্পর্ক আছে। পয়সোঁ বিন্যাসের ক্ষেত্রে ধরে নেওয়া হয়, ছোট ছোট এমন অনেকগুলো উপব্যাপ্তি আছে যাতে একটি ঘটনা দুইবার ঘটার সম্ভাবনা নগণ্য। এই অনুমানের সাহায্যে দ্বিপদী বিন্যাস থেকে পয়সোঁ বিন্যাস তৈরি করা যায়। এর জন্য শুধু প্রয়োজন পূর্ণ ব্যাপ্তিতে মোট ঘটনার প্রত্যাশিত মান। ধরা যাক, এই সমষ্টি হলো $\lambda$। সম্পূর্ণ ব্যাপ্তিকে এবার $n$ সংখ্যক সমান আকারের উপব্যাপ্তিতে বিভক্ত করা হলো। এরা হলো $I_1,\dots ,I_n$। এখানে $n$ > $\lambda$ হতে হবে (আমরা ব্যাপ্তির খুব সামান্য অংশ নিয়ে কাজ করছি বলে এই অনুমান অর্থবহই বটে)। এর অর্থ হলো প্রতিটি $i$ এর জন্য কোনো ব্যাপ্তি $I_i$-এ প্রত্যাশিত ঘটনার সংখ্যা $\lambda /n$। এবার আমরা ধরে নেব, সম্পূর্ণ ব্যাপ্তিতে কোনো ঘটনার সংঘটন একটি বার্নুলি চেষ্টা, যেখানে $i$তম চেষ্টা হলো $I_i$ উপব্যাপ্তিতে ঘটনা ঘটছে কি না তা দেখা, যার সম্ভাবনা $\lambda /n$। এমন $n$ চেষ্টায় প্রত্যাশিত ঘটনার সংখ্যা হবে $\lambda$, যা সম্পূর্ণ ব্যাপ্তিতে মোট ঘটনার প্রত্যাশিত মান। অতএব ব্যাপ্তির সবগুলো বিভক্ত অংশের জন্য আমরা ঘটনা ঘটাকে $\text{B}(n,\lambda /n)$ আকারের বার্নুলি প্রক্রিয়া দিয়ে আসন্নীকৃত করেছি। আগেও বলা হয়েছে, আমরা খুব ছোট উপব্যাপ্তি নিয়ে কাজ করছি। অতএব আমাদের $n$ এর সীমা হবে অসীমের দিকে। এক্ষেত্রে পয়সোঁ সীমা উপপাদ্যের মাধ্যমে দ্বিপদী বিন্যাস পয়সোঁ বিন্যাসের আসন্ন মান প্রদান করবে।

একটি নির্দিষ্ট ডিএনএ এর ক্রমের পরিব্যাপ্তির সংখ্যাসহ উপরের বেশ কিছু উদাহরণে গণনাকৃত ঘটনা প্রকৃতপক্ষে বিচ্ছিন্ন চেষ্টার ফলাফল। ফলে এদেরকে দ্বিপদী বিন্যাস দিয়ে মডেল করলেই বেশি নিখুঁত ফলাফল পাওয়া যাবে। অর্থাৎ,

$X\sim \text{B}(n,p).$

এক্ষেত্রে n খুব বড় আর p খুব ছোট (ফলে np হবে মাঝামাঝি মানের)। তাহলে বিন্যাসটিকে অপেক্ষাকৃত সহজ উপায়ে পয়সোঁ বিন্যাস দিয়ে আসন্নীকৃত করা যাবে।

$X\sim \text{Pois}(np).$

এই আসন্নীকরণকে অনেকসময় দুর্লভ ঘটনার বিধি বলা হয়, কেননা প্রতিটি আলাদা n বার্নুলি ঘটনা এক একটি দুর্লভ ঘটনা। নামটি কিছুটা বিভ্রান্তিকর, কারণ np ছোট হলে পয়সোঁ প্রক্রিয়ার মোট সফলতার সংখ্যাকে দুর্লভ হতে হবে না। যেমন এক ঘণ্টায় একটি ব্যস্ত সুইচবোর্ডে আসা টেলিফোন কলের সংখ্যা পয়সোঁ বিন্যাস মেনে চলবে। এখানে অপারেটরের কাছে মনে হবে খুব ঘন ঘন কল আসছে। তবে একজন সাধারণ মানুষের কাছে একে দুর্লভ ঘটনা মনে হবে, কারণ ঐ নির্দিষ্ট ঘণ্টায় তিনি ঐ সুইচবোর্ড থেকে কল করবেন তার সম্ভাবনা খুব কম।

অনেক সময় বিধি কথাটিকে সম্ভাবনা বিন্যাসের প্রতিশব্দ হিসেবে ব্যবহার করা হয়। আর বিধির অভিসার বলতে বিন্যাসের অভিসার বোঝানো হয়। এ কারণে পয়সোঁ বিন্যাসকে অনেক সময় ছোট সংখ্যার বিধিও বলা হয়। কারণ এটি দুর্লভ ঘটনার সংখ্যার সম্ভাবনা বিন্যাস যেখানে ঘটনা অনেকভাবে ঘটতে পারে। লাদিসলাউস বরতিকিউইসজ ১৮৯৮ সালে ল অব স্মল নাম্বারস (ছোট সংখ্যার বিধি) নামে পয়সোঁ বিন্যাস নিয়ে একটি বই লিখেছেন।

পয়সোঁ বিন্দু পক্রিয়া

কোনো সসীম অঞ্চলে অবস্থিত কোনো পয়সোঁ বিন্দু প্রক্রিয়ার বিন্দু সংখ্যা থেকে পয়সোঁ বিন্যাস পাওয়া যায়। আরও নির্দিষ্ট করে বললে, D যদি কোনো স্থান হয়, যেমন ইউক্লিডীয় স্থান R^d, যাতে ক্ষেত্রফল, আয়তন বা আরও সার্বিক দৃষ্টিকোণ থেকে বললে লেবেসগ পরিমাপ |D| সসীম অঞ্চল হয় এবং টেমপ্লেট:Nowrap দ্বারা D-তে বিন্দুর সংখ্যা বোঝানো হলে

$P(N(D)=k)=\frac{(\lambda |D|{)}^k{e}^{-\lambda |D|}}{k!}.$

বিজ্ঞানে অন্যান্য প্রয়োগ

পয়সোঁ প্রক্রিয়ায় পর্যবেক্ষণকৃত ঘটনার সংখ্যা λ থেকে কম-বেশি হয়। আর পরিমিত ব্যবধান হয় ${\sigma }_k=\sqrt{\lambda }$। এই কম-বেশি হওয়াকে বলে পয়সোঁ নয়েজ বা (বিশেষ করে ইলেকট্রনিক্সে) শট নয়েজ।

স্বাধীন বিচ্ছিন্ন সংঘটন পরিমাপের ক্ষেত্রে গড় ও পরিমিত ব্যবধানের সংশ্লেষণ বৈজ্ঞানিক দৃষ্টিকোণ থেকে খুব কার্যকর। গড় সঙ্কেতের আশেপাশে কতটা ওঠা-নামা বা স্পন্দন হয় সেটা লক্ষ করে একটিমাত্র ঘটনার প্রভাব পরিমাপ করা যায়, যদিও সেই প্রভাব সরাসরি লক্ষ করার মতো যথেষ্ট বেশি নাও হয়। যেমন তড়িৎ প্রবাহ ও এর শট নয়েজের সংশ্লেষণ কাজে লাগিয়ে একটি ইলেকট্রনের আধান e পরিমাপ করা যায়। কোনো নির্দিষ্ট t সময়ে N সংখ্যক ইলেকট্রন একটি বিন্দুকে অতিক্রম করলে গড় প্রবাহ হবে $I=eN/t$। যেহেতু প্রবাহের ওঠা-নামা হবে ${\sigma }_I=e\sqrt{N}/t$ (যা পয়সোঁ প্রক্রিয়ার পরিমিত ব্যবধান) ক্রমের, অতএব আধান $e$ $t{\sigma }_I^2/I$ অনুপাত থেকে পরিমাপ করা যাবে।

একটি সাধারণ উদাহরণ হলো কোনো আলোকচিত্রকে বড় করা হলে যে কণা-প্রবণতা চোখে পড়ে। এ কণা-প্রবণতার কারণ কণারা নিজেরা নয়, বরং কারণ হলো রূপার কণায় পয়সোঁ স্পন্দন কমে যাওয়া। কণা-প্রবণতা ও প্রসারণের মাত্রার সংশ্লেষণ ব্যবহার করে প্রতিটি কণার প্রভাব পরিমাপ করা সম্ভব (এত ছোট এ প্রভাব খালি চোখে দেখা যায় না)। পয়সোঁ নয়েজের আরও অনেক আণবিক প্রয়োগ আবিষ্কৃত হয়েছে। যেমন, কোষ ঝিল্লিতে গ্রাহক অণুর সংখ্যা ঘনত্ব।

$\mathrm{Pr}(N_t=k)=f(k;\lambda t)=\frac{e^{-\lambda t}(\lambda t{)}^k}{k!}.$

কার্যকারণ সেট তত্ত্বে স্থানকালের বিচ্ছিন্ন উপাদান আয়তনের মধ্যে পয়সোঁ বিন্যাস মেনে চলে।

পয়সোঁ দৈব চলক উৎপাদন

জনাব নুথ পয়সোঁ বিন্যাস থেকে সংখ্যা উৎপাদনের (ছদ্ম-দৈব সংখ্যা নমুনায়ন) একটি সরল অ্যালগোরিদম প্রদান করেছেন:

algorithm poisson random number (Knuth):
    init:
         Let L ← e−λ, k ← 0 and p ← 1.
    do:
         k ← k + 1.
         Generate uniform random number u in [0,1] and let p ← p × u.
    while p > L.
    return k − 1.

প্রাপ্ত মান k এর জটিলতা রৈখিক, যার গড় মান λ। একে আরও উন্নত করে অনেকগুলো অ্যালগোরিদম তৈরি করা হয়েছে।

λ বড় হলে L = e^−λ এর মান অনেক বেশি ছোট হয়ে যায়। অ্যালগোরিদমে সামান্য পরিবর্তন এনে এ সমস্যার সমাধান করা যায়। এ জন্য নতুন একটি পরামিতি STEP নিয়ে আসা হয় যাতে e^−STEP এর মান আগের মতো ছোট হয়ে যায় না।

algorithm poisson random number (Junhao, based on Knuth):
    init:
         Let λLeft ← λ, k ← 0 and p ← 1.
    do:
         k ← k + 1.
         Generate uniform random number u in (0,1) and let p ← p × u.
         while p < 1 and λLeft > 0:
              if λLeft > STEP:
                   p ← p × eSTEP
                   λLeft ← λLeft − STEP
              else:
                   p ← p × eλLeft
                   λLeft ← 0
    while p > 1.
    return k − 1.

STEP এর মান নির্ভর করে প্রাথমিকভাবে সর্বোচ্চ কত বড় মান বাছাই করা হবে তার ওপর। ডাবল-প্রিসিশন ফ্লোটিং পয়েন্ট ফরম্যাটের ক্ষেত্রে প্রাথমিক মান e⁷⁰⁰ এর কাছাকাছি। অতএব STEP এর মান ৫০০ নেওয়া নিরাপদ।

λ এর বড় মানের ক্ষেত্রে অন্য সমাধানের মধ্যে রয়েছে প্রত্যাখ্যান নমুনায়ন ও গাউসীয় আসন্নীকরণ।

λ এর ছোট মানের ক্ষেত্রে বিপরীত রূপান্তর নমুনায়ন খুব সরল ও কার্যকর। এক্ষেত্রে প্রতিটি নমুনার জন্য শুধু একটি করে সুষম দৈব সংখ্যা u প্রয়োজন হয়। ক্রমযোজিত সম্ভাবনা u এর বেশি হওয়া পর্যন্ত নেওয়া হতে থাকে।

algorithm Poisson generator based upon the inversion by sequential search:
    init:
         Let x ← 0, p ← e−λ, s ← p.
         Generate uniform random number u in [0,1].
    while u > s do:
         x ← x + 1.
         p ← p * λ/x.
         s ← s + p.
    return x.

পরামিতির পরিমাপ

সর্বোচ্চ সম্ভাব্যতা

i = 1, ..., n এর জন্য n সংখ্যক পরিমাপকৃত মান $k_i\in \{0,1,...\}$ দেওয়া থাকলে আমরা যে পয়সোঁ সমগ্রক থেকে নমুনা নেওয়া হয়েছে তার λ পরামিতির মান পরিমাপ করতে পারব। সর্বোচ্চ সম্ভাব্যতা পরিমাপ হলো

${\hat{\lambda }}_{\mathrm{M}\mathrm{L}\mathrm{E}}=\frac{1}{n}\sum _{i=1}^n{k}_i.$

প্রতিটি মানের প্রত্যাশিত মান λ হওয়ায় এই গড়ের প্রত্যাশিত মানও λ। ফলে সর্বোচ্চ সম্ভাব্যতা পরিমাপ λ এর একটি নিরপেক্ষ পরিমাপক হবে। এছাড়াও এটি হবে একটি সূক্ষ্ম পরিমাপক। অর্থাৎ, এর পরিমাপকৃত ভেদাঙ্ক ক্র্যামার-রাও নিম্ন সীমা (CRLB) অর্জন করবে। অতএব, এটি হবে ন্যূনতম ভেদাঙ্কের নিরপেক্ষ পরিমাপ। এছাড়াও দেখানো যাবে যে এর সমষ্টি (এবং নমুনা গড়, কারণ এটি সমষ্টির এক-এক অপেক্ষক) λ এর একটি পূর্ণাঙ্গ ও পর্যাপ্ত পরিসংখ্যা।

পর্যাপ্ততা প্রমাণের জন্য আমরা উৎপাদকায়ন উপপাদ্য ব্যবহার করতে পারি। এজন্য আমরা নমুনার যুক্ত পয়সোঁ বিন্যাসের সম্ভাবনা ভর অপেক্ষককে দুটি অংশ ভাগ করি: একটি অংশ কেবলমাত্র নমুনা $𝐱$ এর ওপর নির্ভর করে (যাকে $h(𝐱)$ বলা হয়), আর আরেকটি নির্ভর করে λ ও $T(𝐱)$ অপেক্ষকের মাধ্যমে নমুনা $𝐱$ এর ওপর। অতএব $T(𝐱)$ হবে λ এর পর্যাপ্ত পরিসংখ্যা।

$P(𝐱)=\prod _{i=1}^n\frac{{\lambda }^{x_i}{e}^{-\lambda }}{x_i!}=\frac{1}{\prod _{i=1}^n{x}_i!}\times {\lambda }^{\sum _{i=1}^n{x}_i}{e}^{-n\lambda }$

এখানে প্রথম পদ $h(𝐱)$ শুধু $𝐱$ এর ওপর নির্ভর করে। দ্বিতীয় পদ $g(T(𝐱)|\lambda )$ শুধু $T(𝐱)=\sum _{i=1}^n{x}_i$ এর মাধ্যমে নমুনার ওপর নির্ভর করে। অতএব, $T(𝐱)$ পর্যাপ্ত।

সম্ভাবনা অপেক্ষককে সর্বোচ্চ মান প্রদানকারী λ এর মান বের করতে আমরা সম্ভাব্যতা অপেক্ষকের অ্যালগোরিদম ব্যবহার করতে পারি:

$\begin{array}{cc}\ell (\lambda )& =\ln \prod _{i=1}^{n}f(k_i\mid \lambda )\\ & =\sum _{i=1}^n\ln \left(\frac{e^{-\lambda }{\lambda }^{k_i}}{k_i!}\right)\\ & =-n\lambda +\left(\sum _{i=1}^n{k}_i\right)\ln (\lambda )-\sum _{i=1}^n\ln (k_i!).\end{array}$

$\ell$ কে λ এর সাপেক্ষে অন্তরীকরণ করে ০ এর সাথে তুলনা করি:

$\frac{\mathrm{d}}{\mathrm{d}\lambda }\ell (\lambda )=0⟺-n+\left(\sum _{i=1}^n{k}_i\right)\frac{1}{\lambda }=0.$

λ এর জন্য এখান থেকে একটি স্থির বিন্দু পাওয়া যায়

$\lambda =\frac{\sum _{i=1}^n{k}_i}{n}$

অতএব, k_i মানগুলোর গড় হলো λ। স্থির বিন্দুতে L এর দ্বিতীয় অন্তরকের চিহ্ন দেখে জানা যাবে λ কেমন চরম মান।

$\frac{{\partial }^2\ell }{\partial {\lambda }^2}=-{\lambda }^{-2}\sum _{i=1}^n{k}_i$

স্থির বিন্দুতে দ্বিতীয় অন্তরকের মান নির্ণয় করলে হবে:

$\frac{{\partial }^2\ell }{\partial {\lambda }^2}=-\frac{n^2}{\sum _{i=1}^n{k}_i}$

যা k_i মানগুলোর গড়ের বিপরীত সংখ্যার n গুণের ঋণাত্মক সংখ্যা। গড় ধনাত্মক হলে এটি ঋণাত্মক হয়। এই শর্ত পূরণ হলে সম্ভাবনা অপেক্ষক এই স্থির বিন্দুতে সর্বোচ্চ মান প্রদান করে।

অন্য দিকে, একটি বিন্যাস গুচ্ছকে পূর্ণাঙ্গ বলা হয় যদি এবং কেবল যদি সকল $\lambda$ এর জন্য $E(g(T))=0$ থেকে $P_{\lambda }(g(T)=0)=1$ হয়। যদি $X_i$ গুলো সুষম ও স্বাধীন বিন্যাস $\mathrm{P}\mathrm{o}(\lambda )$ মেনে চলে, তাহলে $T(𝐱)=\sum _{i=1}^n{X}_i\sim \mathrm{P}\mathrm{o}(n\lambda )$। আমাদের কাঙ্ক্ষিত বিন্যাস জানা থাকায় সহজেই দেখায় যায়, এই পরিসংখ্যা পূর্ণাঙ্গ।

$E(g(T))=\sum _{t=0}^{\mathrm{\infty }}g(t)\frac{(n\lambda {)}^t{e}^{-n\lambda }}{t!}=0$

এই সমীকরণকে সত্য হতে হলে $g(t)$ এর মান ০ হতে হবে। কারণ সমষ্টির সকল $t$ ও $\lambda$ এর সকল সম্ভাব্য মানের জন্য অন্য কোনো পদই ০ হবে না। অতএব, সকল $\lambda$ এর জন্য $E(g(T))=0$ থেকে $P_{\lambda }(g(T)=0)=1$। অতএব, পরিসংখ্যাটিকে পূর্ণাঙ্গ হিসেবে প্রমাণ করা গেল।

আস্থা ব্যাপ্তি

পয়সোঁ ও কাই-বর্গ বিন্যাসের ক্রমযোজিত বিন্যাস অপেক্ষকের সম্পর্কের মাধ্যমে পয়সোঁ বিন্যাসের গড়ের আস্থা ব্যাপ্তি বের করা যায়। কাই-বর্গ বিন্যাসের সাথে আবার গ্যামা বিন্যাসের নিবিড় সম্পর্ক রয়েছে। এটা থেকে একটি বিকল্প প্রকাশও পাওয়া যায়। μ গড় বিশিষ্ট পয়সোঁ বিন্যাসের একটি মান k দেওয়া থাকলে টেমপ্লেট:Math আস্থা স্তরে μ এর আস্থা ব্যাপ্তি হবে

$\frac{1}{2}{\chi }^2(\alpha /2;2k)\le \mu \le \frac{1}{2}{\chi }^2(1-\alpha /2;2k+2),$

অথবা সমতুল্যরূপে,

$F^{-1}(\alpha /2;k,1)\le \mu \le F^{-1}(1-\alpha /2;k+1,1),$

যেখানে ${\chi }^2(p;n)$ হলো n স্বাধীনতার মাত্রার কাই-বর্গ বিন্যাসের কোয়ান্টাইল অপেক্ষক (p এর নিম্ন প্রান্তে), আর $F^{-1}(p;n,1)$ হলো গ্যামা বিন্যাসের কোয়ান্টাইল অপেক্ষক যেখানে আকৃতি পরামিতি n ও মাপনী পরামিতি ১।টেমপ্লেট:R এই ব্যাপ্তিটি এই অর্থে প্রকৃত যে এর কাভারেজ সম্ভবনা কখনও টেমপ্লেট:Math এর কম হয় না।

গ্যামা বিন্যাসের কোয়ান্টাইল পাওয়া না গেলেও (উইলসন-হিলফেরটি রূপান্তরের ভিত্তিতে) প্রকৃত ব্যাপ্তির একটি নিখুঁত আসন্ন মান প্রস্তাব করা হয়েছে:

$k{(1-\frac{1}{9k}-\frac{z_{\alpha /2}}{3\sqrt{k}})}^3\le \mu \le (k+1){(1-\frac{1}{9(k+1)}+\frac{z_{\alpha /2}}{3\sqrt{k+1}})}^3,$

যেখানে $z_{\alpha /2}$ হলো আদর্শ পরিমিত বিন্যাসের মান, যার উর্ধ অংশের ক্ষেত্রফল টেমপ্লেট:Math।

উপরের মতো একই প্রসঙ্গে (λ গড়ের পয়সোঁ বিন্যাসের n সংখ্যক k_i মানের একটি নমুনা দেওয়া থাকলে) এই সূত্রগুলোর প্রয়োগের জন্য নিচের সমীকরণটি কাজে লাগানো হয়:

$k=\sum _{i=1}^n{k}_i,$

এখান থেকে μ = nλ এর ব্যাপ্তি হিসেব করে বের করা হয় এবং পরিশেষে λ এর ব্যাপ্তি বের করা হয়।

বায়েসীয় অনুমিতি

বায়েসীয় অনুমিতিতে পয়সোঁ বিন্যাসের হার পরামিতি λ এর কনজুগেট প্রাইঅর হলো গ্যামা বিন্যাস। ধরা যাক

$\lambda \sim \mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(\alpha ,\beta )$

দিয়ে বোঝানো হচ্ছে λ একটি গ্যামা সম্ভাবনা বিন্যাস g মেনে চলছে, যার আকৃতি পরামিতি α ও বিপরীত মাপনী পরামিতি β:

$g(\lambda \mid \alpha ,\beta )=\frac{{\beta }^{\alpha }}{\mathrm{\Gamma }(\alpha )}{\lambda }^{\alpha -1}{e}^{-\beta \lambda }\text{ for }\lambda >0.$

তাহলে আগের মতো n আকারের পরিমাপকৃত k_i এর একটি নমুনা ও একটি Gamma(α, β) প্রাইঅর দেওয়া থাকলে পোস্টেরিয়র বা উত্তর হবে:

$\lambda \sim \mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(\alpha +\sum _{i=1}^n{k}_i,\beta +n).$

পোস্টেরিয়র গড় E[λ] সর্বোচ্চ সম্ভাব্যতা পরিমাপ ${\hat{\lambda }}_{\mathrm{M}\mathrm{L}\mathrm{E}}$ এর কাছাকাছি হবে $\alpha \to 0,\beta \to 0$ সীমায়, যা গ্যামা বিন্যাসের গড়ের সাধারণ প্রকাশ থেকে সাথে সাথেই পাওয়া যায়।

একটি বাড়তি মানের জন্য উত্তর পূর্বাভাসমূলক বিন্যাস হবে ঋণাত্মক দ্বিপদী বিন্যাস, যাকে অনেকসময় গ্যামা-পয়সোঁ বিন্যাসও বলা হয়।

বহু পয়সোঁ গড়ের যুগপৎ পরিমাপ

ধরা যাক, $X_1,X_2,\dots ,X_p$ হলো $p$ সংখ্যক পয়সোঁ বিন্যাসের এক গুচ্ছ স্বাধীন দৈব চলক, যাদের প্রতিটির পরামিতি ${\lambda }_i$, যেখানে $i=1,\dots ,p$। আমরা এই পরামিতিগুলো পরিমাপ করতে চাই। ক্লিভেনসন ও জিডেক দেখিয়েছেন যে পরমিতকৃত বর্গ ত্রুটি ক্ষয়ের অধীনে $L(\lambda ,\hat{\lambda })=\sum _{i=1}^p{\lambda }_i^{-1}({\hat{\lambda }}_i-{\lambda }_i{)}^2$, যেখানে $p>1$। তাহলে পরিমিত গড়ের ক্ষেত্রে স্টেইনের উদাহরণের মতোই সর্বোচ্চ সম্ভাব্যতা পরিমাপ (MLE) ${\hat{\lambda }}_i=X_i$ অগ্রহণযোগ্য।

এক্ষেত্রে যে-কোনো $0<c\le 2(p-1)$ ও $b\ge (p-2+p^{-1})$ এর জন্য গুরুলঘু পরিমাপক গুচ্ছ হবে

${\hat{\lambda }}_i=(1-\frac{c}{b+\sum _{i=1}^p{X}_i})X_i,i=1,\dots ,p.$

দ্বিচলক পয়সোঁ বিন্যাস

এই বিন্যাসটিকে সম্প্রসারিত করে দুটি চলকের জন্যও উপযোগী করা হয়েছে। বিন্যাসটির উৎপাদী অপেক্ষক হলো

$g(u,v)=\exp [({\theta }_1-{\theta }_{12})(u-1)+({\theta }_2-{\theta }_{12})(v-1)+{\theta }_{12}(uv-1)]$

যেখানে

${\theta }_1,{\theta }_2>{\theta }_{12}>0$

প্রান্তিক বিন্যাসগুলো হলো Poisson(θ₁) ও Poisson(θ₂) এবং সংশ্লেষণাঙ্ক

$0\le \rho \le \min \{\frac{{\theta }_1}{{\theta }_2},\frac{{\theta }_2}{{\theta }_1}\}$

পরিসরে আবদ্ধ।

দ্বিচলক পয়সোঁ বিন্যাস উৎপাদন করার একটি সরল প্রক্রিয়া হলো তিনটি স্বাধীন পয়সোঁ বিন্যাস $Y_1,Y_2,Y_3$ নেওয়া, যাদের গড় যথাক্রমে ${\lambda }_1,{\lambda }_2,{\lambda }_3$। অতঃপর $X_1=Y_1+Y_3,X_2=Y_2+Y_3$ বসাতে হবে। দ্বিচলক পয়সোঁ বিন্যাসের সম্ভাবনা অপেক্ষক হবে

$\begin{array}{cc}& \mathrm{Pr}(X_1=k_1,X_2=k_2)\\ =& \exp (-{\lambda }_1-{\lambda }_2-{\lambda }_3)\frac{{\lambda }_1^{k_1}}{k_1!}\frac{{\lambda }_2^{k_2}}{k_2!}\sum _{k=0}^{\min (k_1,k_2)}(\genfrac{}{}{0ex}{}{k_1}{k})(\genfrac{}{}{0ex}{}{k_2}{k})k!{\left(\frac{{\lambda }_3}{{\lambda }_1{\lambda }_2}\right)}^k\end{array}$

পয়সোঁ বিন্যাসের জন্য কম্পিউটার সফটওয়্যার

নিবেদিত সফটওয়্যার লাইব্রেরিগুলো দিয়ে পয়সোঁ বিন্যাসের মূলত দুটি কাজ করা হয়: পয়সোঁ বিন্যাস $P(k;\lambda )$ এর মান নির্ণয় ও এ বিন্যাস থেকে দৈব চলক উৎপাদন।

পয়সোঁ বিন্যাসের মান নির্ণয়

নির্দিষ্ট $k$ ও $\lambda$ এর জন্য $P(k;\lambda )$ এর মান সহজেই বের করা যায়। এটা করা হয় সূচকীয়, ঘাত ও ফ্যাক্টোরিয়াল অপেক্ষকের ভিত্তিতে $P(k;\lambda )$ এর আদর্শ সংজ্ঞা কাজে লাগিয়ে। অবশ্য $k$ এর বড় মানের ক্ষেত্রে নাকচ হয়ে যাওয়ার ঝুঁকি থেকেই যায়। এ ঝুঁকি এড়াতে আদর্শ লাইব্রেরি math.h এর টেমপ্লেট:Anchorlgamma ব্যবহার করা যেতে পারে।

কিছু কম্পিউটিং লাইব্রেরিতে সহজাতভাবেই পয়সোঁ বিন্যাস থাকে। যেমন:

• R প্রোগ্রামিং ল্যাংগুয়েজ: ফাংশন dpois(x, lambda);
• এক্সেল: ফাংশন POISSON( x, mean, cumulative), যাতে ক্রমযোজিত বিন্যাস উল্লেখ করারও ব্যবস্থা আছে;
• ম্যাথম্যাটিকা: এক চলকের পয়সোঁ বিন্যাস PoissonDistribution[$\lambda$], দ্বিচলক পয়সোঁ বিন্যাস MultivariatePoissonDistribution[${\theta }_{12}$,{ ${\theta }_1-{\theta }_{12}$, ${\theta }_2-{\theta }_{12}$}],.

পয়সোঁ বিন্যাস থেকে দৈব চয়ন

নির্দিষ্ট $\lambda$ এর জন্য পয়সোঁ বিন্যাস থেকে দৈব পূর্ণ সংখ্যা উৎপাদন করা আরও সহজ। এ কাজের জন্য আছে:

• জিএনইউ সায়েন্টিফিক লাইব্রেরি (GSL): ফাংশন gsl_ran_poisson
• R প্রোগ্রামিং ল্যাংগুয়েজ: ফাংশন rpois(n, lambda);

আরও দেখুন

• সম্ভাবনা বিন্যাস
• সম্ভাবনা তত্ত্ব
• সম্ভাবনার বিধিমালা
• পরিসংখ্যান