问题
I use ggplot to plot my variable and attributes. I am using ggplot and factor by using the following code:
require(ggplot2)
require(reshape2)
df <- data.frame(HMn25_30$avg,HMn25_30$h)
df[3] = c("Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Normal",
"Outlier",
"Outlier",
"Outlier",
"Outlier",
"Outlier",
"Outlier",
"Normal",
"Outlier",
"Outlier",
"Normal",
"Normal",
"Outlier",
"Outlier",
"Normal",
"Normal"
)
names(df)[1] <- 'Node 25'
names(df)[3] <-'Results'
df.m <- melt(df, names(df)[2:3], names(df)[1])
df.m$Results <- factor(df.m$Results)
df.m$HMn25_30.h <- strptime(as.character(df.m$HMn25_30.h), format = "%Y-%m-%d %H:%M:%S")
p <- ggplot(df.m, aes(x = HMn25_30.h, y = value, group = variable, color = variable))
p <- p + scale_shape_manual(values=c(20,22))
p <- p + geom_point(aes(shape = Results), cex=9, color= "blue3")
p <- p + theme(axis.text.x = element_text(angle = 90, hjust = 1, size=13,color="darkred"))
p <- p + scale_color_manual(values=c("Red"))
p <- p + ylim(-1,8)
p <- p + theme_bw()
p <- p + xlab('Date and Time')
p <- p + ylab('Temprature')
p <- p + ggtitle("Temporal Outliers of Node 25 ") + theme(plot.title = element_text(lineheight=3, face="bold", color="black", size=29))
p <- p + theme(legend.text = element_text(colour="darkred", size = 25))
p <- p + theme(legend.title = element_text(colour="brown", size=25))
p <- p + theme(axis.title.x = element_text(face="bold", colour="darkred", size=16),axis.text.x = element_text(angle=90, vjust=0.5, size=26))
p <- p + theme(axis.title.x = element_text(face="bold", colour="darkred", size=14),axis.text.y = element_text(angle=00, vjust=0.5, size=20))
p <- p + labs(x = "Date-Time [UTC] \ 2007-09-30 ", y = "Temprature ")
p <- p + theme(axis.title.y = element_text(size = rel(2.1), angle = 90))
p <- p + theme(axis.title.x = element_text(size = rel(2.1), angle = 00))
p <- p + geom_line(size=1.9)
p
result of the code:
I would like to present only the 'Outlier' and do omit the 'Normal' factor from plot.
Sample data:
Node 25 HMn25_30.h Results
1 0.26000000 2007-09-29 23:00:00 Normal
2 0.01500000 2007-09-30 00:00:00 Normal
3 -0.35333333 2007-09-30 01:00:00 Normal
4 -0.42333333 2007-09-30 02:00:00 Normal
5 -0.73333333 2007-09-30 03:00:00 Normal
6 -0.65000000 2007-09-30 04:00:00 Normal
7 -0.40000000 2007-09-30 05:00:00 Normal
8 -0.09166667 2007-09-30 06:00:00 Normal
9 0.19000000 2007-09-30 07:00:00 Normal
10 0.63500000 2007-09-30 08:00:00 Normal
11 1.05500000 2007-09-30 09:00:00 Outlier
12 1.26833333 2007-09-30 10:00:00 Outlier
13 2.28166667 2007-09-30 11:00:00 Outlier
14 4.17000000 2007-09-30 12:00:00 Outlier
15 6.34000000 2007-09-30 13:00:00 Outlier
16 6.56666667 2007-09-30 14:00:00 Outlier
17 6.74666667 2007-09-30 15:00:00 Normal
18 5.82833333 2007-09-30 16:00:00 Outlier
19 6.36500000 2007-09-30 17:00:00 Outlier
20 4.60333333 2007-09-30 18:00:00 Normal
21 4.98000000 2007-09-30 19:00:00 Normal
22 2.65666667 2007-09-30 20:00:00 Outlier
23 4.90833333 2007-09-30 21:00:00 Outlier
24 5.05000000 2007-09-30 22:00:00 Normal
25 4.56500000 2007-09-30 23:00:00 Normal
回答1:
You can add a subset argument to your call to geom_point and use the plyr .() function to define the subset.
eg
p + geom_point(aes(x = HMn25_30.h, y = value, colour = variable), subset = .(Results == 'Outlier'))
A nice small reproducible example
DF <- data.frame(a = letters[1:4], b = 1:10)
library(plyr) # must be explicitly loaded
ggplot(DF, aes(x = b, y = b)) +
geom_point(subset = .(a == 'a'), colour = 'blue') +
geom_point(subset = .(a == 'c'), colour = 'green') +
geom_line()
回答2:
Use the subset of the data frame where Results == "Outlier" for the geom_point command:
p <- p + geom_point(data = df.m[df.m$Results == "Outlier",],
cex=9, color= "blue3", shape = 22)
Then, there's no need for the scale_shape_manual command.
The complete code:
library(ggplot2)
p <- ggplot(df.m,
aes(x = HMn25_30.h, y = value, group = variable, color = variable))
#p <- p + scale_shape_manual(values=c(20,22)) # command removed
p <- p + geom_point(data = df.m[df.m$Results == "Outlier",],
cex=9, color= "blue3", shape = 22) # command modified
p <- p + theme(axis.text.x = element_text(angle = 90, hjust = 1, size=13,
color="darkred"))
p <- p + scale_color_manual(values=c("Red"))
p <- p + ylim(-1,8)
p <- p + theme_bw()
p <- p + xlab('Date and Time')
p <- p + ylab('Temprature')
p <- p + ggtitle("Temporal Outliers of Node 25 ") +
theme(plot.title = element_text(lineheight=3, face="bold",
color="black", size=29))
p <- p + theme(legend.text = element_text(colour="darkred", size = 25))
p <- p + theme(legend.title = element_text(colour="brown", size=25))
p <- p + theme(axis.title.x = element_text(face="bold",colour="darkred" size=16),
axis.text.x = element_text(angle=90, vjust=0.5, size=26))
p <- p + theme(axis.title.x = element_text(face="bold",colour="darkred",size=14),
axis.text.y = element_text(angle=00, vjust=0.5, size=20))
p <- p + labs(x = "Date-Time [UTC] \ 2007-09-30 ", y = "Temprature ")
p <- p + theme(axis.title.y = element_text(size = rel(2.1), angle = 90))
p <- p + theme(axis.title.x = element_text(size = rel(2.1), angle = 00))
p <- p + geom_line(size=1.9)
p
回答3:
You have not supplied reproducible code that includes your data, but it seems likely you can achieve this with simple subsetting. Try changing this line:
p <- ggplot(df.m, aes(x = HMn25_30.h, y = value, group = variable,
color = variable))
to this:
p <- ggplot(df.m[df.m$Result == "Outlier", ], aes(x = HMn25_30.h, y = value,
group = variable, color = variable))
If I understand the structure of your data frame correctly that should plot only the rows for which the Results column has the value Outlier.
来源：https://stackoverflow.com/questions/14431246/ignoring-one-of-the-factors-of-ggplot-in-r