"Неверное количество измерений" при запуске Zelig "арима" на вмененных данных

Question

Я получаю сообщение об ошибке при попытке запустить arima модель с zelig пакет. Я использую данные МИ с 20 вменениями, которые были созданы с помощью Амелии. Вот краткое резюме моих переменных id и response:

 $ imp20:'data.frame':  442 obs. of  50 variables:
  ..$ region                     : Factor w/ 4 levels "Central Africa",..: 3     3 3 3 3 3 3 3 3 3 ...
  ..$ subregionid                : Factor w/ 4 levels "FC","FE","FS",..: 3 3     3 3 3 3 3 3 3 3 ...
  ..$ country                    : Factor w/ 34 levels "Angola","Benin",..: 1 1 1 1 1 1 1 1 1 1 ...
  ..$ ISO2                       : Factor w/ 34 levels "AO","BF","BJ",..: 1 1 1 1 1 1 1 1 1 1 ...
  ..$ ISO3                       : Factor w/ 34 levels "AGO","BEN","BFA",..: 1 1 1 1 1 1 1 1 1 1 ...
  ..$ year                       : num [1:442] 2002 2003 2004 2005 2006 ...
  ..$ cap.lat                    : num [1:442] -8.5 -8.5 -8.5 -8.5 -8.5 -8.5 -8.5 -8.5 -8.5 -8.5 ...
  ..$ cap.long                   : num [1:442] 13.2 13.2 13.2 13.2 13.2 ...
  ..$ NGDP_RPCH                  : num [1:442] 14.53 5.25 10.88 18.26 20.73 ...
  ..$ NGDPD                      : num [1:442] 3.18 3.31 3.38 3.44 3.48 ...
  ..$ NGDPDPC                    : num [1:442] 2.68 2.69 2.72 2.75 2.78 ...
  ..$ NGSD_NGDP                  : num [1:442] 10.62 7.77 12.63 26.98 40.94     
...
 ..$ PIKE.regional              : num [1:442] 0.225 0.295 0.287 0.358 0.357 ...
  ..$ Definite.Probable          : num [1:442] 36 36 36 36 36.1 ...
  ..$ Elephant.range             : num [1:442] 406006 433613 511662 456046 459418 ...
  ..$ Change.by.year             : num [1:442] 0.000463 0.000463 0.000463 0.000463 0.000463 ...
  ..$ Diff.from.expected         : num [1:442] -0.0415 -0.0415 -0.0415 -0.0415 -0.0415 ...

Diff.from.expected моя переменная ответа И вот код, который я запустил вместе с ошибкой, которую я получаю.

z1 <- zarima$new()
> z1$zelig(Diff.from.expected~GNI, order=c(1,0,1), model="arima",

+                data = a.coVarsTrans.more, ts="year", cs="country")

Error in data[, cs] : incorrect number of dimensions

Так что мне кажется, что есть проблема с cs='country' позвоните, но я не уверен, в чем проблема. Я планирую добавить больше независимых переменных, но хочу убедиться, что базовая модель работает в первую очередь, что явно не работает.

Вот ссылка на мой сохраненный файл Амелии.Rdata.