data(data1)
# Replicating Table 5 (excess return)
y=data1$ret.nyse.vw*100 -data1$tbill*100
x=cbind(log(data1$dy.nyse), data1$tbill*100); k=ncol(x)
p=PR.order(x,y,pmax=10)$p.bic; # AR(1)
Run the code above in your browser using DataCamp Workspace