data(syg_data)
un <- c('ca_u','jp_u','us_u')
ip <- c('ca_ip','jp_ip','us_ip')
cnames <- c('Canada','Japan','United States')
case <- c(5,5,5)
p_lags <- c(3,4,5)
q_lags <- c(3,4,5)
########################################
out_nardl <- lapply(1:length(un), function(i)
nardl_uecm(x = syg_data,
decomp = un[[i]],
control =NULL,
c_q_order = c(2),
p_order =p_lags[[i]],
q_order =q_lags[[i]],
dep_var = ip[[i]],
graph_save = FALSE,
case = 3))
names(out_nardl) <- cnames
out_nardl$Canada$UECM
output <- output_ren(x=out_nardl,listn = un,
D.patern="D.ca_ip|D.jp_ip|D.us_ip",D.repl="D.output",
l.patern ="ca_ip|jp_ip|us_ip",l.repl="output")
output <- output_ren(x=output,listn = un,
D.patern="D.ca_u|D.jp_u|D.us_u",D.repl="D.unemployment",
l.patern ="ca_u|jp_u|us_u",l.repl="unemployment")
output$Canada$UECM
Run the code above in your browser using DataLab