График испытаний плацебо по всей стране для синтетического контроля в R

Question

Я делаю плацебо по всей стране для своего синтетического дизайна управления.Поскольку у меня 17 стран-доноров, я разработал цикл, который помогает мне сэкономить много времени.Тем не менее, мне нужно сохранить некоторые значения, чтобы позже создать график, подобный следующему:

Текущий код - это (это gaps_birth_i, где я хочу сохранить значения):

birthprep <- rep(NA,16)
i <- 3

for(i in 2:17) {

Synthetic control preparation
birthprep_i <- dataprep(foo = main_donors,
                          predictors = c("gdp","population","labor_women",
                                         "labor_men","enrollment","productivity",
                                         "fertility","marriage","divorce","life",
                                         "gross_income","unemployment"),
                          dependent = "birthrate",
                          unit.variable = "code", 
                          time.variable = "year",
                          treatment.identifier = i, # Select country i as the treated country
                          controls.identifier = unique(main_donors$code)[unique(main_donors$code)!=i],                         time.predictors.prior = c(1992:2001),
                          time.optimize.ssr = c(1992:2001),
                          unit.names.variable = "id",
                          time.plot = 1992:2016)

birthsynth_i <- synth(birthprep_i)

gaps_birth_i <- birthprep_i$Y1plot-(
  birthprep_i$Y0plot%*%birthsynth_i$solution.w )

}

Буду очень признателен, если кто-нибудь поможет мне с этим!

Спасибо.