Использование функции Forex () для объекта plm дает предиктивные значения для каждой сущности в модели с фиксированными эффектами. 1003 *
Приведен пример:
library(plm)
data("Produc", package = "plm")
zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
data = Produc, index = c("state","year"))
yy=predict(zz)
# this returns only the first element of the ARIZONA entity prediction.
yy['ARIZONA']
#Output from above code
> yy['ARIZONA']
ARIZONA
-0.4264009
#Expected output:
> yy
ALABAMA ALABAMA ALABAMA ALABAMA ALABAMA
-0.2361898231 -0.2189895892 -0.1696057318 -0.1115933089 -0.0837388837
ALABAMA ALABAMA ALABAMA ALABAMA ALABAMA
-0.0927599459 -0.0162445537 0.0236004000 0.0771475623 0.0973429273
ALABAMA ALABAMA ALABAMA ALABAMA ALABAMA
0.0956672093 0.0784827429 0.0501708544 0.0640040527 0.1157573807
ALABAMA ALABAMA ARIZONA ARIZONA ARIZONA
0.1541848055 0.1727639012 -0.4264009370 -0.3666204604 -0.2707038121
ARIZONA ARIZONA ARIZONA ARIZONA ARIZONA
-0.1809125087 -0.1410211123 -0.1802185768 -0.1477399971 -0.0839823046
ARIZONA ARIZONA ARIZONA ARIZONA ARIZONA
0.0138358070 0.0985224005 0.1273115184 0.1711627777 0.1495094189
ARIZONA ARIZONA ARIZONA ARIZONA
0.1919410264 0.2873534379 0.3458664549 0.4120968673
Мои вопросы: как мне манипулировать объектом pred (), чтобы получить полный набор прогнозируемых значений только для одной из моих сущностей?