Тест стационарности с использованием R

Question

Используя тест на стационарность, я пытаюсь извлечь результат, соответствующий строке "p-value for T", которая в данном случае равна 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