[R] extract p-value from urppTest

Pfaff, Bernhard Dr. Bernhard_Pfaff at fra.invesco.com
Mon May 8 10:20:19 CEST 2006

Hello Anne,

the function 'ur.pp' contained in the package 'urca' has been ported
into the package 'fSeries' (see the documentaion of ?urppTest). 
help("ur.pp-class", package="urca") will tell you that p-values are not
part of this class, i.e., these are not computed but critical values are
Now, inside of the function 'urppTest', 'ur.pp' is called and an object
of class 'fHTEST' is returned (you can see this, by simply typing
'urppTest'). The class 'fHTEST' does contain the slot 'test' (it is a
list). You can see this by spotting at getClass("fHTEST"). However, what
is stored into test is the outcome of 'show(your ur.pp-object)', i.e.
the print-out of the test statistic.
In contrast, by typing 'adfTest' you will easily see, that the p-values
are calculated within this function and are part of the slot 'test'.


>Dear List,
>How do I pick the p-value out of the urppTest result?
>For adfTest the p-value can be extracted by
>  A2 at test$p.value
>  A2 <- adfTest(myData[,i], lags=2, type=c("c"))
>What do I do for urppTest? The above doesn't seem to work. 
>There is a slot @test with
>$output, which is a list of various test results that didn't 
>want to give away only the
>p-value (I'm fairly new to R, so this might well be my fault). 
>I thought I had found a way around this by
> # dissolve testresult into lines
>	listasvectors <- unlist(P1 at test$output)
> # pick the line containing p-value
>      getpvalue <- unlist(strsplit(listasvectors[17], " "))
> # isolate the p-value (number only)
>	getpvalue[14]
>However, I'm doing this in a loop and it seems for each test 
>result there is a different
>number of entries or seperators " ". For example, "<" 
>sometimes ends up being displayed as
>the p-value.
>I'd be very happy about any help on how to extract the p-value.
>R-help at stat.math.ethz.ch mailing list
>PLEASE do read the posting guide! 
Confidentiality Note: The information contained in this mess...{{dropped}}

More information about the R-help mailing list