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