Игнорирование одного из факторов ggplot в R

Я использую ggplot для построения моей переменной и атрибутов. Я использую ggplot и factor, используя следующий код:

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

результат кода: введите здесь описание изображения

Я хотел бы представить только «Выброс» и исключить из сюжета фактор «Нормальный».

Образец данных:

       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

person Hamed Footohi    schedule 21.01.2013    source источник


Ответы (3)


Вы можете добавить аргумент subset к вашему вызову geom_point и использовать функцию plyr .() для определения подмножества.

eg

p + geom_point(aes(x = HMn25_30.h, y = value, colour = variable), subset = .(Results == 'Outlier'))

Хороший небольшой воспроизводимый пример

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()

введите здесь описание изображения

person mnel    schedule 21.01.2013
comment
Спасибо - это именно то, что я искал :) - person Josh Bode; 22.01.2013

Используйте подмножество фрейма данных, где Results == "Outlier" для команды geom_point:

p <- p + geom_point(data = df.m[df.m$Results == "Outlier",], 
                    cex=9, color= "blue3", shape = 22) 

Тогда нет необходимости в команде scale_shape_manual.

введите здесь описание изображения

Полный код:

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
person Sven Hohenstein    schedule 21.01.2013

Вы не предоставили воспроизводимый код, включающий ваши данные, но вполне вероятно, что вы можете добиться этого с помощью простого подмножества. Попробуйте изменить эту строку:

p <- ggplot(df.m, aes(x = HMn25_30.h, y = value, group = variable,
                                         color = variable))

к этому:

p <- ggplot(df.m[df.m$Result == "Outlier", ], aes(x = HMn25_30.h, y = value,
                                    group = variable, color = variable))

Если я правильно понимаю структуру вашего фрейма данных, это должно отображать только строки, для которых столбец Results имеет значение Outlier.

person SlowLearner    schedule 21.01.2013
comment
Проблема с вашим кодом в том, что он игнорирует места, которые не являются выбросами. Я хотел бы иметь полную строку, но игнорировать только символ нормального - person Hamed Footohi; 21.01.2013
comment
Итак, вы хотите, чтобы Outlier данные отображались в виде линии, но вам не нужен символ для этих точек на линии? - person SlowLearner; 21.01.2013
comment
Нет, я хочу, чтобы и выбросы, и нормальные данные отображались в виде линии. Представление только выбросов через символы. - person Hamed Footohi; 21.01.2013
comment
Возможно, изменить цвет нормальных точек на белый или прозрачный? - person Dennis; 21.01.2013