তবুও আরেকটি কেন্দ্রীয় সীমাবদ্ধ তাত্ত্বিক প্রশ্ন


11

যাক সঙ্গে স্বাধীন বের্নুলির র্যান্ডম ভেরিয়েবল একটা ক্রম হতে সেট দেখান যে distribution স্ট্যান্ডার্ড স্বাভাবিক ভেরিয়েবল তে বিতরণে রূপান্তরিত হয় কারণ অসীমের দিকে ঝোঁক।পি { এক্স = 1 } = 1 - পি { এক্স = 0 } = 1{Xn:n1}এসএন= এন কে=1(এক্সকে-1

P{Xk=1}=1P{Xk=0}=1k.
এসএন
Sn=k=1n(Xk1k), Bn2=k=1nk1k2
জেডএনSnBnZn

আমার প্রয়াস লায়াপুনভ সিএলটি ব্যবহার করার জন্য, সুতরাং আমাদের সেখানে একটি উপস্থিতি দেখাতে হবে যে, δ>0

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

সুতরাং সেট δ=1 \ যোগ_ {কে = 1} ^ {n} ই \ বাম | এক্স_ কে-কে ^ {- 1} \ ডান | ^ {3} = \ যোগ_ {কে = 1} ^ {n} \ বাম (rac frac {1} {k} - \ frac {3} {k ^ 2} + \ frac {4} {k ^ 3} - rac frac {2} {k ^ 4} \ ডান)

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

কম্পিউটারে বৃহত্তর এন এর জন্য মূল্যায়ন করে এটি দেখায় যে কীভাবে উভয় {যোগ_ {কে = 1} {}} n} ই | এক্স_কি-কে ^ {- 1} | ^ {3} \ থেকে \ ইনফটিk=1nE|Xkk1|3 এবং Bn3 যেমন n । তবে B_n ^ 3 B_n than 2 এরBn3 চেয়ে দ্রুত বৃদ্ধি পায় তাই \ frac {\ যোগ_ {কে = 1} {{n} ই | এক্স_কি-কে ^ {- 1} | ^ {3}} {বি_এন ^ 3} \ থেকে 0 । কেউ আমাকে এই রূপান্তর ধারণ করে প্রমাণ করতে সহায়তা করতে পারে?Bn2k=1nE|Xkk1|3Bn30


7
এটি প্যাট্রিক বিলিংসলে দ্বারা প্রবণতা এবং পরিমাপের 27.3 উদাহরণ ।
ঝাংসিওনগ

উত্তর:


10

প্রথম নীতি এবং মৌলিক ফলাফলগুলি থেকে এই ফলাফলটি প্রদর্শন করা নির্দেশক হতে পারে , ক্রমবর্ধমান উত্পাদনের কার্যাদি (কেন্দ্রীয় সীমাবদ্ধ তত্ত্বের আদর্শ প্রমাণ হিসাবে) ব্যবহার করে properties এটা তোলে আয় বৃদ্ধির হার বুঝতে আমাদের প্রয়োজন সাধারণ সুরেলা সংখ্যার জন্য এই বৃদ্ধির হারগুলি সুপরিচিত এবং সহজেই অখণ্ড তুলনা করে প্রাপ্ত হয় : তারা জন্য রূপান্তর করে এবং অন্যথায় জন্য ।

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

আসুন এবং । সংজ্ঞা অনুসারে, (এক্স_ এর কোমল্যান্ট জেনারেটিং ফাংশন (সিজিএফ) হয়n21kn(Xk1/k)/Bn

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

কাছাকাছি এর প্রসার থেকে প্রাপ্ত ডান হাতের সিরিজ সম্প্রসারণটি রূপটি গ্রহণ করেlog(1+z)z=0

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

ভগ্নাংশের numerators এবং মধ্যে polynomials হয় নেতৃস্থানীয় শব্দটি সঙ্গে । কারণ লগ সম্প্রসারণ সম্পূর্ণরূপে , যখন এই সম্প্রসারণটি একেবারে রূপান্তরিত হয়kkj1|1+exp(t/Bn)k|<1

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

( ক্ষেত্রে এটি সর্বত্র রূপান্তরিত হয়।) স্থির এবং ক্রমবর্ধমান মানের জন্য , এর (স্পষ্ট) বিচ্যুতিটি বোঝায় নিরঙ্কুশ রূপান্তরটির ডোমেনটি নির্বিচারে বড় হয়। সুতরাং, কোনও স্থির এবং পর্যাপ্ত পরিমাণে বড় এই সম্প্রসারণটি একেবারে রূপান্তরিত হয়।k=1knBntn

ভালোই বড় জন্য , তারপর, সুতরাং আমরা পৃথক যোগফল পারে উপর শব্দটি শব্দটি দ্বারা ক্ষমতা মধ্যে এর cgf প্রাপ্ত ,nψk,nktSn/Bn

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

একবারে বেশি পরিমাণে শর্তাদি গ্রহণের জন্য আমাদের আনুপাতিক এক্সপ্রেশনগুলির মূল্যায়ন করতে হবেk

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

জন্য এবং । পরিচিতিতে উল্লিখিত জেনারেলাইজড হারমোনিক সংখ্যার অ্যাসিম্পটোটিকগুলি ব্যবহার করে, এটি সহজেই অনুসরণ করেj3s=1,2,,j

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

যে

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

এবং ( )s>1

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

হিসাবে বড় হয়। ফলে সম্প্রসারণ সমস্ত পদ পরলোক শুন্যতে মিলিত, কোথা এগোয় করার কোন মানের জন্য । যেহেতু সিজিএফের রূপান্তরটি বৈশিষ্ট্যযুক্ত ফাংশনটির রূপান্তরকে বোঝায়, আমরা লেভি কন্টিনিউটি উপপাদ্য থেকে উপসংহারে যে একটি র্যান্ডম ভেরিয়েবলের কাছে পৌঁছায় যার সিজিএফ 2/2 : এটি স্ট্যান্ডার্ড নরমাল ভেরিয়েবল, কিউইডিnψn(t)t2ψn(t)t2/2tSn/Bnt2/2


এই বিশ্লেষণের খুঁজে পান ঠিক কিভাবে উপাদেয় অভিসৃতি হল: যেহেতু সহগ কেন্দ্রীয় সীমা উপপাদ্য অনেক সংস্করণ মধ্যে হয় (জন্য ), এখানে সহগ হয় কেবল : রূপান্তরটি অনেক ধীর গতি সম্পন্ন। এই অর্থে মানকযুক্ত ভেরিয়েবলগুলির ক্রম "সবেমাত্র সবে" স্বাভাবিক হয়ে যায়।tjO(n1j/2)j3O(((log(n))1j/2)

আমরা এই ধীর অভিযাত্রাটি বিভিন্ন সিমুলেশনগুলিতে দেখতে পাই। হিস্টোগ্রামগুলি চারটি মানের জন্য স্বতন্ত্র পুনরাবৃত্তি প্রদর্শন করে । লাল রেখাচিত্রগুলি ভিজ্যুয়াল রেফারেন্সের জন্য স্ট্যান্ডার্ড সাধারণ ঘনত্ব ফাংশনের গ্রাফ। যদিও স্পষ্টতই স্বাভাবিকতার দিকে ধীরে ধীরে প্রবণতা রয়েছে, এমনকি (যেখানে এখনও বড় আকারের) সেখানে প্রশংসনীয় অ-স্বাভাবিকতা রয়ে গেছে, যেমন স্কিউনেসে প্রমাণিত ( এই নমুনায় সমান )। (এই হিস্টগ্রামের সঙ্কোচনের to এর খুব কাছাকাছি হওয়া কোনও আশ্চর্যের বিষয় নয় , কারণ এটি সিজিএফ-তে term টার্মটি হ'ল তাই)105nn=1000(log(n))1/20.380.35(log(n))1/2t3

চিত্র: এন = 30, 100, 300, 1000 এর জন্য হিস্টোগ্রাম

Rযারা আরও পরীক্ষা-নিরীক্ষা করতে চান তাদের জন্য এখানে কোড।

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

আপনার ইতিমধ্যে একটি দুর্দান্ত উত্তর আছে। আপনি যদি নিজের প্রমাণটিও সম্পূর্ণ করতে চান তবে আপনি নীচের মতো তর্ক করতে পারেন:

যেহেতু সমস্ত জন্য রূপান্তর এবং ( এখানে ) এর জন্য বিভক্ত হয় , তাই আমরা লিখতে পারিk=1n1/kii>1i=1

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

একই যুক্তি দ্বারা,

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

ফলস্বরূপ, এবং এইভাবে,S(n)/Bn2=O(1)

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

যা আমরা দেখাতে চেয়েছিলাম


2

ডিস্ট্রিবিউশনগুলি ;) এর উপর নির্ভর করে যদি প্রথমে আপনার এলোমেলো ভেরিয়েবলগুলি অভিন্নভাবে বিতরণ করা হয় না )k

এছাড়াও আমি আপনার সংকেতটি ব্যবহার করব না :Bn

  • মূলধন বর্ণগুলি সাধারণত এলোমেলো ভেরিয়েবলের জন্য সংরক্ষিত থাকে।
  • এটি কেবল রূপগুলির সমষ্টি তাই আমি এই স্পষ্ট করে তুলতে একটি প্রতীক জড়িত একটি স্বরলিপি ব্যবহার করব ।σ

তারপরে প্রশ্নটি সম্পর্কে আমি জানি না এটি কোনও অনুশীলন বা গবেষণা এবং কোন সরঞ্জামগুলি আপনাকে ব্যবহার করার অনুমতি দিয়েছে তা আমি জানি না। যদি আপনি পরিচিত উপপাদাগুলি পুনরায় প্রমাণ করার চেষ্টা না করে থাকেন, তবে আমি কেবলমাত্র এটিই বলব যে এটি স্বতন্ত্রভাবে অ-পরিচয়যুক্ত বিতরণ করা হলেও সমানভাবে বেঁধে দেওয়া আরভির জন্য একটি কেন্দ্রীয় সীমাবদ্ধ উপপাদ্য এবং এটি একটি দিন কল করে। আমার হাতে কোনও ভাল উত্স নেই তবে এটি খুঁজে পাওয়া খুব বেশি কঠিন হওয়া উচিত নয়, উদাহরণস্বরূপ /mathpro/29508/is-there-a-central-limit-theorem- এ দেখুন সীমাবদ্ধ-অ-অভিন্ন-বিতরণ-এলোমেলো জন্য

সম্পাদনা করুন: আমার খারাপ, অবশ্যই অভিন্ন শর্তের অবস্থা যথেষ্ট নয়,

k=1nσk2
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.