Một câu trả lời tốt, nghiêm ngặt, thanh lịch đã được đăng. Mục đích của việc này là để lấy được kết quả tương tự trong một cách mà có thể là một chút tiết lộ nhiều về cấu trúc cơ bản của . Nó cho thấy tại sao hàm mật độ xác suất (pdf) phải là số ít ở 0 .XY0
Nhiều thứ có thể được thực hiện bằng cách tập trung vào các hình thức phân phối thành phần :
là hai lần mộtbiến ngẫu nhiên U ( 0 , 1 ) . U ( 0 , 1 ) là một dạng tiêu chuẩn, "đẹp" đặc trưng của tất cả các phân phối đồng đều.XU(0,1)U(0,1)
là mười lần mộtbiến ngẫu nhiên U ( 0 , 1 ) .|Y|U(0,1)
Các dấu hiệu của tuân theo phân phối Rademacher: nó bằng - 1 hoặc 1 , mỗi với xác suất 1 / 2 .Y−111/2
(Bước cuối cùng này chuyển đổi một biến không âm thành phân phối đối xứng quanh , cả hai đuôi của chúng trông giống như phân phối ban đầu.)0
Do đó (a) đối xứng khoảng 0 và (b) giá trị tuyệt đối của nó là 2 × 10 = 20 lần tích của hai biến ngẫu nhiên U ( 0 , 1 ) độc lập .XY02×10=20U(0,1)
Sản phẩm thường được đơn giản hóa bằng cách lấy logarit. Thật vậy, ai cũng biết rằng nhật ký âm của biến có phân phối mũ (vì đây là cách đơn giản nhất để tạo ra các biến thiên hàm mũ ngẫu nhiên), từ đó nhật ký âm của sản phẩm của hai trong số đó có phân phối tổng của hai số mũ. Các Exponential là một Γ ( 1 , 1 ) phân phối. Các bản phân phối gamma có cùng tham số tỷ lệ rất dễ thêm: bạn chỉ cần thêm tham số hình dạng của chúng. Một Γ ( 1 , 1 ) cộng với một Γ ( 1U(0,1)Γ(1,1)Γ(1,1) do đó cóphân phối Γ ( 2 , 1 ) . hậu quả làΓ(1,1)Γ(2,1)
Biến ngẫu nhiên là phiên bản symmetrized của 20 lần so với hàm mũ của tiêu cực của một Γ ( 2 , 1 ) biến.XY20Γ(2,1)
Việc xây dựng PDF của từ đó của một U ( 0 , 1 ) phân phối được thể hiện từ trái sang phải, tiến hành trên đồng phục, đến mũ, với Γ ( 2 , 1 ) , với hàm mũ của tiêu cực của nó , với điều tương tự được nhân rộng bằng 20 , và cuối cùng là phiên bản đối xứng của điều đó. PDF của nó là vô hạn ở 0 , xác nhận sự gián đoạn ở đó.XYU(0,1)Γ(2,1)200
Chúng tôi có thể hài lòng để dừng lại ở đây. Chẳng hạn, việc mô tả đặc điểm này cho chúng ta một cách để tạo ra các nhận thức trực tiếp về , như trong biểu thức này :XYR
n <- 1; 20 * exp(-rgamma(n, 2, scale=1)) * ifelse(runif(n) < 1/2, -1, 1)
Phân tích này cũng cho thấy lý do tại sao pdf tăng lên . 0 Kỳ dị mà lần đầu tiên xuất hiện khi chúng ta coi là mũ của (các tiêu cực của) một phân phối, tương ứng với nhân một U ( 0 , 1 ) variate bằng nhau. Giá trị trong (nói) ε của 0 nảy sinh bằng nhiều cách, bao gồm (nhưng không giới hạn) khi (a) một trong những yếu tố nhỏ hơn ε hoặc (b) tất cả những yếu tố nhỏ hơn √Γ(2,1)U(0,1)ε0ε . Đó là căn bậc hai là vô cùng lớn hơnεbản thân khiεgần0. Điều này buộc rất nhiều xác suất, với số tiền lớn hơn √ε√εε0 , được vắt vào một khoảng thời gian dàiε. Để điều này có thể xảy ra, mật độ của sản phẩm phải trở nên lớn tùy ý ở mức0. Các thao tác tiếp theo - thay đổi kích thước theo hệ số20và đối xứng - rõ ràng sẽ không loại bỏ được điểm kỳ dị đó.ε√ε020
XYΓ(2,1)
f(t)dt=te−tdt, 0<t<∞.
t=−log(z)dt=−d(log(z))=−dz/z0<z<1tz
f(t)dt=−(−log(z)e−(−log(z))(−dz/z))=−log(z)dz, 0<z<1.
20
−log(z/20)d(z/20)=−120log(z/20)dz, 0<z<20.
z|z|−20202(−20,0)(0,20)
fXY(z)dzfXY(z)dz=−12120log(|z|/20), −20<z<20;=0 otherwise.