Мне нужно создать простой график линии с группами, используя следующие данные:
test = data.frame(x = rep(1:3, each = 2),
group = rep(c("Group 1","Group 2"),3),
groupcd= rep(c(1,2),3),
y= c(22,8,11,4,7,5)
)
Я легко могу сделать это с помощью GGPLOT:
library(ggplot2)
#GGPLOT
qplot(x=x, y=y,
data=test,
colour=group,
main="GGPLOT line plot with groups") +
geom_line()
Я также могу сделать это с помощью TRELLIS:
library(lattice)
xyplot(y~x,
type="b",
group=group,
data=test,
main="TRELLIS line plot with groups",
auto.key =list(
lines = TRUE)
)
Однако, я немного неохотно использую GGPLOT или TRELLIS прямо сейчас. Я хотел бы иметь возможность создать этот график с базой R. Единственный способ, которым я могу получить этот график для работы в Base R, - это использовать для цикла:
# set up empty plot
plot(test$y ~test$x, ylab="y", xlab="x", type="n", main="Base R line plot with groups")
colors<-c("red","blue")
#plot each group in the for loop
number_of_groups <- as.numeric(max(unique(test$groupcd))) #calculate number of groups
for (i in 1:number_of_groups)
{
temp <- subset(test, groupcd==i )
lines(temp$x, temp$y, col=colors[i])
points(temp$x, temp$y, col=colors[i])
}
legend("top", legend=unique(test$group), text.col =colors )
Этот подход кажется довольно запутанным. Есть ли более простой способ сделать это в базе R? Есть ли опция группы в базовой функции R-графика? Большое вам спасибо.