У меня есть 40 предметов, двух групп, более чем 15 недель, с некоторой измеренной переменной (Y).
Я хочу иметь график где: x = время, y = T, строки предметами и цветами группами.
Я нашел, что это может быть сделано как это:
TIME <- paste("week",5:20)
ID <- 1:40
GROUP <- sample(c("a","b"),length(ID), replace = T)
group.id <- data.frame(GROUP, ID)
a <- expand.grid(TIME, ID)
colnames(a) <-c("TIME", "ID")
group.id.time <- merge(a, group.id)
Y <- rnorm(dim(group.id.time)[1], mean = ifelse(group.id.time$GROUP =="a",1,3) )
DATA <- cbind(group.id.time, Y)
qplot(data = DATA,
x=TIME, y=Y,
group=ID,
geom = c("line"),colour = GROUP)
Но теперь я хочу добавить к графику что-то для показа различия между этими двумя группами (например, линия тренда для каждой группы, с некоторым CI shadelines) - как это может быть сделано?
Я помню однажды наблюдение, что ggplot2 может (легко) сделать это с geom_smooth, но я пропускаю что-то о том, как заставить его работать.
Кроме того, я размышлял над, возможно, наличием строк быть похожим на коробчатую диаграмму для каждой группы (со строкой для различных квантилей и заборов и так далее). Но я предполагаю отвечать на первый вопрос, помог бы мне разрешить второе.
Спасибо.
p <- ggplot(data=DATA, aes(x=TIME, y=Y, group=ID)) +
geom_line(aes(colour=GROUP)) +
geom_smooth(aes(group=GROUP))
geom_smooth plot http://img143.imageshack.us/img143/7678/geomsmooth.png