Я хочу иметь возможность генерировать таблицу коррелограмм в стиле Stata, например, созданную командой corrgram
, но в R: ![Stata output](https://i.stack.imgur.com/Lr6ZH.png)
Мне удалось сгенерировать что-токоторый аппроксимирует первые четыре столбца, используя этот код:
max_lags = 20
acf_vec <- acf(z, plot=FALSE, lag.max=max_lags)$acf
pacf_vec <- pacf(z, plot=FALSE, lag.max=max_lags)$acf
Q_stats <- c()
Q_pvals <- c()
for (i in 1:max_lags) {
Q = Box.test(z, lag=i)
Q_stats[[i]] <- Q$statistic
Q_pvals[[i]] <- Q$p.value
}
corrgram <- cbind(LAG=seq(1,max_lags), ACF=acf_vec, PAC=pacf_vec, Q=Q_stats, "Prob>Q"=Q_pvals)
kable(corrgram)
Это приводит к таблице ниже.Проблема в том, что я не могу понять, как создавать мини-сюжеты в последних двух столбцах.
![enter image description here](https://i.stack.imgur.com/Zejqb.png)
У кого-нибудь есть какие-либо предложения о том, как создавать столбцы, подобные этим R?Или еще лучше, если есть эквивалентная функция?
Примечание 1: В R есть пакет с именем corrgram
, но он создает цветные графики, а не простую таблицу.
Примечание 2: Этот пост имеет тот же вопрос, но нет удовлетворительного ответа.Я не хочу делать графики для ACF и PACF, я хочу видеть саму таблицу.Любая помощь могла бы быть полезна.