如何从 r svydesign 包中保存加权数据？

Question

有没有办法保存加权数据，以便下次我直接加载它？我有一个大型调查数据集，并使用 R 的 survey 包。我使用feather 包加载数据，但应用svydesign 需要相当长的时间。这是一个可重现的示例：

df <- data.frame(col1 = rnorm(20, 0, 1), col2 = rnorm(20, 2, 2), w = rnorm(20, 1, .2))
df.w <- svydesign(id = ~1, data = df, weights = ~w)

我想保存 df.w 并将其用于将来的分析。有什么办法吗？

Answer 1

您可以像这样使用saveRDS 和readRDS 保存/读取单个R 对象。

library(survey)

df <- data.frame(col1 = rnorm(20, 0, 1), col2 = rnorm(20, 2, 2), w = rnorm(20, 1, .2))
df.w <- svydesign(id = ~1, data = df, weights = ~w)

####### save to file ##########
storage_file <- tempfile()
#storage_file <- "mydesign.rds"   ## uncomment here to use a local file
saveRDS(df.w, storage_file)

######## clear workspace #########
rm(df, df.w)

######### load the data ###########
df.w.loaded <- readRDS(storage_file)

df.w.loaded
## Independent Sampling design (with replacement)
## svydesign(id = ~1, data = df, weights = ~w)

######## delete storage file ######
file.remove(storage_file)

如果您想将多个对象保存在一个文件中，请查看?save 和?load。