if ( gqtest(paired.ds$unemployment ~1)$p.value < 0.1) {

print ("**********P-VALUE VALUE LESS THAN 0.05 WE CAN SAY THAT WE CAN REJECT NON CAUSALITY, WE MAY HAVE CAUSALITY *************")

#### predicting next 5 months, *****  here is my doubt! How is the predict calculated?
##### how to know two time series in better than one?
p<- predict(var, n.ahead=5, ci=0.95)
fcst = forecast(var, h = 5)
print( accuracy( fcst, testdata))

print ("**********P-VALUE VALUE LESS THAN 0.05 WE CAN SAY THAT WE CAN REJECT NON CAUSALITY, WE MAY HAVE CAUSALITY *************")


#### predicting next 5 months, *****  here is my doubt! How is the predict calculated?
##### how to know two time series in better than one?
p<- predict(var, n.ahead=5, ci=0.95)
fcst = forecast(var, h = 5)
print( accuracy(fcst, testdata))
plot(var)
plot(p)