我遇到麻烦ggplot尝试在同一图表上绘制2个不完整的时间序列,其中y数据在x轴(年)上没有相同的值 - 因此,NA存在于某些年份:< / p>
test<-structure(list(YEAR = c(1937, 1938, 1942, 1943, 1947, 1948, 1952,
1953, 1957, 1958, 1962, 1963, 1967, 1968, 1972, 1973, 1977, 1978,
1982, 1983, 1986.5, 1987, 1993.5), A1 = c(NA, 24, NA, 32, 32,
NA, 34, NA, NA, 18, 12, NA, 10, NA, 11, NA, 15, NA, 24, NA, NA,
25, 26), A2 = c(40, NA, 38, NA, 25, NA, 26, NA, 20, NA, 17,
17, 17, NA, 16, 18, 21, 18, 17, 25, NA, NA, 26)), .Names = c("YEAR", "A1",
"A2"), row.names = c(NA, -23L), class = "data.frame")
我尝试下面的代码输出一个脱节的混乱:
ggplot(test, aes(x=YEAR)) +
geom_line(aes(y = A1), size=0.43, colour="red") +
geom_line(aes(y = A2), size=0.43, colour="green") +
xlab("Year") + ylab("Percent") +
scale_x_continuous(limits=c(1935, 1995), breaks = seq(1935, 1995, 5),
expand = c(0, 0)) +
scale_y_continuous(limits=c(0,50), breaks=seq(0, 50, 10), expand = c(0, 0))
我该如何解决这个问题?
答案 0 :(得分:9)
我首选的解决方案是将其重塑为长格式。那么你只需要1 geom_line
个电话。特别是如果你有很多系列,那就更整洁了。与LyzandeR的第二张图相同。
library(ggplot2)
library(reshape2)
test2 <- melt(test, id.var='YEAR')
test2 <- na.omit(test2)
ggplot(test2, aes(x=YEAR, y=value, color=variable)) +
geom_line() +
scale_color_manual(values=c('red', 'green')) +
xlab("Year") + ylab("Percent") +
scale_x_continuous(limits=c(1935, 1995), breaks = seq(1935, 1995, 5),
expand = c(0, 0)) +
scale_y_continuous(limits=c(0,50), breaks=seq(0, 50, 10), expand = c(0, 0))
除了该行之外,您可以考虑添加geom_point()
来电,因此可以清楚哪些点是实际值,哪些点是缺失的。长格式的另一个优点是额外的geom每个只需要1个调用,而不是每个系列1个。
答案 1 :(得分:5)
您可以使用na.omit
删除它们:
library(ggplot2)
#use na.omit below
ggplot(na.omit(test), aes(x=YEAR)) +
geom_line(aes(y = A1), size=0.43, colour="red") +
geom_line(aes(y = A2), size=0.43, colour="green") +
xlab("Year") + ylab("Percent") +
scale_x_continuous(limits=c(1935, 1995), breaks = seq(1935, 1995, 5),
expand = c(0, 0)) +
scale_y_continuous(limits=c(0,50), breaks=seq(0, 50, 10), expand = c(0, 0))
修改
使用2个单独的data.frames与na.omit
:
#test1 and test2 need to have the same column names
test1 <- test[1:2]
test2 <- tes[c(1,3)]
colnames(test2) <- c('YEAR','A1')
library(ggplot2)
ggplot(NULL, aes(y = A1, x = YEAR)) +
geom_line(data = na.omit(test1), size=0.43, colour="red") +
geom_line(data = na.omit(test2), size=0.43, colour="green") +
xlab("Year") + ylab("Percent") +
scale_x_continuous(limits=c(1935, 1995), breaks = seq(1935, 1995, 5),
expand = c(0, 0)) +
scale_y_continuous(limits=c(0,50), breaks=seq(0, 50, 10), expand = c(0, 0))
答案 2 :(得分:0)