Making a nice table out of similar object elements in R?

Question

Not sure this has been asked before exactly but I have a bunch of coefficient tests for many different models:

coeftest(arima(detrend, order = c(0,0,2))) #MA2
coeftest(arima(detrend, order = c(0,0,3))) #MA3
coeftest(arima(detrend, order = c(0,0,4))) #MA4
coeftest(arima(detrend, order = c(1,0,0))) #AR1
coeftest(arima(detrend, order = c(2,0,0))) #AR2
coeftest(arima(detrend, order = c(1,0,1))) #ARMA(1,1)
coeftest(arima(detrend, order = c(2,0,1))) #ARMA(2,1) 
coeftest(arima(detrend, order = c(2,0,1))) #ARMA(1,2)

How would I make a nice table where I can have the model/object name if I store those tests in objects and the resulting time series estimates with their p-values? Something like:

Model  AR Terms MA Terms P-value
MA1      NA       1       0.006
AR2      1        NA      0.0042
         2        NA      0.0020
ARMA11   1        0       0.0031
         0        1       0.0005

Any arima model will do for a sample, so I will provide this ar2 process:


phi.1 <- 1.2
phi.2 <- -0.6
sigma.e <- 9
Y <- rnorm(2, 0, sigma.e)
n <- 200

for(t in 3:(n+2)){
    new <- phi.1*(Y[t-1]) + phi.2*(Y[t-2]) + rnorm(1, 0, sigma.e)
    Y <- c(Y, new)
}

ar2 <- as.ts(Y)

please provide a sample of your data. In order to help you we have to be able to reproduce the lines of code above. — Omry Atia
Okay, I added some model generation code. Any time series model will really do for this. — Michael

Humpelstielzchen Humpelstielzchen · Accepted Answer · 2019-05-14T05:19:37

Can you clarify your desired output a little bit more? The models have several AR terms and several MA terms resulting in many P-values.

UPDATE: So like this?

library(tidyverse)

test <- list(
  "MA2" = coeftest(arima(ar2 , order = c(0,0,2))),
  "MA3" = coeftest(arima(ar2 , order = c(0,0,3))), 
  "MA4" = coeftest(arima(ar2 , order = c(0,0,4))), 
  "AR1" = coeftest(arima(ar2 , order = c(1,0,0))), 
  "AR2" = coeftest(arima(ar2 , order = c(2,0,0))), 
  "ARMA(1,1)" = coeftest(arima(ar2 , order = c(1,0,1))),
  "ARMA(2,1)" = coeftest(arima(ar2 , order = c(2,0,1))), 
  "ARMA(1,2)" = coeftest(arima(ar2 , order = c(2,0,1))))

map_df(test, tidy, .id = "model") %>%
  select(-std.error, -statistic) %>%
  complete(model, term, fill = list(NA)) %>%
  nest(estimate, p.value, .key = 'stats') %>%
  spread(key = term, value = stats) %>%
  unnest(.sep = '_') -> res


# A tibble: 8 x 15
  model ar1_estimate ar1_p.value ar2_estimate ar2_p.value intercept_estim~
  <chr>        <dbl>       <dbl>        <dbl>       <dbl>            <dbl>
1 AR1          0.750   8.92e- 60       NA       NA                   1.04 
2 AR2          1.29    1.54e-152       -0.714    2.28e-48            1.00 
3 ARMA~        0.625   1.03e- 26       NA       NA                   1.05 
4 ARMA~        1.30    5.92e- 77       -0.722    3.33e-32            1.00 
5 ARMA~        1.30    5.92e- 77       -0.722    3.33e-32            1.00 
6 MA2         NA      NA               NA       NA                   0.982
7 MA3         NA      NA               NA       NA                   1.00 
8 MA4         NA      NA               NA       NA                   1.07 
# ... with 9 more variables: intercept_p.value <dbl>, ma1_estimate <dbl>,
#   ma1_p.value <dbl>, ma2_estimate <dbl>, ma2_p.value <dbl>, ma3_estimate <dbl>,
#   ma3_p.value <dbl>, ma4_estimate <dbl>, ma4_p.value <dbl>

Where would you wanna go from here?

Making a nice table out of similar object elements in R?

1 Answers