Множественное вменение с Амелией - улучшение результатов - PullRequest
0 голосов
/ 16 апреля 2020

Я вменяю значения АН в Амелию, и результаты не велики. Мне нужно вменять НС в отношении иностранных инвестиций (ПИИ) для некоторых регионов.

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

head(df_nominal)

Набор данных

Когда я включаю границы от 0 до бесконечности, результаты завышаются. Если я не включу их, это даст мне отрицательные значения. См., Например, результаты для «La Rioja»:

bds <- matrix(c(4, 0, Inf), nrow = 1, ncol = 3)

a.out <- amelia(df_nominal, m=5, ts="Year", cs="Region", p2s=0, polytime = 1, intercs = TRUE,lags=c("Price","FDI"), leads=c("Price","FDI"),bounds = bds)
a.out1 <- amelia(df_nominal, m=5, ts="Year", cs="Region", p2s=0, polytime = 1, intercs = TRUE,lags=c("Price","FDI"), leads=c("Price","FDI"))

tscsPlot(a.out, cs = "La Rioja", main = "La Rioja", var = "FDI")
tscsPlot(a.out1, cs = "La Rioja", main = "La Rioja", var = "FDI")

Результаты с границами Результаты без границ

Я добавил несколько приоров с максимальными значениями с уровнем достоверности почти 1, но результаты все еще превышают максимальное значение моих априоров

prs <- matrix(c(13,31,49,67,121,139,157,4,4,4,4,4,4,4,0,0,0,0,0,0,0,1499,1499,1499,1499,1499,1499,1499,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999,0.9999),nrow=7, ncol=5)

bds <- matrix(c(4, 0, Inf), nrow = 1, ncol = 3)

a.out <- amelia(df_nominal, m=5, ts="Year", cs="Region", p2s=0, polytime = 1, intercs = TRUE,lags=c("Price","FDI"), leads=c("Price","FDI"),bounds = bds, priors=prs)

tscsPlot(a.out, cs = "La Rioja", main = "La Rioja", var = "FDI")

Как вы можете видеть на графике, результаты превышают определенный мной диапазон и являются вряд ли это будет реалистично c. Результаты с оценками и априорами

Можете ли вы дать мне несколько советов о том, как улучшить результаты?

Заранее спасибо!

...