İkinci derece durağanlık katı durağanlıktan daha zayıftır. İkinci dereceden durağanlık, birinci ve ikinci dereceden momentlerin (ortalama, varyans ve kovaryanslar) zaman içinde sabit olmasını ve dolayısıyla sürecin gözlemlendiği zamana bağlı olmamasını gerektirir. Özellikle, dediğin gibi, kovaryans sadece gecikme sırasına ( bağlıdır , ancak ölçülme zamanına bağlı değildir, C o v ( x t , x t - k ) = C o v ( x t + h , x t + h - k ) herkes içinkCov(xt,xt−k)=Cov(xt+h,xt+h−k) .t
Dediğiniz gibi sıkı bir durağanlık sürecinde, tüm siparişlerin anları, yani zaman boyunca ortak dağıtım sabiti kalır ortak dağılımı ile aynı olan X- t 1 + K + x t 2 + k + . . . + X T m + K tüm t 1 , t 2 , . . .Xt 1, Xt 2, . . . , Xt mXt 1 + k+ Xt 2 + k+ . . . + Xt m + k ve k .t 1 , t 2 , . . . , t mk
Bu nedenle, katı durağanlık ikinci dereceden durağanlığı içerir, ancak tersi doğru değildir.
Düzenle (@ whuber'ın yorumuna cevap olarak düzenlendi)
Önceki ifade, zayıf ve güçlü durağanlığın genel olarak anlaşılmasıdır. Zayıf anlamdaki durağanlığın daha güçlü bir anlamda durağanlık anlamına gelmediği düşüncesi sezgi ile anlaşabilse de, aşağıdaki yorumda whuber'ın işaret ettiği gibi ispat etmek o kadar kolay olmayabilir. Bu yorumda önerildiği gibi fikri göstermek faydalı olabilir.
İkinci derece durağan (zaman içinde ortalama, varyans ve kovaryans sabiti) olan bir süreci nasıl tanımlayabiliriz, ancak katı anlamda sabit değildir (daha yüksek dereceli anlar zamana bağlıdır)?
@Whuber tarafından önerildiği gibi (doğru anladıysam), farklı dağılımlardan gelen gözlem gruplarını birleştirebiliriz. Sadece bu dağılımların aynı ortalamaya ve varyansa sahip olmasına dikkat etmeliyiz (bu noktada birbirlerinden bağımsız olarak örneklendiklerini düşünelim). Bir yandan, biz örneğin Student gelen gözlemlerini üretebilir ile -Dağıtım 5 serbestlik derecesine. Ortalama sıfır ve varyans 5 / ( 5 - 2 ) = 5 / 3'tür . Başka taraftan, sıfır ortalama ve varyans ile Gauss dağılımı alabilir 5 / 3 .t55 / ( 5 - 2 ) = 5 / 35 / 3
Her iki dağılımlar aynı ortalamaya (sıfır) ve varyans (paylaşan ). Dolayısıyla, bu dağılımdan rastgele değerlerin birleştirilmesi en azından ikinci dereceden durağan olacaktır. Bununla birlikte, Gauss dağılımı ile yönetilen bu noktalardaki basıklık 3 , verinin Student t- dağılımından geldiği noktalarda 3 + 6 / ( 5 - 4 ) = 9 olacaktır . Bu nedenle, bu şekilde üretilen veriler tam anlamıyla durağan değildir, çünkü dördüncü dereceden momentler sabit değildir.5 / 33t3 + 6 / ( 5 - 4 ) = 9
Kovaryanslar da sabittir ve sıfıra eşittir, çünkü bağımsız gözlemleri düşündük. Bu önemsiz görünebilir, bu nedenle aşağıdaki otoregresif modele göre gözlemler arasında bir miktar bağımlılık yaratabiliriz.
ile
ε t ~ { N ( 0 , σ 2 = 5 / 3 )
yt= ϕ yt - 1+ ϵt,| ϕ | < 1,t = 1 , 2 , . . . , 120
εt∼ { N( 0 , σ2= 5 / 3 )t5Eğert ∈ [ 0 , 20 ] , [ 41 , 60 ] , [ 81 , 100 ]Eğert ∈ [ 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))
}
Sonuçlar beklediğim gibi değil:
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
t20