从R中的推算数据绘制多个箱线图

Question

我有一个正在分析的估算数据集，我正在尝试绘制箱线图，但我无法理解正确的过程。

我的数据（一个样本，原始每个插补有 20 个观察值，每组有 13 个变量，所有值范围从 0 到 25）：

.imp  .id   FTE_RM  FTE_PD  OMZ_RM  OMZ_PD
1     1     25      25      24      24
1     2     4       0       2       6
1     3     11      5       3       2
1     4     12      3       3       3
2     1     20      15      15      15
2     2     4       1       2       3
2     3     0       0       0       6
2     4     20      0       0       0

.imp 表示插补轮次，.id 表示每次观察的标识符。

我想在一个图中绘制所有 FTE_* 变量（以及另一个图中的 `OMZ_*），但想知道如何处理所有插补，我可以只包含所有值吗？插补数据现在有 500 个观测值。例如，使用方差分析，我需要将方差分析结果平均 5 以返回 20 个观察值。但这是否也需要用于箱线图，因为我只处理中位数、均值、最大值。和分钟？

如：

data_melt <- melt(df[grep("^FTE_", colnames(df))])
ggplot(data_melt, aes(x=variable, y=value))+geom_boxplot()

我玩过几次 ggplot，但我认为自己是个新手。

Answer 1

我假设你想在融化后保留 .imp 和 .id 的标识符，所以宁愿放：

data_melt <- melt(df,c(".imp",".id"))

为了数据框的完整性，引入一个标识类型的列可能会有所帮助 - FTE 与 OMZ：

data_melt$type <- ifelse(grepl("FTE",data_melt$variable),"FTE","OMZ")

有了这个data.frame，你可以例如对类型进行分面（或者你可以在data_melt上使用一个简单的过滤器语句来限制一个类型）：

ggplot(data_melt, aes(x=variable, y=value))+geom_boxplot()+facet_wrap(~type,scales="free_x")

这看起来像这样。 编辑：修复数据混乱