দ্বিতীয় আদেশের স্থিরতা কড়া স্টেরারিটির চেয়ে দুর্বল। দ্বিতীয় আদেশের স্থিতিস্থানের প্রয়োজন যে প্রথম এবং দ্বিতীয় ক্রমের মুহুর্তগুলি (অর্থাত্, বিবর্তন এবং সমবায়িকাগুলি) পুরো সময় স্থির থাকে এবং তাই প্রক্রিয়াটি যে সময়ের মধ্যে পর্যবেক্ষণ করা হয় তার উপর নির্ভর করে না। বিশেষত, যেমনটি আপনি বলেছেন, সমবায়ুতা কেবলমাত্র পিছনের ক্রমের উপর নির্ভর করে, , তবে এটি যে সময় পরিমাপ করা হয় তা নয়, সি ও ভি ( x টি , এক্স টি - কে ) = সি ও ভি ( এক্স টি + এইচ) , x t + h - k ) সবার জন্যkCov(xt,xt−k)=Cov(xt+h,xt+h−k) ।t
একটি কঠোর স্টেশনারিটি প্রক্রিয়াতে, সমস্ত আদেশের মুহুর্তগুলি পুরো সময় জুড়ে স্থির থাকে, যেমন আপনি বলেছিলেন, এর যৌথ বন্টন । । । , এক্স টি এম এক্স টি 1 + কে + এক্স টি 2 + কে + এর যৌথ বিতরণ হিসাবে সমান । । । + + এক্স টি মিটার + + ট সবার জন্য টি 1 , T 2 , । । ।Xt1,Xt2,...,Xtmএক্সt 1 + কে+ এক্সটি 2 + কে+ । । । + এক্সt মি + কে এবং কে ।T 1 , T 2 , । । । , t মিট
অতএব, কঠোরভাবে স্থিতিশীলতার সাথে দ্বিতীয় ক্রম স্থিরত্ব জড়িত তবে কনভার্সটি সত্য নয়।
সম্পাদনা করুন (@ whuber এর মন্তব্যের জবাব হিসাবে সম্পাদিত)
পূর্ববর্তী বিবৃতিটি দুর্বল এবং দৃ strong় স্থিরতার সাধারণ উপলব্ধি। যদিও দৃ weak় অর্থে স্থিরতা বোঝায় না এই ধারণাটি দৃu় অর্থে স্বজ্ঞাততার সাথে একমত হতে পারে তবে নীচের মন্তব্যে হুবহু দ্বারা চিহ্নিত করা যেমন প্রমাণের পক্ষে এতটা সোজা নাও হতে পারে। মন্তব্যটিতে পরামর্শ মতো ধারণাটি চিত্রিত করতে এটি সহায়ক হতে পারে।
দ্বিতীয় অর্ডারের স্থিতিশীল এমন কোনও প্রক্রিয়াটিকে আমরা কীভাবে সংজ্ঞায়িত করতে পারি (মানে, পুরো সময় জুড়ে বৈকল্পিকতা এবং কোভেরিয়েন্স ধ্রুবক) তবে এটি কঠোর অর্থে স্থির নয় (উচ্চতর আদেশের মুহুর্তগুলি সময়ের উপর নির্ভর করে)?
@ হুবারের পরামর্শ অনুসারে (যদি আমি সঠিকভাবে বুঝতে পারি) তবে আমরা বিভিন্ন বিতরণ থেকে আসা পর্যবেক্ষণগুলির ব্যাচগুলিকে একত্রিত করতে পারি। আমাদের কেবল সাবধান হওয়া দরকার যে সেই বিতরণগুলির একই গড় এবং বৈচিত্র রয়েছে (এই মুহুর্তে আসুন বিবেচনা করুন যে তারা একে অপরের থেকে স্বাধীনভাবে নমুনা নিচ্ছেন)। একদিকে আমরা উদাহরণস্বরূপ 5 ডিগ্রি স্বাধীনতার সাথে শিক্ষার্থীর বিতরণ থেকে পর্যবেক্ষণগুলি তৈরি করতে পারি । গড় শূন্য হয় এবং ভ্যারিয়েন্স হয় 5 / ( 5 - 2 ) = 5 / 3 । অন্য দিকে, আমরা শূন্য গড় এবং ভ্যারিয়েন্স সঙ্গে গসিয়ান বন্টন নিতে পারেন 5 / 3 ।টি55 / ( 5 - 2 ) = 5 / 35 / 3
উভয় ডিস্ট্রিবিউশন একই গড় (শূন্য) এবং ভ্যারিয়েন্স (ভাগ )। সুতরাং, এই বিতরণ থেকে এলোমেলো মানগুলির উপসংহারটি কমপক্ষে, দ্বিতীয়-ক্রমের স্থিতিশীল হবে। তবে গাউসীয় বিতরণ দ্বারা পরিচালিত সেই পয়েন্টগুলিতে কুর্তোসিসটি 3 হবে , যখন সেই সময় পয়েন্টগুলিতে যেখানে শিক্ষার্থীর টি- বিতরণ থেকে ডেটা আসে এটি 3 + 6 / ( 5 - 4 ) = 9 হবে । সুতরাং, এই উপায়ে উত্পন্ন ডেটা কঠোর অর্থে স্থিতিশীল নয় কারণ চতুর্থ ক্রমের মুহূর্তগুলি স্থির হয় না।5 / 33টি3 + 6 / ( 5 - 4 ) = 9
সমবায়ীরাও স্থির এবং শূন্যের সমান, যেহেতু আমরা স্বাধীন পর্যবেক্ষণ বিবেচনা করি। এটি তুচ্ছ মনে হতে পারে, তাই আমরা নিম্নলিখিত স্বতঃসংশ্লিষ্ট মডেল অনুসারে পর্যবেক্ষণের মধ্যে কিছুটা নির্ভরতা তৈরি করতে পারি।
সঙ্গে
ε টি ~ { এন ( 0 , σ 2 = 5 / 3 )
Yটি= ϕ yটি - 1+ + εটি,| ϕ | < 1,T = 1 , 2 , । । । , 120
εটি∼ { এন( 0 , σ)2= 5 / 3 )টি5যদিটি ∈ [ 0 , 20 ] , [ 41 , 60 ] , [ 81 , 100 ]যদিটি ∈ [ 21 , 40 ] , [ 61 , 80 ] , [ 101 , 120 ]।
| ϕ | < 1
20ϕ = 0.8n = 240
# this function is required below
kurtosis <- function(x)
{
n <- length(x)
m1 <- sum(x)/n
m2 <- sum((x - m1)^2)/n
m3 <- sum((x - m1)^3)/n
m4 <- sum((x - m1)^4)/n
b1 <- (m3/m2^(3/2))^2
(m4/m2^2)
}
# begin simulation
set.seed(123)
n <- 240
Mmeans <- Mvars <- Mcovs <- Mkurts <- matrix(nrow = 1000, ncol = n/20)
for (i in seq(nrow(Mmeans)))
{
eps1 <- rnorm(n = n/2, sd = sqrt(5/3))
eps2 <- rt(n = n/2, df = 5)
eps <- c(eps1[1:20], eps2[1:20], eps1[21:40], eps2[21:40], eps1[41:60], eps2[41:60],
eps1[61:80], eps2[61:80], eps1[81:100], eps2[81:100], eps1[101:120], eps2[101:120])
y <- arima.sim(n = n, model = list(order = c(1,0,0), ar = 0.8), innov = eps)
ly <- split(y, gl(n/20, 20))
Mmeans[i,] <- unlist(lapply(ly, mean))
Mvars[i,] <- unlist(lapply(ly, var))
Mcovs[i,] <- unlist(lapply(ly, function(x)
acf(x, lag.max = 1, type = "cov", plot = FALSE)$acf[2,,1]))
Mkurts[i,] <- unlist(lapply(ly, kurtosis))
}
ফলাফলগুলি আমার প্রত্যাশা মতো নয়:
round(colMeans(Mmeans), 4)
# [1] 0.0549 -0.0102 -0.0077 -0.0624 -0.0355 -0.0120 0.0191 0.0094 -0.0384
# [10] 0.0390 -0.0056 -0.0236
round(colMeans(Mvars), 4)
# [1] 3.0430 3.0769 3.1963 3.1102 3.1551 3.2853 3.1344 3.2351 3.2053 3.1714
# [11] 3.1115 3.2148
round(colMeans(Mcovs), 4)
# [1] 1.8417 1.8675 1.9571 1.8940 1.9175 2.0123 1.8905 1.9863 1.9653 1.9313
# [11] 1.8820 1.9491
round(colMeans(Mkurts), 4)
# [1] 2.4603 2.5800 2.4576 2.5927 2.5048 2.6269 2.5251 2.5340 2.4762 2.5731
# [11] 2.5001 2.6279
টি20