如何从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。