[英]Stationarity test using R
Using the stationarity test I am try to extract the result corresponding to the "p-value for T" line which in this case is 0, but I cannot seem to do it using R. 使用平稳性测试,我尝试提取与“ 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
where lg.day.ret can just be any real vector. 其中lg.day.ret可以是任何实向量。
Thanks in advance. 提前致谢。
Taking a guess at what you're looking for... 猜一猜您要寻找的东西...
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
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.