R中的分类变量子集

Question

我如何根据“湖泊”对我的数据进行子集，然后在同一个图上绘制所有三个湖泊的营养物(P，N，C)与日期的关系。另一件事是，我的数据中有空格(NA)，我不想将其替换为零。因为它会显示化学物质的浓度为零，这不是真的。有没有办法使R不为NA值绘制任何点？看起来我在R中得到了看起来奇怪的nA值的图。但excel没有这个问题。(excel根本不会为空白绘制任何内容)

lake                     date    depth        P          N           C
East Long Lake       9/5/1994     0.7      21.9      254.8
East Long Lake       9/5/1994     1.0      30.1     1190.0       257.0
East Long Lake       9/5/1994     1.5      20.5      256.6
East Long Lake       9/5/1994     2.5      22.1      249.0
East Long Lake       9/5/1994    12.0     212.5     2011.6      1090.6
Central Long Lake   6/30/1995     0.6      22.9       91.1
Central Long Lake   6/30/1995     4.0
Peter Lake           7/6/1994     6.1      41.9      527.2        29.6
Peter Lake           7/6/1994    12.0     138.8     1994.0      1409.6
Peter Lake          7/13/1994     0.0      19.1      746.7        22.6 
Peter Lake          7/13/1994     0.7      19.2      21.3

Answer 1

我想这就是你要找的。

library(data.table)
library(ggplot2)
dt = fread("tmp.csv")
dt$date = as.Date(dt$date, "%m/%d/%Y")
dt[,"depth" := NULL]
d <- melt(dt, id.vars=c("lake","date"))

ggplot(d=subset(d, !is.na(value)), aes(date, value, col = variable)) + geom_point() + facet_wrap(~ lake)

这可能就是您要查找的图表：

ggplot(d=subset(d, !is.na(value)), aes(date, value, col = lake)) + geom_point() + facet_wrap(~ variable)