r - 使用 R 进行平稳性检验

Question

使用平稳性检验，我尝试提取与“T 的 p 值”行相对应的结果，在这种情况下为 0，但我似乎无法使用 R 来做到这一点。

    > stationarity(lg.day.ret)
1
2
3
4
5
6
N = 2609, nblock = 11, n_block_max = 238, dt =     1.0000
7
8
9
10
11
12
13
14
15
16
17
18
Priestley-Subba Rao stationarity Test for lg.day.ret
----------------------------------------------------
Samples used              : 2609 
Samples available         : 2607 
Sampling interval         : 1 
SDF estimator             : Multitaper 
  Number of (sine) tapers : 5 
  Centered                : TRUE 
  Recentered              : FALSE 
Number of blocks          : 11 
Block size                : 237 
Number of blocks          : 11 
p-value for T             : 0 
p-value for I+R           : 0.08981712 
p-value for T+I+R         : 0 

其中 lg.day.ret 可以是任何真实向量。

提前致谢。

Answer 1

猜猜你在找什么……

library(fractal) # Don't forget to include where the function comes from
# in your example

results <- stationarity(ecgrr, n.block=8) #using sample data from
# fractal package

str(results) # Tells you more about the object
pvalues <- attr(results,"pvals")

##           T          I+R        T+I+R 
## 0.0001217759 0.0592977070 0.0013647583 

pvalues[1]
##           T 
## 0.0001217759