具有时间点的分箱纵向数据的意大利面条图

Question

我有一个包含纵向数据的数据框，如下所示

测试

    Names     hr1    hr2    hr3    hr4      workhr_bin
      41     80     76     70     60               7
      42     80     74     75     NA               8
      43     85     NA     60     65               6
      44     NA     NA     NA     60               3
      45     80     70     NA     NA               8
      46     NA     NA     NA     60               3

hr1、hr2、hr3、的变量hr4包括“姓名”列下的受试者以重复的时间间隔报告的服务时间。“workhr_bin”列包括使用分位数函数获得的 bin。总共有 10 个垃圾箱，1:10。

我正在尝试生成多个由垃圾箱组成的小时意大利面图。本质上它应该产生 10 个图，一个图用于 bin 1 中的数据，另一个图用于数据 bin 2 等。

我尝试这样做：

head(melt(test[,c(2:6)]))

但我最终得到的workhr_bin变量消失了，取而代之的是一个像这样的文件：

  variable value
1   hours1    80
2   hours1    80
3   hours1    85
4   hours1    NA
5   hours1    80
6   hours1    NA

我也尝试过

melt(test, id.var = "Names")

我最终变量workhr_bin也消失了

   Names variable value
    41   hr1      80.00
    42   hr1      80.00
    43   hr1      85.00

我尝试使用latticeandggplot2但由于某种原因，我无法将数据转换为正确的格式来生成代表每个箱内样本的意大利面 10 个图。

本质上我需要一个数据框：

Names    variable     value    workhr_bin
    41   hr1      80.00             7
    42   hr1      80.00             8
    43   hr1      85.00             6

然后我希望能够创建一个 bin-faceted 多色意大利面条图，其中 x 轴上有“变量”（由hr1, hr2, hr3,组成hr4），y 轴上有相应的“值”。

Answer 1

您应该研究 的文档reshape2::melt。

DF <- read.table(text="    Names     hr1    hr2    hr3    hr4      workhr_bin
      41     80     76     70     60               7
      42     80     74     75     NA               8
      43     85     NA     60     65               6
      44     NA     NA     NA     60               3
      45     80     70     NA     NA               8
      46     NA     NA     NA     60               3", header=TRUE)

library(reshape2)
DF_melt <- melt(DF, id.vars=c("Names", "workhr_bin"))
#make time numeric
DF_melt$variable <- as.numeric(gsub("hr", "", DF_melt$variable))

library(ggplot2)
p <- ggplot(DF_melt, aes(x=variable, y=value, color=factor(Names))) +
  geom_line() +
  geom_point() +
  facet_grid(workhr_bin ~ .)

print(p)

（这会导致 ggplot2 发出警告，因为没有足够的数据。我假设您的真实数据集不会出现此问题。）