คำสั่งที่สองคือคำสั่งที่อ่อนแอกว่าคำสั่งที่เข้มงวด ลำดับที่สอง stationarity ต้องการช่วงเวลาลำดับที่หนึ่งและสอง (หมายถึงความแปรปรวนและความแปรปรวนร่วม) เป็นค่าคงที่ตลอดเวลาและดังนั้นจึงไม่ขึ้นอยู่กับเวลาที่กระบวนการสังเกต โดยเฉพาะอย่างยิ่งเมื่อคุณพูดความแปรปรวนร่วมจะขึ้นอยู่กับลำดับความล่าช้าเท่านั้น แต่ไม่ใช่ในเวลาที่วัดC o v ( x t , x t - k ) = C o v ( x t + h , x t + h - k )สำหรับทุกคนkCov(xt,xt−k)=Cov(xt+h,xt+h−k) .t
ในกระบวนการ stationarity เข้มงวดช่วงเวลาของการสั่งซื้อทั้งหมดที่คงที่ตลอดเวลากล่าวคือในขณะที่คุณพูดว่าการกระจายร่วมกันของเป็นเช่นเดียวกับการจัดจำหน่ายร่วมกันของเอ็กซ์ที1 + k + X เสื้อ2 + k + . . + X ตันเมตร+ kสำหรับทุกเสื้อ1 , T 2 , . .Xt1,Xt2,...,XtmXt 1 + k+ Xt 2 + k+ . . . + Xt m + kและkเสื้อ1 , T 2 , . . , เสื้อเมตรk
ดังนั้นความนิ่งแบบเข้มงวดเกี่ยวข้องกับความนิ่งลำดับที่สอง แต่การสนทนาไม่เป็นความจริง
แก้ไข (แก้ไขเป็นคำตอบสำหรับความคิดเห็นของ @ whuber)
ข้อความก่อนหน้าคือความเข้าใจโดยทั่วไปเกี่ยวกับความอ่อนแอและความแข็งแกร่งโดยทั่วไป แม้ว่าความคิดที่ว่าการคงอยู่ในความรู้สึกที่อ่อนแอไม่ได้หมายความถึงความคงที่ในความรู้สึกที่แข็งแกร่งอาจเห็นด้วยกับสัญชาตญาณ แต่ก็อาจไม่ตรงไปตรงมาพิสูจน์ มันจะมีประโยชน์ในการแสดงความคิดตามที่แนะนำในความคิดเห็นนั้น
เราจะกำหนดกระบวนการที่คงที่ลำดับที่สองได้อย่างไร (ค่าเฉลี่ยความแปรปรวนและความแปรปรวนร่วมคงที่ตลอดเวลา) แต่มันไม่คงที่ในความหมายที่เข้มงวด (ช่วงเวลาของลำดับที่สูงขึ้นกับเวลา)
ตามที่แนะนำโดย @whuber (ถ้าฉันเข้าใจถูกต้อง) เราสามารถเชื่อมโยงการสังเกตที่มาจากการแจกแจงที่ต่างกัน เราแค่ต้องระวังว่าการแจกแจงนั้นมีค่าเฉลี่ยและความแปรปรวนเท่ากัน (ณ จุดนี้ลองพิจารณาว่าพวกมันสุ่มตัวอย่างเป็นอิสระจากกัน) ในอีกด้านหนึ่งเราสามารถสร้างการสังเกตจากการแจกแจงแบบของนักเรียนด้วย5องศาอิสระ ค่าเฉลี่ยเป็นศูนย์และความแปรปรวนเป็น5 / ( 5 - 2 ) = 5 / 3 บนมืออื่นเราสามารถใช้การกระจายเสียนกับศูนย์ค่าเฉลี่ยและความแปรปรวน5 / 3เสื้อ55 / ( 5 - 2 ) = 5 / 35 / 3
ทั้งการกระจายแบ่งปันค่าเฉลี่ยเท่ากัน (ศูนย์) และความแปรปรวน ( ) ดังนั้นการต่อค่าสุ่มเข้าด้วยกันจากการแจกแจงเหล่านี้จะเป็นลำดับที่สองที่อยู่นิ่ง อย่างไรก็ตามโด่งที่จุดเหล่านั้นควบคุมโดยเสียนกระจายจะเป็น3ในขณะที่บรรดาจุดเวลาที่ข้อมูลมาจากนักศึกษาเสื้อ -distribution มันจะเป็น3 + 6 / ( 5 - 4 ) = 9 ดังนั้นข้อมูลที่สร้างขึ้นในลักษณะนี้จึงไม่คงที่อย่างเคร่งครัดเนื่องจากช่วงเวลาของคำสั่งที่สี่ไม่คงที่5 / 33เสื้อ3 + 6 / ( 5 - 4 ) = 9
ความแปรปรวนร่วมนั้นยังคงที่และมีค่าเท่ากับศูนย์เนื่องจากเราพิจารณาการสังเกตอย่างอิสระ สิ่งนี้อาจดูเล็กน้อยดังนั้นเราสามารถสร้างการพึ่งพาระหว่างการสังเกตตามแบบจำลองการตอบโต้อัตโนมัติต่อไปนี้
กับ
ε T ~ { N ( 0 , σ 2 = 5 / 3 )
Yเสื้อ= ϕ yt - 1+ ϵเสื้อ,| ϕ | < 1,T = 1 , 2 , . . , 120
εเสื้อ∼ { N( 0 , σ2= 5 / 3 )เสื้อ5ถ้าt ∈ [ 0 , 20 ] , [ 41 , 60 ] , [ 81 , 100 ]ถ้าt ∈ [ 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