Một câu hỏi định lý giới hạn trung tâm khác


11

Đặt là một chuỗi các biến ngẫu nhiên Bernoulli độc lập với Đặt Cho thấy hội tụ trong phân phối cho biến thông thường tiêu chuẩn khi có xu hướng vô cùng.{Xn:n1}

P{Xk=1}=1P{Xk=0}=1k.
Sn=k=1n(Xk1k), Bn2=k=1nk1k2
SnBnZn

Nỗ lực của tôi là sử dụng Lyapunov CLT, do đó chúng tôi cần chứng minh có tồn tại một sao cho δ>0

limn1Bn2+δk=1nE[|Xk1k|2+δ]=0.

Vì vậy, đặt δ=1

k=1nE|Xkk1|3=k=1n(1k3k2+4k32k4)
Bn3=(k=1n1k1k2)(k=1n1k1k2)

Bằng cách đánh giá n lớn trên máy tính, nó cho thấy cả k=1nE|Xkk1|3Bn3 như n . Nhưng Bn3 tăng nhanh hơn Bn2 nên k=1nE|Xkk1|3Bn30 . Ai đó có thể giúp tôi chứng minh sự hội tụ này không?


7
Đây là ví dụ 27.3 về Xác suất và Đo lường của Patrick Billingsley.
Zhanxiong

Câu trả lời:


10

Có thể được hướng dẫn để chứng minh kết quả này từ các nguyên tắc đầu tiên và kết quả cơ bản , khai thác các thuộc tính của các hàm tạo tích lũy (chính xác như trong các bằng chứng tiêu chuẩn của Định lý giới hạn trung tâm). Nó đòi hỏi chúng ta phải hiểu tốc độ tăng trưởng của các số hài tổng quát cho Các tốc độ tăng trưởng này nổi tiếng và dễ dàng đạt được khi so sánh với các tích phân : chúng hội tụ cho và nếu không thì phân kỳ logarit cho . s = 1 , 2 , 210 . n 1 x - s d x s > 1 s = 1

H(n,s)=k=1nks
s=1,2,.1nxsdxs>1s=1

Đặt và . Theo định nghĩa, hàm tạo tích lũy (cgf) của là1 k n ( X k - 1 / k ) / B nn21kn(Xk1/k)/Bn

ψk,n(t)=logE(exp(Xk1/kBnt))=tkBn+log(1+1+exp(t/Bn)k).

Sự mở rộng chuỗi của phía bên tay phải, thu được từ việc mở rộng quanh , có dạngz = 0log(1+z)z=0

ψk,n(t)=(k1)2k2Bn2t2+k23k+26k3Bn3t3++kj1±(j1)!j!kjBnjtj+.

Các tử số của các phân số là đa thức tính bằng với số hạng đầu . Bởi vì mở rộng nhật ký hội tụ hoàn toàn cho , bản mở rộng này hội tụ hoàn toàn khik j - 1 | - 1 + điểm kinh nghiệm ( t / B n )kkj1|1+exp(t/Bn)k|<1

|exp(t/Bn)1|<k.

(Trong trường hợp nó hội tụ ở mọi nơi.) Đối với giá trị cố định và tăng của , độ phân kỳ (rõ ràng) của ngụ ý miền hội tụ tuyệt đối tăng lên tùy ý. Do đó, đối với bất kỳ cố định đủ lớn , sự mở rộng này hội tụ hoàn toàn.k n B n t nk=1knBntn

Do đó, đối với đủ lớn , chúng tôi có thể tính tổng theo kỳ hạn theo quyền hạn của để có được cgf của ,ψ k , n k t S n / B nnψk,nktSn/Bn

ψn(t)=k=1nψk,n(t)=12t2++1Bnj(k=1n(k1±(j1)!kj))tjj+.

Lấy các số hạng trong tổng số trên một lần yêu cầu chúng ta đánh giá các biểu thức tỷ lệ vớik

b(s,j)=1Bnjk=1nks

cho và . Sử dụng sự không triệu chứng của các số hài tổng quát được đề cập trong phần giới thiệu, nó dễ dàng theo saus = 1 , 2 , ... , jj3s=1,2,,j

Bn2=H(n,1)H(n,2)log(n)

cái đó

b(1,j)(log(n))1j/20

và (đối với )s>1

b(s,j)(log(n))j/20

khi phát triển lớn. Do đó, tất cả các điều khoản trong việc mở rộng ngoài hội tụ về 0, khi đó hội tụ đến cho bất kỳ giá trị nào của . Do sự hội tụ của cgf ngụ ý sự hội tụ của hàm đặc trưng, ​​chúng tôi kết luận từ Định lý liên tục Levy rằng tiếp cận một biến ngẫu nhiên có cgf là : đó là biến Bình thường tiêu chuẩn, QED .ψ n ( t ) t 2 ψ n ( t ) t 2 / 2 t S n / B n t 2 / 2nψn(t)t2ψn(t)t2/2tSn/Bnt2/2


Phân tích này chỉ ra mức độ hội tụ tinh tế như thế nào: trong khi trong nhiều phiên bản của Định lý giới hạn trung tâm, hệ số của là (đối với ), ở đây hệ số là chỉ : sự hội tụ chậm hơn nhiều. Theo nghĩa này, chuỗi các biến được tiêu chuẩn hóa "chỉ vừa đủ" trở thành Bình thường. O ( n 1 - j / 2 ) j 3 O ( ( ( log ( n ) ) 1 - j / 2 )tjO(n1j/2)j3O(((log(n))1j/2)

Chúng ta có thể thấy sự hội tụ chậm này trong một loạt mô phỏng. Các biểu đồ hiển thị lần lặp độc lập cho bốn giá trị của . Các đường cong màu đỏ là đồ thị của các hàm mật độ chuẩn thông thường để tham khảo trực quan. Mặc dù rõ ràng là có xu hướng dần dần về tính quy tắc, ngay cả ở (trong đó vẫn còn khá lớn) vẫn có sự phi quy tắc đáng chú ý, được chứng minh trong độ lệch (bằng trong mẫu này). (Không có gì ngạc nhiên khi độ lệch của biểu đồ này gần với , vì đó chính xác là thuật ngữ trong cgf là gì.) n n = 1000 ( log ( n ) ) - 1 / 20,38 0,35 (105nn=1000(log(n))1/20.380.35 t 3(log(n))1/2t3

Hình: biểu đồ cho n = 30, 100, 300, 1000

Đây là Rmã cho những người muốn thử nghiệm thêm.

set.seed(17)
par(mfrow=c(1,4))
n.iter <- 1e5
for(n in c(30, 100, 300, 1000)) {
  B.n <- sqrt(sum(rev((((1:n)-1) / (1:n)^2))))
  x <- matrix(rbinom(n*n.iter, 1, 1/(1:n)), nrow=n, byrow=FALSE)
  z <- colSums(x - 1/(1:n)) / B.n
  hist(z, main=paste("n =", n), freq=FALSE, ylim=c(0, 1/2))
  curve(dnorm(x), add=TRUE, col="Red", lwd=2)
}

6

Bạn đã có một câu trả lời tuyệt vời rồi. Nếu bạn cũng muốn hoàn thành bằng chứng của riêng mình, bạn có thể lập luận như sau:

Vì hội tụ cho tất cả và phân kỳ cho ( ở đây ), chúng tôi có thể viếtk=1n1/kii>1i=1

S(n):=k=1n(1k3k2+4k33k4)=k=1n1k+O(1).

Cùng một lập luận,

Bn2=k=1n1k+O(1).

Do đó, và, do đó,S(n)/Bn2=O(1)

S(n)/Bn3=O(1)(Bn2)1/20,

đó là những gì chúng tôi muốn thể hiện.


2

Đầu tiên các biến ngẫu nhiên của bạn không được phân phối giống hệt nhau nếu các phân phối phụ thuộc vào ;)k

Ngoài ra tôi sẽ không sử dụng ký hiệu của bạn như:Bn

  • chữ in hoa thường được dành cho các biến ngẫu nhiên.
  • nó chỉ là tổng của các phương sai nên tôi sẽ sử dụng một ký hiệu liên quan đến biểu tượng để làm cho điều này trở nên rõ ràng.σ

Sau đó, liên quan đến câu hỏi tôi không biết đây là bài tập hay nghiên cứu và công cụ nào bạn được phép sử dụng. Nếu bạn không cố gắng chứng minh lại các định lý đã biết, tôi chỉ nói rằng đó là một định lý giới hạn trung tâm cho RV không phân phối độc lập nhưng bị ràng buộc thống nhất và gọi nó là một ngày. Tôi không có nguồn tốt trong tay nhưng không quá khó để tìm một nguồn, ví dụ: hãy xem /mathpro/29508/is-there-a-central-limit-theorem- for-ràng buộc-không xác định-phân phối-ngẫu nhiên .

Chỉnh sửa: Điều xấu của tôi, tất nhiên điều kiện giới hạn đồng đều là không đủ, bạn cũng cần

k=1nσk2
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.