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

อ่าน 27 นาที

การแจกแจงแบบเอกซ์โปเนนเชียล

ในทฤษฎีความน่าจะเป็นและสถิติการแจกแจงเอกซ์โพเนนเชียลหรือการแจกแจงเอกซ์โพเนนเชียลเชิงลบคือการแจกแจงความน่าจะเป็นของระยะห่างระหว่างเหตุการณ์ในกระบวนการจุดปัวซง กล่าวคือ

การแจกแจงแบบเอกซ์โปเนนเชียล

เลขชี้กำลัง
ฟังก์ชันความหนาแน่นความน่าจะเป็น
กราฟแสดงฟังก์ชันความหนาแน่นความน่าจะเป็นของการแจกแจงแบบเอกซ์โปเนนเชียล
ฟังก์ชันการกระจายสะสม
ฟังก์ชันการกระจายสะสม
พารามิเตอร์λ>0,{\displaystyle \lambda >0,}อัตรา หรือมาตราส่วน ผกผัน
สนับสนุนx[0,){\displaystyle x\in [0,\infty )}
พีดีλอีλx{\displaystyle \lambda e^{-\lambda x}}
ซีดีเอฟ1อีλx{\displaystyle 1-e^{-\lambda x}}
ควอนไทล์ln(1พี)λ{\displaystyle -{\frac {\ln(1-p)}{\แลมบ์ดา }}}
หมายถึง1λ{\displaystyle {\frac {1}{\lambda }}}
ค่ามัธยฐานln2λ{\displaystyle {\frac {\ln 2}{\lambda }}}
โหมด0{\displaystyle 0}
ความแปรปรวน1λ2{\displaystyle {\frac {1}{\lambda ^{2}}}}
ความเบี่ยงเบน2{\displaystyle 2}
ความโค้งส่วนเกิน6{\displaystyle 6}
เอนโทรปี1lnλ{\displaystyle 1-\ln \lambda }
เอ็มจีเอฟλλที, สำหรับ ที<λ{\displaystyle {\frac {\lambda }{\lambda -t}},{\text{ for }}t<\lambda }
ซีเอฟλλฉันที{\displaystyle {\frac {\lambda }{\lambda -it}}}
ข้อมูลของฟิชเชอร์1λ2{\displaystyle {\frac {1}{\lambda ^{2}}}}
ความแตกต่าง Kullback–Leiblerlnλ0λ+λλ01{\displaystyle \ln {\frac {\lambda _{0}}{\lambda }}+{\frac {\lambda }{\lambda _{0}}}-1}

ในทฤษฎีความน่าจะเป็นและสถิติการแจกแจงเอกซ์โพเนนเชียลหรือการแจกแจงเอกซ์โพเนนเชียลเชิงลบคือการแจกแจงความน่าจะเป็นของระยะห่างระหว่างเหตุการณ์ในกระบวนการจุดปัวซง กล่าวคือ กระบวนการที่เหตุการณ์เกิดขึ้นอย่างต่อเนื่องและเป็นอิสระในอัตราเฉลี่ยคงที่ พารามิเตอร์ระยะห่างอาจเป็นการวัดมิติเดียวที่มีความหมายของกระบวนการ เช่น เวลาระหว่างข้อผิดพลาดในการผลิต หรือความยาวตามม้วนผ้าในกระบวนการผลิตการทอ[ 1 ]เป็นการแจกแจงแกมมา แบบเฉพาะเจาะจง เป็นอนาล็อกแบบต่อเนื่องของการแจกแจงเรขาคณิตและมีคุณสมบัติสำคัญคือไม่มีความจำ[ 2 ]นอกจากการใช้สำหรับการวิเคราะห์กระบวนการจุดปัวซงแล้ว ยังพบได้ในบริบทอื่นๆ อีกมากมาย[ 3 ]

การแจกแจงแบบเอกซ์โพเนนเชียลไม่เหมือนกับกลุ่ม การแจกแจง แบบตระกูลเอกซ์โพเนนเชีย ล ซึ่งเป็นกลุ่มการแจกแจงความน่าจะเป็นขนาดใหญ่ที่รวมการแจกแจงแบบเอกซ์โพเน นเชียลไว้เป็นหนึ่งในสมาชิก แต่ยังรวมถึงการแจกแจงอื่นๆ อีกมากมาย เช่นการแจกแจงแบบปกติแบบทวินาม แบบแกมมา และแบบปัวซง[ 3 ]

คำจำกัดความ

ฟังก์ชันความหนาแน่นความน่าจะเป็น

ฟังก์ชันความหนาแน่นความน่าจะเป็น (pdf) ของการแจกแจงแบบเอกซ์โปเนนเชียลคือ

เอฟ(x;λ)={λอีλxx0,0x<0.{\displaystyle f(x;\lambda )={\begin{cases}\lambda e^{-\lambda x}&x\geq 0,\\0&x<0.\end{cases}}}

ในที่นี้λ > 0 คือพารามิเตอร์ของการแจกแจง ซึ่งมักเรียกว่าพารามิเตอร์อัตราการแจกแจงนี้ครอบคลุมช่วง[ 0, ∞)หากตัวแปรสุ่มXมีการแจกแจงนี้ เราจะเขียนว่าX ~ Exp( λ )  

การแจกแจงแบบเอกซ์โปเนนเชียลแสดงให้เห็นถึงการหารลงตัวอย่างไม่มีที่สิ้นสุด

ฟังก์ชันการกระจายสะสม

ฟังก์ชันการกระจายสะสมกำหนดโดย

เอฟ(x;λ)={1อีλxx0,0x<0.{\displaystyle F(x;\lambda )={\begin{cases}1-e^{-\lambda x}&x\geq 0,\\0&x<0.\end{cases}}}

การกำหนดพารามิเตอร์ทางเลือก

การแจกแจงแบบเอกซ์โปเนนเชียลบางครั้งจะถูกกำหนดพารามิเตอร์โดยใช้พารามิเตอร์มาตราส่วนβ = 1/ λซึ่งก็คือค่าเฉลี่ยด้วยเช่นกัน: เอฟ(x;เบต้า)={1เบต้าอีx/เบต้าx0,0x<0.เอฟ(x;เบต้า)={1อีx/เบต้าx0,0x<0.{\displaystyle f(x;\beta )={\begin{cases}{\frac {1}{\beta }}e^{-x/\beta }&x\geq 0,\\0&x<0.\end{cases}}\qquad \qquad F(x;\beta )={\begin{cases}1-e^{-x/\beta }&x\geq 0,\\0&x<0.\end{cases}}}

คุณสมบัติ

ค่าเฉลี่ย ความแปรปรวน โมเมนต์ และมัธยฐาน

ค่าเฉลี่ยคือจุดศูนย์กลางมวลความน่าจะเป็น ซึ่งก็คือโมเมนต์แรก
ค่ามัธยฐานคือภาพก่อนหน้าF −1 (1/2)

ค่าเฉลี่ยหรือค่าคาดหวังของตัวแปรสุ่มX ที่มีการแจกแจงแบบเอกซ์โปเนนเชียล โดยมีพารามิเตอร์อัตราλจะกำหนดโดย อี[X]=1λ.{\displaystyle \operatorname {E} [X]={\frac {1}{\lambda }}.}

จากตัวอย่างที่ยกมาด้านล่างนี้ ทำให้เข้าใจได้ คนที่ได้รับโทรศัพท์โดยเฉลี่ยสองสายต่อชั่วโมง สามารถคาดหวังได้ว่าช่วงเวลาระหว่างการโทรแต่ละครั้งจะอยู่ที่ 0.5 ชั่วโมง หรือ 30 นาที

ค่าความแปรปรวนของXกำหนดโดย วาร์[X]=1λ2,{\displaystyle \operatorname {Var} [X]={\frac {1}{\lambda ^{2}}},} ดังนั้นค่าเบี่ยงเบนมาตรฐานจึงเท่ากับค่าเฉลี่ย

ช่วงเวลาของXสำหรับnเอ็น{\displaystyle n\in \mathbb {N} }ได้รับจาก อี[Xn]=n!λn.{\displaystyle \operatorname {E} \left[X^{n}\right]={\frac {n!}{\lambda ^{n}}}.}

ช่วงเวลาสำคัญของXสำหรับnเอ็น{\displaystyle n\in \mathbb {N} }ได้รับจาก μn=!nλn=n!λnเค=0n(1)เคเค!.{\displaystyle \mu _{n}={\frac {!n}{\lambda ^{n}}}={\frac {n!}{\lambda ^{n}}}\sum _{k=0}^{n}{\frac {(-1)^{k}}{k!}}.} โดยที่ ! nคือ แฟกทอเรี ยลย่อยของn

ค่ามัธยฐานของXหาได้จากสูตร [X]=ln(2)λ<อี[X],{\displaystyle \operatorname {m} [X]={\frac {\ln(2)}{\lambda }}<\operatorname {E} [X],} โดยที่lnหมายถึงลอการิทึมธรรมชาติดังนั้นผลต่างสัมบูรณ์ระหว่างค่าเฉลี่ยและค่ามัธยฐานคือ |อี[X][X]|=1ln(2)λ<1λ=σ[X],{\displaystyle \left|\operatorname {E} \left[X\right]-\operatorname {m} \left[X\right]\right|={\frac {1-\ln(2)}{\lambda }}<{\frac {1}{\lambda }}=\operatorname {\sigma } [X],}

ตามความไม่เท่าเทียมกันระหว่างค่ามัธยฐานและค่าเฉลี่ย

คุณสมบัติไร้ความทรงจำของตัวแปรสุ่มเอกซ์โพเนนเชียล

ตัวแปรสุ่มT ที่มีการแจกแจงแบบเอกซ์โปเนนเชียล เป็นไปตามความสัมพันธ์ ดังต่อไปนี้ปร.(ที>+ทีที>)=ปร.(ที>ที),,ที0.{\displaystyle \Pr \left(T>s+t\mid T>s\right)=\Pr(T>t),\qquad \forall s,t\geq 0.}

สามารถเห็นได้จากการพิจารณาฟังก์ชันการกระจายสะสมเสริม : ปร.(ที>+ทีที>)=ปร.(ที>+ทีที>)ปร.(ที>)=ปร.(ที>+ที)ปร.(ที>)=อีλ(+ที)อีλ=อีλที=ปร.(ที>ที).{\displaystyle {\begin{aligned}\Pr \left(T>s+t\mid T>s\right)&={\frac {\Pr \left(T>s+t\cap T>s\right)}{\Pr \left(T>s\right)}}\\[4pt]&={\frac {\Pr \left(T>s+t\right)}{\Pr \left(T>s\right)}}\\[4pt]&={\frac {e^{-\lambda (s+t)}}{e^{-\lambda s}}}\\[4pt]&=e^{-\lambda t}\\[4pt]&=\Pr(T>t).\end{aligned}}}

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

การแจกแจงแบบเอกซ์โปเนนเชียลและการแจกแจงแบบเรขาคณิตเป็นเพียงการแจกแจงความน่าจะเป็นที่ไม่มีความจำเพียงสองแบบเท่านั้น

ดังนั้น การแจกแจงแบบเอกซ์โปเนนเชียลจึงเป็นการแจกแจงความน่าจะเป็นแบบต่อเนื่องเพียงอย่างเดียวที่มีอัตราความล้มเหลว คงที่ด้วยเช่น กัน

ควอนไทล์

เกณฑ์ความผิดปกติของ Tukey สำหรับฟังก์ชันการกระจายความน่าจะเป็นแบบเอกซ์โปเนนเชียล
เกณฑ์ของ Tukey สำหรับความผิดปกติ

ฟังก์ชันควอนไทล์ (ฟังก์ชันการกระจายสะสมผกผัน) สำหรับ Exp( λ ) คือ เอฟ1(พี;λ)=ln(1พี)λ,0พี<1{\displaystyle F^{-1}(p;\lambda )={\frac {-\ln(1-p)}{\lambda }},\qquad 0\leq p<1}

ดังนั้น ควาร์ไทล์จึงเป็นดังนี้:

และผลที่ตามมาคือช่วงควาร์ไทล์คือln(3)/ λ

มูลค่าความเสี่ยงแบบมีเงื่อนไข (การขาดทุนที่คาดหวัง)

ค่าความเสี่ยงตามเงื่อนไข (CVaR) หรือที่รู้จักกันในชื่อการขาดทุนที่คาดหวังหรือซูเปอร์ควอนไทล์สำหรับ Exp( λ ) ได้มาดังนี้: [ 4 ]

q¯α(X)=11αα1qพี(X)พี=1(1α)α1ln(1พี)λพี=1λ(1α)1α0ln(y)y=1λ(1α)01αln(y)y=1λ(1α)[(1α)ln(1α)(1α)]=ln(1α)+1λ{\displaystyle {\begin{aligned}{\bar {q}}_{\alpha }(X)&={\frac {1}{1-\alpha }}\int _{\alpha }^{1}q_{p}(X)dp\\&={\frac {1}{(1-\alpha )}}\int _{\alpha }^{1}{\frac {-\ln(1-p)}{\lambda }}dp\\&={\frac {-1}{\lambda (1-\alpha )}}\int _{1-\alpha }^{0}-\ln(y)dy\\&={\frac {-1}{\lambda (1-\alpha )}}\int _{0}^{1-\alpha }\ln(y)dy\\&={\frac {-1}{\lambda (1-\alpha )}}[(1-\alpha )\ln(1-\alpha )-(1-\alpha )]\\&={\frac {-\ln(1-\alpha )+1}{\lambda }}\\\end{aligned}}}

ความน่าจะเป็นเกินค่าที่ปรับลดแล้ว (bPOE)

ความน่าจะเป็นที่เกินค่าที่ปรับแก้แล้วคือหนึ่งลบด้วยระดับความน่าจะเป็นที่ค่า CVaR เท่ากับค่าเกณฑ์x{\displaystyle x}ได้มาดังนี้: [ 4 ]

พี¯x(X)={1α|q¯α(X)=x}={1α|ln(1α)+1λ=x}={1α|ln(1α)=1λx}={1α|อีln(1α)=อี1λx}={1α|1α=อี1λx}=อี1λx{\displaystyle {\begin{aligned}{\bar {p}}_{x}(X)&=\{1-\alpha |{\bar {q}}_{\alpha }(X)=x\}\\&=\{1-\alpha |{\frac {-\ln(1-\alpha )+1}{\lambda }}=x\}\\&=\{1-\alpha |\ln(1-\alpha )=1-\lambda x\}\\&=\{1-\alpha |e^{\ln(1-\alpha )}=e^{1-\lambda x}\}=\{1-\alpha |1-\alpha =e^{1-\lambda x}\}=e^{1-\lambda x}\end{aligned}}}

ความแตกต่าง Kullback–Leibler

ความแตกต่าง แบบKullback–Leibler ที่มีทิศทาง ในnatsของอีλ{\displaystyle e^{\lambda }}(การแจกแจงโดยประมาณ) จากอีλ0{\displaystyle e^{\lambda _{0}}}(การกระจาย 'ที่แท้จริง') กำหนดโดย Δ(λ0λ)=อีλ0(บันทึกพีλ0(x)พีλ(x))=อีλ0(บันทึกλ0อีλ0xλอีλx)=บันทึก(λ0)บันทึก(λ)(λ0λ)อีλ0(x)=บันทึก(λ0)บันทึก(λ)+λλ01.{\displaystyle {\begin{aligned}\Delta (\lambda _{0}\parallel \lambda )&=\mathbb {E} _{\lambda _{0}}\left(\log {\frac {p_{\lambda _{0}}(x)}{p_{\lambda }(x)}}\right)\\&=\mathbb {E} _{\lambda _{0}}\left(\log {\frac {\lambda _{0}e^{\lambda _{0}x}}{\lambda e^{\lambda x}}}\right)\\&=\log(\lambda _{0})-\log(\lambda )-(\lambda _{0}-\lambda )E_{\lambda _{0}}(x)\\&=\log(\lambda _{0})-\log(\lambda )+{\frac {\lambda }{\lambda _{0}}}-1.\end{aligned}}}

การกระจายเอนโทรปีสูงสุด

ในบรรดาการแจกแจงความน่าจะเป็นต่อเนื่องทั้งหมดที่มีช่วง[ 0, ∞)และค่าเฉลี่ยμการแจกแจงเอกซ์โพเนนเชียลที่มีλ = 1/ μ มี เอนโทรปีเชิงอนุพันธ์มากที่สุดกล่าวอีกนัยหนึ่งคือเป็นการแจกแจงความน่าจะเป็นที่มีเอนโทรปีสูงสุดสำหรับตัวแปรสุ่มXซึ่งมากกว่าหรือเท่ากับศูนย์และสำหรับ E[ X ] คงที่[ 5 ]

การกระจายของค่าต่ำสุดของตัวแปรสุ่มเอกซ์โปเนนเชียล

ให้X , ..., X เป็น ตัวแปรสุ่ม อิสระที่มีการแจกแจงแบบเอกซ์โปเนนเชียล โดยมีพารามิเตอร์อัตราλ , ..., λ แล้ว นาที{X1,,Xn}{\displaystyle \min \left\{X_{1},\dotsc ,X_{n}\right\}} นอกจากนี้ยังมีการแจกแจงแบบเอกซ์โปเนนเชียลด้วย โดยมีพารามิเตอร์ λ=λ1++λn.{\displaystyle \lambda =\lambda _{1}+\dotsb +\lambda _{n}.}

สามารถเห็นได้จากการพิจารณาฟังก์ชันการกระจายสะสมเสริม : ปร.(นาที{X1,,Xn}>x)=ปร.(X1>x,,Xn>x)=ฉัน=1nปร.(Xฉัน>x)=ฉัน=1nเอ็กซ์(xλฉัน)=เอ็กซ์(xฉัน=1nλฉัน).{\displaystyle {\begin{aligned}&\Pr \left(\min\{X_{1},\dotsc ,X_{n}\}>x\right)\\={}&\Pr \left(X_{1}>x,\dotsc ,X_{n}>x\right)\\={}&\prod _{i=1}^{n}\Pr \left(X_{i}>x\right)\\={}&\prod _{i=1}^{n}\exp \left(-x\lambda _{i}\right)=\exp \left(-x\sum _{i=1}^{n}\lambda _{i}\right).\end{aligned}}}

ดัชนีของตัวแปรที่ได้ค่าต่ำสุดจะมีการกระจายตัวตามการแจกแจงเชิงหมวดหมู่ ปร.(Xเค=นาที{X1,,Xn})=λเคλ1++λn.{\displaystyle \Pr \left(X_{k}=\min\{X_{1},\dotsc ,X_{n}\}\right)={\frac {\lambda _{k}}{\lambda _{1}+\dotsb +\lambda _{n}}}.}

สามารถเห็นหลักฐานได้โดยการปล่อยให้ฉัน=อาร์กมินฉัน{1,,n}{X1,,Xn}{\displaystyle I=\operatorname {argmin} _{i\in \{1,\dotsb ,n\}}\{X_{1},\dotsc ,X_{n}\}}. แล้ว, ปร.(ฉัน=เค)=0ปร.(Xเค=x)ปร.(ฉันเคXฉัน>x)x=0λเคอีλเคx(ฉัน=1,ฉันเคnอีλฉันx)x=λเค0อี(λ1++λn)xx=λเคλ1++λn.{\displaystyle {\begin{aligned}\Pr(I=k)&=\int _{0}^{\infty }\Pr(X_{k}=x)\Pr(\forall _{i\neq k}X_{i}>x)\,dx\\&=\int _{0}^{\infty }\lambda _{k}e^{-\lambda _{k}x}\left(\prod _{i=1,i\neq k}^{n}e^{-\lambda _{i}x}\right)dx\\&=\lambda _{k}\int _{0}^{\infty }e^{-\left(\lambda _{1}+\dotsb +\lambda _{n}\right)x}dx\\&={\frac {\lambda _{k}}{\lambda _{1}+\dotsb +\lambda _{n}}}.\end{aligned}}}

โปรดทราบว่า สูงสุด{X1,,Xn}{\displaystyle \max\{X_{1},\dotsc ,X_{n}\}} ไม่มีการแจกแจงแบบเอกซ์โปเนนเชียล หากX , ..., X ไม่มีพารามิเตอร์เป็น 0 ทั้งหมด[ 6 ]

โมเมนต์ร่วมของสถิติอันดับเลขชี้กำลังแบบอิสระและมีการกระจายเหมือนกัน

อนุญาตX1,,Xn{\displaystyle X_{1},\dotsc ,X_{n}}เป็นn{\displaystyle n}ตัวแปรสุ่มเอกซ์โพเนน เชียลอิสระและมีการกระจายเหมือนกันโดยมีพารามิเตอร์อัตราλให้X(1),,X(n){\displaystyle X_{(1)},\dotsc ,X_{(n)}}แทนค่าสถิติลำดับ ที่สอดคล้องกัน สำหรับฉัน<เจ{\displaystyle i<j}โมเมนต์ร่วมอี[X(ฉัน)X(เจ)]{\displaystyle \operatorname {E} \left[X_{(i)}X_{(j)}\right]}ของสถิติลำดับX(ฉัน){\displaystyle X_{(i)}}และX(เจ){\displaystyle X_{(j)}}ได้รับจาก อี[X(ฉัน)X(เจ)]=เค=0เจ11(nเค)λอี[X(ฉัน)]+อี[X(ฉัน)2]=เค=0เจ11(nเค)λเค=0ฉัน11(nเค)λ+เค=0ฉัน11((nเค)λ)2+(เค=0ฉัน11(nเค)λ)2.{\displaystyle {\begin{aligned}\operatorname {E} \left[X_{(i)}X_{(j)}\right]&=\sum _{k=0}^{j-1}{\frac {1}{(n-k)\lambda }}\operatorname {E} \left[X_{(i)}\right]+\operatorname {E} \left[X_{(i)}^{2}\right]\\&=\sum _{k=0}^{j-1}{\frac {1}{(n-k)\lambda }}\sum _{k=0}^{i-1}{\frac {1}{(n-k)\lambda }}+\sum _{k=0}^{i-1}{\frac {1}{((n-k)\lambda )^{2}}}+\left(\sum _{k=0}^{i-1}{\frac {1}{(n-k)\lambda }}\right)^{2}.\end{aligned}}}

สิ่งนี้สามารถเห็นได้จากการอ้างอิงกฎแห่งความคาดหวังโดยรวมและคุณสมบัติไร้ความทรงจำ: อี[X(ฉัน)X(เจ)]=0อี[X(ฉัน)X(เจ)X(ฉัน)=x]เอฟX(ฉัน)(x)x=x=0xอี[X(เจ)X(เจ)x]เอฟX(ฉัน)(x)x(เนื่องจาก X(ฉัน)=xX(เจ)x)=x=0x[อี[X(เจ)]+x]เอฟX(ฉัน)(x)x(โดยทรัพย์สินที่ไร้ความทรงจำ)=เค=0เจ11(nเค)λอี[X(ฉัน)]+อี[X(ฉัน)2].{\displaystyle {\begin{aligned}\operatorname {E} \left[X_{(i)}X_{(j)}\right]&=\int _{0}^{\infty }\operatorname {E} \left[X_{(i)}X_{(j)}\mid X_{(i)}=x\right]f_{X_{(i)}}(x)\,dx\\&=\int _{x=0}^{\infty }x\operatorname {E} \left[X_{(j)}\mid X_{(j)}\geq x\right]f_{X_{(i)}}(x)\,dx&&\left({\textrm {since}}~X_{(i)}=x\implies X_{(j)}\geq x\right)\\&=\int _{x=0}^{\infty }x\left[\operatorname {E} \left[X_{(j)}\right]+x\right]f_{X_{(i)}}(x)\,dx&&\left({\text{by the memoryless property}}\right)\\&=\sum _{k=0}^{j-1}{\frac {1}{(n-k)\lambda }}\operatorname {E} \left[X_{(i)}\right]+\operatorname {E} \left[X_{(i)}^{2}\right].\end{aligned}}}

สมการแรกได้มาจากกฎของความคาดหวังโดยรวมสมการที่สองใช้ประโยชน์จากข้อเท็จจริงที่ว่า เมื่อเรากำหนดเงื่อนไขบนX(ฉัน)=x{\displaystyle X_{(i)}=x}จึงต้องสรุปได้ว่าX(เจ)x{\displaystyle X_{(j)}\geq x}สมการที่สามอาศัยคุณสมบัติไร้ความทรงจำในการแทนที่อี[X(เจ)X(เจ)x]{\displaystyle \operatorname {E} \left[X_{(j)}\mid X_{(j)}\geq x\right]}กับอี[X(เจ)]+x{\displaystyle \operatorname {E} \left[X_{(j)}\right]+x}.

ผลรวมของตัวแปรสุ่มเอกซ์โปเนนเชียลอิสระสองตัว

ฟังก์ชันการกระจายความน่าจะเป็น (PDF) ของผลรวมของตัวแปรสุ่มอิสระสองตัว คือการสังเคราะห์ (convolution) ของ PDF ของตัวแปรสุ่มแต่ละตัวถ้าX1{\displaystyle X_{1}}และX2{\displaystyle X_{2}}เป็นตัวแปรสุ่มเอกซ์โพเนนเชียลอิสระที่มีพารามิเตอร์อัตราตามลำดับλ1{\displaystyle \lambda _{1}}และλ2,{\displaystyle \lambda _{2},}จากนั้นความหนาแน่นของความน่าจะเป็นของ=X1+X2{\displaystyle Z=X_{1}+X_{2}}ได้รับจาก เอฟ(z)=เอฟX1(x1)เอฟX2(zx1)x1=0zλ1อีλ1x1λ2อีλ2(zx1)x1=λ1λ2อีλ2z0zอี(λ2λ1)x1x1={λ1λ2λ2λ1(อีλ1zอีλ2z) ถ้า λ1λ2λ2zอีλz ถ้า λ1=λ2=λ.{\displaystyle {\begin{aligned}f_{Z}(z)&=\int _{-\infty }^{\infty }f_{X_{1}}(x_{1})f_{X_{2}}(z-x_{1})\,dx_{1}\\&=\int _{0}^{z}\lambda _{1}e^{-\lambda _{1}x_{1}}\lambda _{2}e^{-\lambda _{2}(z-x_{1})}\,dx_{1}\\&=\lambda _{1}\lambda _{2}e^{-\lambda _{2}z}\int _{0}^{z}e^{(\lambda _{2}-\lambda _{1})x_{1}}\,dx_{1}\\&={\begin{cases}{\dfrac {\lambda _{1}\lambda _{2}}{\lambda _{2}-\lambda _{1}}}\left(e^{-\lambda _{1}z}-e^{-\lambda _{2}z}\right)&{\text{ if }}\lambda _{1}\neq \lambda _{2}\\[4pt]\lambda ^{2}ze^{-\lambda z}&{\text{ if }}\lambda _{1}=\lambda _{2}=\lambda .\end{cases}}\end{aligned}}} เอนโทรปีของการแจกแจงนี้สามารถหาได้ในรูปแบบปิด: โดยสมมติว่าλ1>λ2{\displaystyle \lambda _{1}>\lambda _{2}}(โดยไม่เสียความเป็นทั่วไป) แล้ว ชม()=1+γ+ln(λ1λ2λ1λ2)+ψ(λ1λ1λ2),{\displaystyle {\begin{aligned}H(Z)&=1+\gamma +\ln \left({\frac {\lambda _{1}-\lambda _{2}}{\lambda _{1}\lambda _{2}}}\right)+\psi \left({\frac {\lambda _{1}}{\lambda _{1}-\lambda _{2}}}\right),\end{aligned}}} ที่ไหนγ{\displaystyle \gamma }คือค่าคงที่ออยเลอร์-มาสเชโรนีและψ(){\displaystyle \psi (\cdot )}คือฟังก์ชันไดแกมมา[ 7 ]

ในกรณีที่พารามิเตอร์อัตราเท่ากัน ผลลัพธ์ที่ได้คือการแจกแจงแบบ Erlangที่มีรูปร่าง 2 และพารามิเตอร์λ,{\displaystyle \lambda ,}ซึ่งในทางกลับกันก็เป็นกรณีพิเศษของ การ แจกแจงแกมมา

ผลรวมของตัวแปรสุ่มเอกซ์โพเนนเชียลอิสระ n ตัว Exp( λ) มี การแจกแจงแบบแกมมา(n, λ)

  • ถ้าX ~ Laplace(μ, β −1 )แล้ว | X − μ| ~ ประสบการณ์(β) [ 8 ]
  • ถ้าX ~ U (0, 1)แล้ว −log( X ) ~ Exp(1)
  • ถ้าX ~ Pareto(1, λ)แล้ว log( X ) ~ Exp(λ) [ 8 ]
  • ถ้าX ~ SkewLogistic(θ)แล้วบันทึก(1+อีX)~เอ็กซ์พี(θ){\displaystyle \log \left(1+e^{-X}\right)\sim \operatorname {Exp} (\theta )}.
  • ถ้าX ~ U (0, 1)แล้วลิมnnนาที(X1,,Xn)~เอ็กซ์พี(1){\displaystyle \lim _{n\to \infty }n\min \left(X_{1},\ldots ,X_{n}\right)\sim \operatorname {Exp} (1)}
  • การแจกแจงแบบเอกซ์โปเนนเชียลเป็นลิมิตของการแจกแจงแบบเบตา ที่ปรับขนาดแล้ว :ลิมnnเบต้า(1,n)=เอ็กซ์พี(1).{\displaystyle \lim _{n\to \infty }n\operatorname {Beta} (1,n)=\operatorname {Exp} (1).}
  • การแจกแจงแบบเอกซ์โปเน นเชียลเป็นกรณีพิเศษของการแจกแจงแบบเพียร์สัน ประเภทที่ 3
  • การแจกแจงแบบเอกซ์โพเนนเชียลเป็นกรณีพิเศษของการแจกแจงแบบแกมมาที่มีพารามิเตอร์รูปร่างเท่ากับ 1 [ 8 ]
  • อีxพี(1/2)=χ22{\displaystyle \mathrm {Exp} (1/2)=\chi _{2}^{2}}การแจกแจงไคกำลังสองที่มี 2 องศาอิสระ
  • ถ้าX ~ Exp(λ) และX ~ Exp(λ ) แล้ว:
    • เคX~เอ็กซ์พี(λเค){\displaystyle kX\sim \operatorname {Exp} \left({\frac {\lambda }{k}}\right)}การปิดภายใต้การปรับขนาดด้วยปัจจัยบวก
    • 1  + X ~ BenktanderWeibull (λ, 1) ซึ่งลดรูปเป็นการกระจายแบบเอกซ์โปเนนเชียลแบบตัดทอน 
    • ke X ~ พาเรโต ( k , แลม). [ 8 ]
    • e −λX ~ U (0, 1) .
    • e −X ~ Beta (λ, 1). [ 8 ]
    • 1 / k e X ~ PowerLaw( k , λ )
    • X~เรย์ลีย์(12λ){\displaystyle {\sqrt {X}}\sim \operatorname {Rayleigh} \left({\frac {1}{\sqrt {2\lambda }}}\right)}การกระจายแบบเรย์ลี[ 8 ]
    • X~ไวบูลล์(1λ,1){\displaystyle X\sim \operatorname {Weibull} \left({\frac {1}{\lambda }},1\right)}การกระจาย Weibull [ 8 ]
    • X2~ไวบูลล์(1λ2,12){\displaystyle X^{2}\sim \operatorname {Weibull} \left({\frac {1}{\lambda ^{2}}},{\frac {1}{2}}\right)}[ 8 ]
    • μ - β log(แลX ) ∼ Gumbel (μ, β)
    • X~เรขาคณิต(1อีλ){\displaystyle \lfloor X\rfloor \sim \operatorname {Geometric} \left(1-e^{-\lambda }\right)}การแจกแจงทางเรขาคณิตบน 0, 1, 2, 3,...
    • X~เรขาคณิต(1อีλ){\displaystyle \lceil X\rceil \sim \operatorname {Geometric} \left(1-e^{-\lambda }\right)}การแจกแจงทางเรขาคณิตบน 1,2,3,4,...
    • ถ้าเช่นนั้นด้วยวาย~{\displaystyle Y\sim }แกมมา ( n , แลมบ์ดา) = เออร์แลง ( n , แลมบ์ดา) ดังนั้นXวาย+1~พาเรโต(1,n){\displaystyle {\frac {X}{Y}}+1\sim \operatorname {Pareto} (1,n)}.
    • ถ้าเป็น แล ~ แกมมา ( k , θ) (รูปร่าง, การกำหนดพารามิเตอร์มาตราส่วน) ดังนั้นการกระจายตัวส่วนขอบของXคือLomax ( k , 1/θ) ซึ่ง เป็นส่วนผสมของแกมมา
    • แลสX − แลมY ~ ลาปลาซ(0, 1) .
    • นาที{ X , ..., X } ~ ประสบการณ์(แล + ... + แลม )
    • ถ้าX เป็นอิสระต่อกันและ λ = λ:
      • X1++Xn=ฉัน=1nXฉัน~{\displaystyle X_{1}+\cdots +X_{n}=\sum _{i=1}^{n}X_{i}\sim }Erlang ( n , λ) = Gamma ( n , λ) โดยมีพารามิเตอร์รูปร่างจำนวนเต็มnและพารามิเตอร์อัตรา λ [ 9 ]
      • ถ้าที=(X1++Xn)=ฉัน=1nXฉัน{\displaystyle T=(X_{1}+\cdots +X_{n})=\sum _{i=1}^{n}X_{i}}, แล้ว2λที~χ2n2{\displaystyle 2\lambda T\sim \chi _{2n}^{2}}.
      • X X ~ Laplace (0, λ −1 ).
      • XX+ฉัน=1nXฉัน~{\displaystyle {\frac {X}{X+\sum _{i=1}^{n}X_{i}}}\sim }เบต้า (1, n ). [ 10 ]โดยเฉพาะอย่างยิ่งXX+X1~{\displaystyle {\frac {X}{X+X_{1}}}\sim }U (0,1).
      • =λฉันXฉันλเจXเจ{\displaystyle Z={\frac {\lambda _{i}X_{i}}{\lambda _{j}X_{j}}}}มีฟังก์ชันความหนาแน่นความน่าจะเป็นเอฟ(z)=1(z+1)2{\displaystyle f_{Z}(z)={\frac {1}{(z+1)^{2}}}}สิ่งนี้สามารถนำไปใช้เพื่อหาช่วงความเชื่อมั่นสำหรับλฉันλเจ{\displaystyle {\frac {\lambda _{i}}{\lambda _{j}}}}.
    • ถ้า λ = 1:
      • μเบต้าบันทึก(อีX1อีX)~โลจิสติกส์(μ,เบต้า){\displaystyle \mu -\beta \log \left({\frac {e^{-X}}{1-e^{-X}}}\right)\sim \operatorname {Logistic} (\mu ,\beta )}การแจกแจงโลจิสติกส์
      • μเบต้าบันทึก(XฉันXเจ)~โลจิสติกส์(μ,เบต้า){\displaystyle \mu -\beta \log \left({\frac {X_{i}}{X_{j}}}\right)\sim \operatorname {Logistic} (\mu ,\beta )}
      • μ − σ บันทึก( X ) ~ GEV (μ, σ, 0)
      • นอกจากนี้ หากวาย~Γ(α,เบต้าα){\displaystyle Y\sim \Gamma \left(\alpha ,{\frac {\beta }{\alpha }}\right)}แล้วXวาย~เค(α,เบต้า){\displaystyle {\sqrt {XY}}\sim \operatorname {K} (\alpha ,\beta )}( การแจกแจง K )
  • ถ้าX~เอ็กซ์พี(1λ){\displaystyle X\sim \operatorname {Exp} \left({\frac {1}{\lambda }}\right)}และวายX{\displaystyle Y\mid X}~ Poisson( X )จากนั้นวาย~เรขาคณิต(11+λ){\displaystyle Y\sim \operatorname {Geometric} \left({\frac {1}{1+\lambda }}\right)} ( การแจกแจงทางเรขาคณิต )
  • การแจกแจงฮอยต์สามารถได้มาจากการแจกแจงเอกซ์โพเนนเชียลและการแจกแจงอาร์คไซน์
  • การแจกแจงแบบเอกซ์โปเนนเชียลเป็นลิมิตของการแจกแจงแบบเอกซ์โปเนนเชียลκในκ=0{\displaystyle \kappa =0}กรณี.
  • การแจกแจงแบบเอกซ์โพเนนเชียลเป็นลิมิตของการแจกแจงแกมมาทั่วไปแบบ κ ในα=1{\displaystyle \alpha =1}และν=1{\displaystyle \nu =1}กรณี:
    ลิม(α,ν)(0,1)พีκ(x)=(1+κν)(2κ)νΓ(12κ+ν2)Γ(12κν2)αλνΓ(ν)xαν1เอ็กซ์κ(λxα)=λอีλx{\displaystyle \lim _{(\alpha ,\nu )\to (0,1)}p_{\kappa }(x)=(1+\kappa \nu )(2\kappa )^{\nu }{\frac {\Gamma {\Big (}{\frac {1}{2\kappa }}+{\frac {\nu }{2}}{\Big )}}{\Gamma {\Big (}{\frac {1}{2\kappa }}-{\frac {\nu }{2}}{\Big )}}}{\frac {\alpha \lambda ^{\nu }}{\Gamma (\nu )}}x^{\alpha \nu -1}\exp _{\kappa }(-\lambda x^{\alpha })=\lambda e^{-\lambda x}}

การแจกแจงอื่นๆ ที่เกี่ยวข้อง:

การอนุมานทางสถิติ

ต่อไปนี้ สมมติว่าตัวแปรสุ่มXมีการแจกแจงแบบเอกซ์โปเนนเชียลโดยมีพารามิเตอร์อัตรา λ และx1,,xn{\displaystyle x_{1},\dotsc ,x_{n}}n ตัวอย่างอิสระจากXโดยมีค่าเฉลี่ยของแต่ละตัวอย่างx¯{\displaystyle {\bar {x}}}.

การประมาณค่าพารามิเตอร์

ตัวประมาณค่าความน่าจะเป็นสูงสุด

ตัว ประมาณค่า ความน่าจะเป็นสูงสุดสำหรับ λ ถูกสร้างขึ้นดังต่อไปนี้

ฟังก์ชันความน่าจะเป็นสำหรับ λ เมื่อกำหนดตัวอย่างx = ( x , ..., x ) ที่เป็นอิสระและมีการกระจายเหมือนกันซึ่งสุ่มมาจากตัวแปร มีดังนี้: แอล(λ)=ฉัน=1nλเอ็กซ์(λxฉัน)=λnเอ็กซ์(λฉัน=1nxฉัน)=λnเอ็กซ์(λnx¯),{\displaystyle L(\lambda )=\prod _{i=1}^{n}\lambda \exp(-\lambda x_{i})=\lambda ^{n}\exp \left(-\lambda \sum _{i=1}^{n}x_{i}\right)=\lambda ^{n}\exp \left(-\lambda n{\overline {x}}\right),}

ที่ไหน: x¯=1nฉัน=1nxฉัน{\displaystyle {\overline {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}} คือค่าเฉลี่ยของกลุ่มตัวอย่าง

อนุพันธ์ของลอการิทึมของฟังก์ชันความน่าจะเป็นคือ: λlnแอล(λ)=λ(nlnλλnx¯)=nλnx¯ {>0,0<λ<1x¯,=0,λ=1x¯,<0,λ>1x¯.{\displaystyle {\frac {d}{d\lambda }}\ln L(\lambda )={\frac {d}{d\lambda }}\left(n\ln \lambda -\lambda n{\overline {x}}\right)={\frac {n}{\lambda }}-n{\overline {x}}\ {\begin{cases}>0,&0<\lambda <{\frac {1}{\overline {x}}},\\[8pt]=0,&\lambda ={\frac {1}{\overline {x}}},\\[8pt]<0,&\lambda >{\frac {1}{\overline {x}}}.\end{cases}}}

ดังนั้น ค่าประมาณ ความน่าจะเป็นสูงสุดสำหรับพารามิเตอร์อัตราจึงเป็นดังนี้: λ^เอ็มแอลอี=1x¯=nฉัน=1nxฉัน{\displaystyle {\widehat {\lambda }}_{\text{mle}}={\frac {1}{\overline {x}}}={\frac {n}{\sum _{i=1}^{n}x_{i}}}}

นี่ไม่ใช่ตัวประมาณค่าที่เป็นกลางของλ,{\displaystyle \lambda ,}แม้ว่าx¯{\displaystyle {\overline {x}}}เป็น MLE ที่ไม่ลำเอียง[ 11 ] [ 12 ]ของ1/λ{\displaystyle 1/\lambda }และค่าเฉลี่ยของการกระจาย มีวิธีการต่างๆ เพื่อลดอคติของλ^เอ็มแอลอี{\displaystyle {\widehat {\lambda }}_{\text{mle}}}ส่งผลให้ได้ตัวประมาณค่าความน่าจะเป็นสูงสุดที่ได้รับการแก้ไขความลำเอียง

สามารถหาค่า ประมาณที่ทำให้ค่าความคลาดเคลื่อนกำลังสองเฉลี่ย (ดูเพิ่มเติม: การแลกเปลี่ยนระหว่างอคติและความแปรปรวน ) มีค่าน้อยที่สุดได้ โดยสมมติว่าขนาดตัวอย่างมากกว่าสอง พร้อมด้วยปัจจัยแก้ไขสำหรับค่าประมาณความน่าจะเป็นสูงสุด (MLE): λ^=(n2n)(1x¯)=n2ฉัน=1nxฉัน{\displaystyle {\widehat {\lambda }}=\left({\frac {n-2}{n}}\right)\left({\frac {1}{\bar {x}}}\right)={\frac {n-2}{\sum _{i=1}^{n}x_{i}}}} ค่านี้ได้มาจากค่าเฉลี่ยและความแปรปรวนของการแจกแจงแบบอินเวอร์สแกมมาอินวี-แกมมา(n,λ){\textstyle {\mbox{Inv-Gamma}}(n,\lambda )}[ 13 ]

UMVUE

ตัวประมาณค่าที่ไม่เอนเอียงที่มีความแปรปรวนต่ำสุดอย่างสม่ำเสมอสำหรับλ{\displaystyle \lambda }เป็น

λ^คุณวีคุณอี=n1ฉัน=1nxฉัน,{\displaystyle {\widehat {\lambda }}_{\mathrm {umvue} }={\frac {n-1}{\sum _{i=1}^{n}x_{i}}},}

ซึ่งแตกต่างจาก MLE ตรงที่โดยนิยามแล้วไม่มีอคติ

ข้อมูลของฟิชเชอร์

ข้อมูลฟิชเชอร์ซึ่งแสดงด้วยสัญลักษณ์ฉัน(λ){\displaystyle {\mathcal {I}}(\lambda )}สำหรับตัวประมาณค่าพารามิเตอร์อัตราλ{\displaystyle \lambda }กำหนดให้เป็นดังนี้: ฉัน(λ)=อี[(λบันทึกเอฟ(x;λ))2|λ]=(λบันทึกเอฟ(x;λ))2เอฟ(x;λ)x{\displaystyle {\mathcal {I}}(\lambda )=\operatorname {E} \left[\left.\left({\frac {\partial }{\partial \lambda }}\log f(x;\lambda )\right)^{2}\right|\lambda \right]=\int \left({\frac {\partial }{\partial \lambda }}\log f(x;\lambda )\right)^{2}f(x;\lambda )\,dx}

เมื่อแทนค่าการแจกแจงและแก้สมการจะได้ผลลัพธ์ดังนี้: ฉัน(λ)=0(λบันทึกλอีλx)2λอีλxx=0(1λx)2λอีλxx=λ2.{\displaystyle {\mathcal {I}}(\lambda )=\int _{0}^{\infty }\left({\frac {\partial }{\partial \lambda }}\log \lambda e^{-\lambda x}\right)^{2}\lambda e^{-\lambda x}\,dx=\int _{0}^{\infty }\left({\frac {1}{\lambda }}-x\right)^{2}\lambda e^{-\lambda x}\,dx=\lambda ^{-2}.}

สิ่งนี้จะกำหนดปริมาณข้อมูลที่แต่ละตัวอย่างอิสระของการแจกแจงแบบเอกซ์โปเนนเชียลมีเกี่ยวกับพารามิเตอร์อัตราที่ไม่ทราบค่าλ{\displaystyle \lambda }.

ช่วงความเชื่อมั่น

ช่วงความเชื่อมั่น 100(1 − α)% ที่แน่นอนสำหรับพารามิเตอร์อัตราของการแจกแจงแบบเอกซ์โพเนนเชียลกำหนดโดย: [ 14 ]2nλ^เอ็มแอลอีχα2,2n2<1λ<2nλ^เอ็มแอลอีχ1α2,2n2,{\displaystyle {\frac {2n}{{\widehat {\lambda }}_{\textrm {mle}}\chi _{{\frac {\alpha }{2}},2n}^{2}}}<{\frac {1}{\lambda }}<{\frac {2n}{{\widehat {\lambda }}_{\textrm {mle}}\chi _{1-{\frac {\alpha }{2}},2n}^{2}}}\,,} ซึ่งเท่ากับเช่นกัน 2nx¯χα2,2n2<1λ<2nx¯χ1α2,2n2,{\displaystyle {\frac {2n{\overline {x}}}{\chi _{{\frac {\alpha }{2}},2n}^{2}}}<{\frac {1}{\lambda }}<{\frac {2n{\overline {x}}}{\chi _{1-{\frac {\alpha }{2}},2n}^{2}}}\,,} โดยที่χ²p คือเปอร์เซ็นไทล์ที่100 ( p ) ของ การแจกแจงไคกำลังสองมีv องศาอิสระ n คือจำนวนการสังเกต และ x-bar คือค่าเฉลี่ยของตัวอย่าง สามารถประมาณค่าจุดปลายของช่วงความเชื่อมั่นได้อย่างง่าย ๆ โดยใช้การประมาณค่าแบบปกติของ การแจกแจง χ²p นี้ให้ค่าต่อไปนี้สำหรับช่วงความเชื่อมั่น 95% :λต่ำกว่า=λ^(11.96n)λด้านบน=λ^(1+1.96n){\displaystyle {\begin{aligned}\lambda _{\text{lower}}&={\widehat {\lambda }}\left(1-{\frac {1.96}{\sqrt {n}}}\right)\\\lambda _{\text{upper}}&={\widehat {\lambda }}\left(1+{\frac {1.96}{\sqrt {n}}}\right)\end{aligned}}}

การประมาณค่านี้อาจยอมรับได้สำหรับตัวอย่างที่มีองค์ประกอบอย่างน้อย 15 ถึง 20 รายการ[ 15 ]

การอนุมานแบบเบย์เซียนด้วยไพรเออร์แบบคอนจูเกต

การแจกแจงความน่าจะเป็นก่อนหน้าแบบสังยุคสำหรับการแจกแจงเอกซ์โพเนนเชียลคือการแจกแจงแกมมา (ซึ่งการแจกแจงเอกซ์โพเนนเชียลเป็นกรณีพิเศษของการแจกแจงแกมมา) การกำหนดพารามิเตอร์ของฟังก์ชันความหนาแน่นความน่าจะเป็นของแกมมาต่อไปนี้มีประโยชน์:

แกมมา(λ;α,เบต้า)=เบต้าαΓ(α)λα1เอ็กซ์(λเบต้า).{\displaystyle \operatorname {Gamma} (\lambda ;\alpha ,\beta )={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\lambda ^{\alpha -1}\exp(-\lambda \beta ).}

การแจกแจงความน่าจะเป็นภายหลังpสามารถแสดงได้ในรูปของฟังก์ชันความน่าจะเป็นที่กำหนดไว้ข้างต้นและค่าความน่าจะเป็นก่อนหน้าแบบแกมมา:

พี(λ)แอล(λ)Γ(λ;α,เบต้า)=λnเอ็กซ์(λnx¯)เบต้าαΓ(α)λα1เอ็กซ์(λเบต้า)λ(α+n)1เอ็กซ์(λ(เบต้า+nx¯)).{\displaystyle {\begin{aligned}p(\lambda )&\propto L(\lambda )\Gamma (\lambda ;\alpha ,\beta )\\&=\lambda ^{n}\exp \left(-\lambda n{\overline {x}}\right){\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\lambda ^{\alpha -1}\exp(-\lambda \beta )\\&\propto \lambda ^{(\alpha +n)-1}\exp(-\lambda \left(\beta +n{\overline {x}}\right)).\end{aligned}}}

ตอนนี้ความหนาแน่นของความน่าจะเป็นภายหลังpได้ถูกกำหนดไว้แล้ว ยกเว้นค่าคงที่การทำให้เป็นมาตรฐานที่ขาดหายไป เนื่องจากมีรูปแบบเป็นฟังก์ชันความหนาแน่นความน่าจะเป็นแบบแกมมา จึงสามารถเติมค่านี้ได้อย่างง่ายดาย และจะได้ผลลัพธ์ดังนี้:

พี(λ)=แกมมา(λ;α+n,เบต้า+nx¯).{\displaystyle p(\lambda )=\operatorname {Gamma} (\lambda ;\alpha +n,\beta +n{\overline {x}}).}

ในที่นี้พารามิเตอร์αสามารถตีความได้ว่าเป็นจำนวนการสังเกตการณ์ก่อนหน้า และβเป็นผลรวมของการสังเกตการณ์ก่อนหน้า ค่าเฉลี่ยภายหลังในที่นี้คือ: α+nเบต้า+nx¯.{\displaystyle {\frac {\alpha +n}{\beta +n{\overline {x}}}}.}

การอนุมานแบบเบย์เซียนด้วยไพรเออร์สำหรับการปรับเทียบ

การแจกแจงเอกซ์โพเนนเชียลเป็นหนึ่งในหลายๆ การแจกแจงทางสถิติที่มี โครงสร้าง แบบกลุ่มเนื่องจากโครงสร้างแบบกลุ่มนี้ การแจกแจงเอกซ์โพเนนเชียลจึงมีมาตรวัดฮาร์ ที่เกี่ยวข้อง ซึ่งก็คือ1/λ.{\displaystyle 1/\lambda .}การใช้มาตรวัด Haarเป็นค่าก่อนหน้า (เรียกว่าค่าก่อนหน้า Haar) ในการทำนายแบบ Bayesian จะให้ค่าความน่าจะเป็นที่ปรับเทียบได้อย่างสมบูรณ์แบบ สำหรับค่าพารามิเตอร์ที่แท้จริงพื้นฐานใดๆ[ 16 ] [ 17 ] [ 18 ]ค่าความน่าจะเป็นที่ปรับเทียบได้อย่างสมบูรณ์แบบมีคุณสมบัติที่ว่าค่าความน่าจะเป็นที่ทำนายไว้จะตรงกับความถี่ของเหตุการณ์นอกตัวอย่างอย่างแม่นยำ สำหรับเลขชี้กำลัง มีนิพจน์ที่แน่นอนสำหรับการทำนายแบบ Bayesian ที่สร้างขึ้นโดยใช้ค่าก่อนหน้า Haar ซึ่งกำหนดโดย

พีชมเอเอพีฉันโอ(xn+1x1,,xn)=nn+1(x¯)n(nx¯+xn+1)n+1.{\displaystyle p_{\rm {Haar-prior}}(x_{n+1}\mid x_{1},\ldots ,x_{n})={\frac {n^{n+1}\left({\overline {x}}\right)^{n}}{\left(n{\overline {x}}+x_{n+1}\right)^{n+1}}}.}

นี่เป็นตัวอย่างของการปรับเทียบการทำนายก่อนหน้า โดยเลือกค่าก่อนหน้าเพื่อปรับปรุงการปรับเทียบ (และในกรณีนี้ เพื่อให้การปรับเทียบสมบูรณ์แบบ) การปรับเทียบการทำนายก่อนหน้าสำหรับฟังก์ชันเอกซ์โปเนนเชียลโดยใช้ค่าก่อนหน้าของ Haar นั้นถูกนำไปใช้ใน แพ็กเกจซอฟต์แวร์ Rที่ชื่อ fitdistcp

การคาดการณ์เดียวกันนี้สามารถได้มาจากมุมมองอื่นๆ อีกหลายประการ ดังที่ได้กล่าวไว้ใน ส่วน การคาดการณ์ด้านล่าง

การเกิดขึ้นและการประยุกต์ใช้

การเกิดขึ้นของเหตุการณ์

การแจกแจงแบบเอกซ์โปเน นเชียลเกิดขึ้นตามธรรมชาติเมื่ออธิบายความยาวของช่วงเวลาระหว่างการมาถึงในกระบวนการปัวซง ที่เป็นเอกพันธ์

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

ในสถานการณ์จริง การสมมติว่าอัตราคงที่ (หรือความน่าจะเป็นต่อหน่วยเวลา) นั้นแทบจะไม่เป็นจริงเลย ตัวอย่างเช่น อัตราการโทรเข้าจะแตกต่างกันไปตามช่วงเวลาของวัน แต่ถ้าเราพิจารณาช่วงเวลาที่อัตราค่อนข้างคงที่ เช่น ตั้งแต่ 2 ถึง 4 โมงเย็นในวันทำงาน การแจกแจงแบบเอกซ์โปเนนเชียลสามารถใช้เป็นแบบจำลองโดยประมาณที่ดีสำหรับเวลาจนกว่าจะมีสายเรียกเข้าครั้งต่อไป ข้อควรระวังที่คล้ายกันนี้ใช้กับตัวอย่างต่อไปนี้ซึ่งให้ตัวแปรที่มีการแจกแจงแบบเอกซ์โปเนนเชียลโดยประมาณ:

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

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

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

ปรับการแจกแจงแบบเอกซ์โปเนนเชียลสะสมให้เข้ากับปริมาณน้ำฝนสูงสุดรายวันต่อปี

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

ในด้านอุทกวิทยาการกระจายแบบเอกซ์โพเนนเชียลใช้ในการวิเคราะห์ค่าสุดขั้วของตัวแปรต่างๆ เช่น ค่าสูงสุดรายเดือนและรายปีของปริมาณน้ำฝนรายวันและปริมาณการไหลของแม่น้ำ[ 19 ]

ภาพสีน้ำเงินแสดงตัวอย่างการปรับใช้การแจกแจงแบบเอกซ์โปเนนเชียลกับปริมาณน้ำฝนสูงสุดในหนึ่งวันต่อปีที่เรียงลำดับแล้ว โดยแสดงแถบความเชื่อมั่น 90% ตามการแจกแจงแบบทวินาม ด้วย ข้อมูลปริมาณน้ำฝนแสดงโดย การพล็ อตตำแหน่งเป็นส่วนหนึ่งของการวิเคราะห์ความถี่สะสม

ในการจัดการห้องผ่าตัด การกระจายระยะเวลาการผ่าตัดสำหรับกลุ่มการผ่าตัดที่ไม่มีเนื้อหางานเฉพาะเจาะจง (เช่น ในห้องฉุกเฉิน ซึ่งครอบคลุมการผ่าตัดทุกประเภท)

การทำนาย

เมื่อได้สังเกตตัวอย่าง ข้อมูล nจุดจาก1การแจกแจงเอกซ์โพเนนเชียลที่ไม่ทราบค่าแล้ว งานทั่วไปอย่างหนึ่งคือการใช้ตัวอย่างเหล่านี้เพื่อทำนายข้อมูลในอนาคตจากแหล่งเดียวกัน การแจกแจงทำนายทั่วไปสำหรับตัวอย่างในอนาคตคือการแจกแจงแบบปลั๊กอิน ซึ่งเกิดจากการแทนค่าประมาณที่เหมาะสมสำหรับพารามิเตอร์อัตราλลงในฟังก์ชันความหนาแน่นเอกซ์โพเนนเชียล ตัวเลือกทั่วไปของค่าประมาณคือค่าที่ได้จากหลักการความน่าจะเป็นสูงสุด และการใช้ค่านี้จะให้ความหนาแน่นของการทำนายสำหรับตัวอย่างในอนาคตx โดยมีเงื่อนไขจากตัวอย่างที่สังเกตได้x = ( x , ..., x ) ที่กำหนดโดย พีเอ็มแอล(xn+1x1,,xn)=(1x¯)เอ็กซ์(xn+1x¯).{\displaystyle p_{\rm {ML}}(x_{n+1}\mid x_{1},\ldots ,x_{n})=\left({\frac {1}{\overline {x}}}\right)\exp \left(-{\frac {x_{n+1}}{\overline {x}}}\right).}

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

การแจกแจงการทำนายที่ปราศจากปัญหาในการเลือกค่าความน่าจะเป็นล่วงหน้า (priors) ที่เกิดขึ้นภายใต้แนวทางแบบเบย์เซียนเชิงอัตวิสัย คือ

พีซีเอ็นเอ็มแอล(xn+1x1,,xn)=nn+1(x¯)n(nx¯+xn+1)n+1,{\displaystyle p_{\rm {CNML}}(x_{n+1}\mid x_{1},\ldots ,x_{n})={\frac {n^{n+1}\left({\overline {x}}\right)^{n}}{\left(n{\overline {x}}+x_{n+1}\right)^{n+1}}},}

ซึ่งสามารถพิจารณาได้ว่าเป็น

  1. การแจกแจงความเชื่อมั่นแบบความถี่ซึ่งได้มาจากการแจกแจงของปริมาณสำคัญxn+1/x¯{\displaystyle {x_{n+1}}/{\overline {x}}}; [ 20 ]
  2. ความน่าจะเป็นในการทำนายโปรไฟล์ที่ได้มาจากการกำจัดพารามิเตอร์λจากความน่าจะเป็นร่วมของx และλโดยการเพิ่มค่าสูงสุด[ 21 ]
  3. การแจกแจงความน่าจะเป็นภายหลังแบบเบย์เซียนเชิงวัตถุประสงค์ที่ได้มาโดยใช้ความน่าจะเป็นก่อนหน้าของเจฟฟรีย์ที่ ไม่ให้ข้อมูล 1/ λซึ่งเท่ากับความน่าจะเป็นก่อนหน้าของฮาร์ด้านขวาในกรณีนี้ การคาดการณ์ที่สร้างขึ้นโดยใช้ความน่าจะเป็นก่อนหน้าของฮาร์ด้านขวาจะรับประกันว่าจะให้ความน่าจะเป็นที่ปรับเทียบอย่างสมบูรณ์แบบ[ 22 ] [ 23 ]
  4. การแจกแจงการทำนายความน่าจะเป็นสูงสุดแบบปกติตามเงื่อนไข (CNML) จากการพิจารณาทางทฤษฎีสารสนเทศ[ 24 ]

ความแม่นยำของการแจกแจงการทำนายสามารถวัดได้โดยใช้ระยะทางหรือความแตกต่างระหว่างการแจกแจงเอกซ์โพเนนเชียลที่แท้จริงที่มีพารามิเตอร์อัตราλ และการแจกแจงการทำนายโดยอิงจากตัวอย่างx ความแตกต่างของ Kullback –Leiblerเป็นมาตรวัดความแตกต่างระหว่างการแจกแจงสองแบบที่ใช้กันทั่วไปและไม่ขึ้นอยู่กับพารามิเตอร์ ให้ Δ( λ || p ) แทนความแตกต่างของ Kullback–Leibler ระหว่างเอกซ์โพเนนเชียลที่มีพารามิเตอร์อัตราλ และการแจกแจงการทำนายpจะสามารถแสดงได้ว่า

อีλ0[Δ(λ0พีเอ็มแอล)]=ψ(n)+1n1บันทึก(n)อีλ0[Δ(λ0พีซีเอ็นเอ็มแอล)]=ψ(n)+1nบันทึก(n){\displaystyle {\begin{aligned}\operatorname {E} _{\lambda _{0}}\left[\Delta (\lambda _{0}\parallel p_{\rm {ML}})\right]&=\psi (n)+{\frac {1}{n-1}}-\log(n)\\\operatorname {E} _{\lambda _{0}}\left[\Delta (\lambda _{0}\parallel p_{\rm {CNML}})\right]&=\psi (n)+{\frac {1}{n}}-\log(n)\end{aligned}}}

โดยที่ค่าคาดหวังจะคำนวณจากค่าการแจกแจงแบบเอกซ์โพเนนเชียลที่มีพารามิเตอร์อัตราλ ∈ (0, ∞)และψ( · )คือฟังก์ชันไดแกมมา เห็นได้ชัดว่าการแจกแจงทำนาย CNML นั้นเหนือกว่าการแจกแจงแบบปลั๊กอินความน่าจะเป็นสูงสุดอย่างชัดเจนในแง่ของค่าเฉลี่ยความแตกต่างของ Kullback–Leibler สำหรับขนาดตัวอย่างทั้งหมดn > 0

การสร้างตัวแปรสุ่ม

วิธีการที่เข้าใจง่ายมากในการสร้างตัวแปร เลขชี้กำลังนั้น ใช้หลักการสุ่มแบบแปลงผกผัน : เมื่อกำหนดตัวแปรสุ่มUที่ได้จากการแจกแจงเอกรูปในช่วงหน่วย(0, 1)ตัวแปรนั้น...

ที=เอฟ1(ยู){\displaystyle T=F^{-1}(U)}

มีการกระจายแบบเอกซ์โปเนนเชียล โดยที่F −1คือฟังก์ชันควอนไทล์ซึ่งกำหนดโดย

เอฟ1(พี)=ln(1พี)λ.{\displaystyle F^{-1}(p)={\frac {-\ln(1-p)}{\lambda }}.}

นอกจากนี้ ถ้าUเป็นค่าเอกซ์โปเนนเชียลบนช่วง (0, 1) แล้ว 1 − U ก็เป็น ค่าเอกซ์โปเนนเชียลเช่นกัน ซึ่งหมายความว่าเราสามารถสร้างตัวแปรเลขชี้กำลังได้ดังนี้:

ที=ln(ยู)λ.{\displaystyle T={\frac {-\ln(U)}{\lambda }}.}

Knuth [ 25 ]และ Devroye [ 26 ]ได้กล่าวถึงวิธีการอื่นๆ ในการสร้างตัวแปรเลขชี้กำลัง

นอกจากนี้ยังมีวิธีที่รวดเร็วในการสร้างชุดตัวแปรเลขชี้กำลังที่เรียงลำดับไว้แล้วโดยไม่ต้องใช้ขั้นตอนการเรียงลำดับ[ 26 ]

ดูเพิ่มเติม

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

สรุปเนื้อหา

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

ข้อมูลสำคัญเกี่ยวกับ การแจกแจงแบบเอกซ์โปเนนเชียล

ในทฤษฎีความน่าจะเป็นและสถิติการแจกแจงเอกซ์โพเนนเชียลหรือการแจกแจงเอกซ์โพเนนเชียลเชิงลบคือการแจกแจงความน่าจะเป็นของระยะห่างระหว่างเหตุการณ์ในกระบวนการจุดปัวซง กล่าวคือ

ฟังก์ชันความหนาแน่นความน่าจะเป็น

ฟังก์ชัน ความหนาแน่นความน่าจะเป็น (pdf) ของการแจกแจงแบบเอกซ์โปเนนเชียลคือ

การกำหนดพารามิเตอร์ทางเลือก

การแจกแจงแบบเอกซ์โปเนนเชียลบางครั้งจะถูกกำหนดพารามิเตอร์โดยใช้ พารามิเตอร์มาตราส่วน β = 1/ λ ซึ่งก็คือค่าเฉลี่ยด้วยเช่นกัน: เอฟ ( x ; เบต้า ) = { 1 เบต้า อี − x / เบต้า x ≥ 0 , 0 x < 0. เอฟ ( x ; เบต้า ) = { 1 − อี − x / เบต้า x ≥ 0 , 0 x < 0.

ค่าเฉลี่ย ความแปรปรวน โมเมนต์ และมัธยฐาน

ค่าเฉลี่ยหรือ ค่าคาดหวัง ของตัวแปรสุ่ม X ที่มีการแจกแจงแบบเอกซ์โปเนนเชียล โดยมีพารามิเตอร์อัตรา λ จะกำหนดโดย อี ⁡ [ X ] = 1 λ . {\displaystyle \operatorname {E} [X]={\frac {1}{\lambda }}.}