Как выполнить тестирование панели Root с использованием функции cipstest функции plm Пакет и тест коинтеграции панели с использованием pedroni99m Функция пакета pco

Question

Привет, сообщество stackoverflow, я пытаюсь рассчитать два теста на данных панели. Однако, следуя документации для различных пакетов и соответствующих функций, я получаю ошибки. Мой поиск в Интернете возможных лекарств не принес результатов.

Во-первых, я хочу выполнить тест панели 10 * * второго поколения для данных панели, используя функцию cipstest в пакете plm, но я получить ошибку: Ошибка в прибл. (nintv, cvals [nintl: девятый, tintl, i], n = max (nintv) -: требуется по крайней мере два значения, отличных от NA, для интерполяции

Код, который я попытался использовать для расчета теста:

Y1CT <- cipstest(AT$y1, lags = 2, type = c("trend"), model = c("cmg"), truncated = FALSE)                

Во-вторых, я хочу выполнить тест на интеграцию панели с теми же данными панели, используя функцию pedroni99m в пакете pco, но я получить ошибку: Ошибка в 2: dim (X) [3]: аргумент NA / NaN

Код, который я пытался использовать для вычисления теста:

pedroni99m(PCOY1)

* За 1020 * до использования приведенного выше кода, pedroni99m(PCOY1), я создал новый набор фреймов данных, соответствующих каждой зависимой переменной и ее независимым и управляющим переменным, используя следующие коды:

PCOY1 <- data.frame(AT$y1, AT$x1, AT$x2, AT$x3, AT$x4, AT$x5, AT$x6, AT$x7, AT$z1, AT$z2, AT$z3, AT$z4)


PCOY2 <- data.frame(AT$y2, AT$x1, AT$x2, AT$x3, AT$x4, AT$x5, AT$x6, AT$x7, AT$z1, AT$z2, AT$z3, AT$z4)


PCOY3 <- data.frame(AT$y3, AT$x1, AT$x2, AT$x3, AT$x4, AT$x5, AT$x6, AT$x7, AT$z1, AT$z2, AT$z3, AT$z4)


PCOY4 <- data.frame(AT$y4, AT$x1, AT$x2, AT$x3, AT$x4, AT$x5, AT$x6, AT$x7, AT$z1, AT$z2, AT$z3, AT$z4)

ссылка на мой набор данных, https://www.dropbox.com/s/wazb9kr8wictwjg/mds.csv?dl=0, (сохраняется как * 102 8 * file).

Описание набора данных выглядит следующим образом:

1. Каждый y является зависимой переменной
2. Каждая x является независимой переменной
3. Каждая z является управляющей переменной

Я объявил данные как панель данные следующим образом:

AT <- pdata.frame(mds,index = c('ccode','year'))

Был бы более чем рад получить ваши вдумчивые соображения. Спасибо.