我有一个要用2个不同的y轴绘制的数据框。一切工作正常,除了类型:第一个图形“ Grundwasserstand”的“线”。但是,选择type ='point'可以将测量结果正确打印在图中。我想问题是由NA造成的。
通常,我需要使type ='line'工作。
我该怎么办?在这种情况下,插值是不可选择的。也不离开Elbewasserstand的数据。
显示了数据框的某些行:
Time Grundwasserstand Elbewasserstand
1 1980-01-01 119.62 110.576
2 1980-01-02 NA 111.426
3 1980-01-03 NA 111.316
4 1980-01-04 NA 111.126
5 1980-01-05 NA 110.956
6 1980-01-06 NA 110.786
7 1980-01-07 NA 110.726
8 1980-01-08 NA 110.866
9 1980-01-09 NA 110.816
10 1980-01-10 NA 110.756
11 1980-01-11 NA 110.656
12 1980-01-12 NA 110.636
13 1980-01-13 NA 110.316
14 1980-01-14 NA 110.176
15 1980-01-15 119.65 110.146
16 1980-01-16 NA 110.166
我的绘图代码是什么样的?´
png(filename = speicherpfad_inkl_namen_png, width = 1000)
par(mar=c(5, 4, 2, 4) + 0.1)
plot(zusm$Time,zusm$Grundwasserstand, type = "b",
col="1", ylab = "Grundwasserstand [m NHN]", xlab = "Time", axes = T)
par(new=T)
plot(zusm$Time, zusm$Elbewasserstand, type = "l",col="blue", lwd=1,
axes=F, ylab = "", xlab = "" )
axis(side = 4, col="black")
mtext(4,text="Elbewasserstand [m NHN]", col = "black", line = 2)
legend("bottomright",title="Legende", legend = c("Grundwasserstand", "Elbewasserstand"),
lwd = 1, col = c(1,4))
dev.off()
更新:
通过以下方式解决:
zusm$Time[!is.na(zusm$Grundwasserstand)],zusm$Grundwasserstand[!is.na(zusm$Grundwasserstand)]
在我的情节代码中,在情节nr中。 1!
png(filename = speicherpfad_inkl_namen_png, width = 1000)
par(mar=c(5, 4, 2, 4) + 0.1)
plot(zusm$Time[!is.na(zusm$Grundwasserstand)],zusm$Grundwasserstand[!is.na(zusm$Grundwasserstand)], type = "l",
col="1", ylab = "Grundwasserstand [m NHN]", xlab = "Time", axes = T)
par(new=T)
plot(zusm$Time, zusm$Elbewasserstand, type = "l",col="blue", lwd=1,
axes=F, ylab = "", xlab = "" )
axis(side = 4, col="black")
mtext(4,text="Elbewasserstand [m NHN]", col = "black", line = 2)
legend("bottomright",title="Legende", legend = c("Grundwasserstand", "Elbewasserstand"),
lwd = 1, col = c(1,4))
dev.off()