Мне нужно запустить тест стационарности для нескольких временных рядов.Я использую функцию ur.df из пакета urca, чтобы сделать их.Затем я сохраняю результаты каждого из этих adf тестов в виде списка списков, поскольку каждый из них является самим списком.
Мне нужно иметь возможность хранить все параметры, как от вывода в кадре данных.Есть ли способ сделать это?
Я знаю, что могу извлечь некоторые параметры, такие как @teststat, @cval и т.п., но как мы можем получить все параметры, как мы можем сlm вывод регрессии с использованием broom
Например, если это выход из ur.df function
test1 <- ur.df(usage_1601_1612, type = "none", lags = 1, selectlags = "AIC")
Содержимое test1 показано ниже
Test regression none
Call:
lm(formula = z.diff ~ z.lag.1 - 1 + z.diff.lag)
Residuals:
Min 1Q Median 3Q Max
-6093.2 -1385.8 -100.9 1414.3 6962.8
Coefficients:
Estimate Std. Error t value Pr(>|t|)
z.lag.1 -0.004212 0.005191 -0.811 0.4177
z.diff.lag -0.126685 0.052161 -2.429 0.0156 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 2351 on 362 degrees of freedom
Multiple R-squared: 0.01838, Adjusted R-squared: 0.01296
F-statistic: 3.39 on 2 and 362 DF, p-value: 0.03479
Value of test-statistic is: -0.8114
Critical values for test statistics:
1pct 5pct 10pct
tau1 -2.58 -1.95 -1.62
Это общий код, который я запускаю
urresultorigobjects <-
lapply(usagextsobjects, function(x) {
summary(ur.df(x, type = "none", lags = 1, selectlags = "AIC"))
})
. Он генерирует список urresultorigobjects, содержащий выходные данные теста стационарности для каждого из xts объектов в списке usagextsobjects.
Я хотел бы иметь возможность собрать все выходные данные и сохранить соответствующую информацию во фрейме данных для справки на основе кода для последующей обработки.
Любые предложения будут с благодарностью.