Построение двух векторов данных на графике рассеяния GGPLOT2 с использованием R

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

в lattice я мог бы сделать это:

xyplot(Predicted_value + Actual_value ~ x_value | State_CD, data=dd)

и это дало бы мне панель для каждого State_CD с каждым столбцом

Я могу сделать один столбец с ggplot2:

pg <- ggplot(dd, aes(x_value, Predicted_value)) + geom_point(shape = 2) 
      + facet_wrap(~ State_CD) + opts(aspect.ratio = 1)
print(pg)

Я не могу понять, как добавить Actual_value в ggplot выше.

РЕДАКТИРОВАТЬ Хэдли указал, что это действительно будет проще с воспроизводимым примером. Вот код, который, кажется, работает. Есть ли лучший или более краткий способ сделать это с помощью ggplot? Почему синтаксис для добавления другого набора точек в ggplot так отличается от добавления первого набора данных?

library(lattice)
library(ggplot2)

#make some example data
dd<-data.frame(matrix(rnorm(108),36,3),c(rep("A",24),rep("B",24),rep("C",24)))
colnames(dd) <- c("Predicted_value", "Actual_value", "x_value", "State_CD")

#plot with lattice
xyplot(Predicted_value + Actual_value ~ x_value | State_CD, data=dd)

#plot with ggplot
pg <- ggplot(dd, aes(x_value, Predicted_value)) + geom_point(shape = 2) + facet_wrap(~ State_CD) + opts(aspect.ratio = 1)
print(pg)

pg + geom_point(data=dd,aes(x_value, Actual_value,group=State_CD), colour="green")

Выход решетки выглядит следующим образом:
_{(источник: cerebralmastication.com)}

и ggplot выглядит так:
_{(источник: cerebralmastication.com)}

Ответ 1

Просто следуя тому, что предложил Ян: для ggplot2 вы действительно хотите, чтобы все элементы оси Y в одном столбце с другим столбцом в качестве фактора, указывающего, как вы хотите его украсить. Это легко сделать с помощью melt. Для остроумия:

qplot(x_value, value, 
      data = melt(dd, measure.vars=c("Predicted_value", "Actual_value")), 
      colour=variable) + facet_wrap(~State_CD)

Вот как это выглядит для меня:
_{(источник: princeton.edu)}

Чтобы понять, что на самом деле делает melt, вот глава:

> head(melt(dd, measure.vars=c("Predicted_value", "Actual_value")))
     x_value State_CD        variable      value
1  1.2898779        A Predicted_value  1.0913712
2  0.1077710        A Predicted_value -2.2337188
3 -0.9430190        A Predicted_value  1.1409515
4  0.3698614        A Predicted_value -1.8260033
5 -0.3949606        A Predicted_value -0.3102753
6 -0.1275037        A Predicted_value -1.2945864

Видите ли, он "плавит" Predicted_value и Actual_value в один столбец с именем value и добавляет еще один столбец с именем variable позволяющий узнать, из какого столбца он изначально пришел.

Ответ 2

Обновить: через несколько лет я почти всегда использую метод Джонатана (через пакет tidyr) с ggplot2. Мой ответ ниже работает в крайнем случае, но получает утомительно быстро, когда у вас есть переменные 3+.

Я уверен, что у Хэдли будет лучший ответ, но - синтаксис отличается, потому что синтаксис ggplot(dd,aes()) (я думаю) в первую очередь предназначен для построения только одной переменной. Для двоих я бы использовал:

ggplot() + 
geom_point(data=dd, aes(x_value, Actual_value, group=State_CD), colour="green") + 
geom_point(data=dd, aes(x_value, Predicted_value, group=State_CD), shape = 2) + 
facet_wrap(~ State_CD) + 
theme(aspect.ratio = 1)

Вытягивание первого набора точек из ggplot() дает ему тот же синтаксис, что и второй. С этим легче справиться, потому что синтаксис один и тот же, и он подчеркивает "Грамматику графики", которая лежит в основе ggplot2.

Ответ 3

вы можете просто немного изменить форму своих данных, так что у вас есть одна переменная оси y с дополнительной факторной переменной, указывающей, является ли она прогнозируемой или фактической переменной.

Это что-то вроде того, что вы пытаетесь сделать?

dd<-data.frame(type=rep(c("Predicted_value","Actual_value"),20),y_value=rnorm(40),
                x_value=rnorm(40),State_CD=rnorm(40)>0)
qplot(x_value,y_value,data=dd,colour=type,facets=.~State_CD)

Ответ 4

после публикации вопроса я наткнулся на этот раздел справки R, который, возможно, помог мне. Похоже, я могу это сделать:

 pg + geom_line(data=dd,aes(x_value, Actual_value,group=State_CD), colour="green")

- это хороший способ сделать что-то? Это странно для меня, потому что добавление второго элемента имеет совершенно другой синтаксис, чем первый.