I want to test multiple possible cointegration relationships with the Johansen cointegration test. I'm currently using the urca package in R with the ca.jo test. I was going to use the Bonferroni correction to correct the familywise error rate. However, ca.jo doesn't report p-values. It only reports test statistics and critical values:
######################
# Johansen-Procedure #
######################
Test type: maximal eigenvalue statistic (lambda max) , with linear trend
Eigenvalues (lambda):
[1] 0.335639191 0.001256000
Values of teststatistic and critical values of test:
test 10pct 5pct 1pct
r <= 1 | 1.26 6.50 8.18 11.65
r = 0 | 408.52 12.91 14.90 19.19
[More ca.jo output snipped]
How can I correct for rejecting the null incorrectly when I perform multiple cointegration tests?