在函数中调用 amelia imputations 以计算多重插补后的 y-hat

Question

我最终想要做的是在多重插补后获得 y-hat 分数，但 Amelia 不提供拟合值。我有使用特定数据集执行此操作的代码，但我正在尝试创建一个无论数据集如何都可以执行此操作的函数。比如：

yhat<-function(a.out,num.obs,num.imp,model.qe){}

其中num.imp 是amelia 或m 中使用的插补数。不过，这就是我卡住的地方：

如果a.out 是amelia 对象并且names(a.out$imputations)[1] 返回"imp1"，为什么a.out$imputations$names(a.out$imputations)[1] 不返回与a.out$imputations$"imp1" 相同的东西？

对于a.out$imputations$names(a.out$imputations)[1]，R 说：Error: attempt to apply non-function。

如何为单个插补和其中的变量创建通用调用？

Answer 1

来自 Amelia 文档的示例

library(Amelia)

data(africa)
a.out <- amelia(x = africa, cs = "country", ts = "year", logs = "gdp_pc")

a.out$imputations 是一个列表

str(a.out$imputations)

列表中的元素可以通过名称或数字索引调用

head(a.out$imputations$imp1)
head(a.out$imputations[["imp1"]])
head(a.out$imputations[[1]])

每个插补都是一个数据框。数据框只是一种特定类型的列表，可以以这种方式调用元素。

head(a.out$imputations[[1]]$country)
head(a.out$imputations[[1]][["country"]])
head(a.out$imputations[[1]][[2]])
head(a.out$imputations[[1]][2])

这些方法中的最后一个方法与其他方法的不同之处在于它返回一个单列数据框，而其他方法返回一个向量。

要完全通过数字索引提取单个值，您可以执行 f.ex

a.out$imputations[[1]][3, 2]

或

a.out$imputations[[1]][[2]][3]