Как отличить 4 фактора в ggplot2?

Question

Как отличить 4 разных фактора (не используя размер)? Можно ли использовать полые и сплошные точки для различения переменной в ggplot2?

test=data.frame(x=runif(12,0,1), 
    y=runif(12,0,1), 
    siteloc=as.factor(c('a','b','a','b','a','b','a','b','a','b','a','b')), 
    modeltype=as.factor(c('q','r','s','q','r','s','q','r','s','q','r','s')), 
    mth=c('Mar','Apr','May','Mar','Apr','May','Mar','Apr','May','Mar','Apr','May'), 
    yr=c(2010,2011,2010,2011,2010,2011,2010,2011,2010,2011,2010,2011)) 

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

Answer 1

Вы можете использовать полые и сплошные точки, but only with certain shapes as described in this answer.

Так, что оставляет вас с fill, colour, shape и alpha как эстетических отображений. Это выглядит уродливо, но вот оно:

ggplot(test, aes(x, y, 
       fill=modeltype, 
       shape=siteloc, 
       colour=mth, 
       alpha=factor(yr) 
       )) + 
geom_point(size = 4) + 
scale_shape_manual(values=21:25) + 
scale_alpha_manual(values=c(0.35,1)) 

Уродливое, но я думаю, это то, о чем вы просили. (Я не потрудился выяснить, что происходит с легендой. - это, очевидно, не отображаются границы правой)

Если вы хотите отобразить переменную в виде заказного эстетике (полые и твердые), вам придется пойти немного дальше:

test$fill.type<-ifelse(test$yr==2010,'other',as.character(test$mth)) 
cols<-c('red','green','blue') 

ggplot(test, aes(x, y, 
       shape=modeltype, 
       alpha=siteloc, 
       colour=mth, 
       fill=fill.type 
)) + 
    geom_point(size = 10) + 
    scale_shape_manual(values=21:25) + 
    scale_alpha_manual(values=c(1,0.5)) + 
    scale_colour_manual(values=cols) + 
    scale_fill_manual(values=c(cols,NA)) 

Тем не менее уродливые, но это работает. Я не знаю более чистого способа сопоставления yr с одним цветом, если это 2010 год, и mth, если нет; Я был бы рад, если бы кто-то показал мне более чистый способ сделать это. И теперь гиды (легенда) совершенно ошибочны, но вы можете исправить это вручную.

Answer 2

Думаю, очень сложно визуально различать/сравнивать значения x и y в соответствии с 4 факторами. Я бы использовал огранку, и я уменьшаю количество факторов, используя, например, interaction.

Вот пример с использованием geom_bar:

set.seed(10) 
library(reshape2) 
test.m <- melt(test,measure.vars=c('x','y')) 
ggplot(test.m)+ 
    geom_bar(aes(x=interaction(yr,mth),y=value, 
       fill=variable),stat='identity',position='dodge')+ 
    facet_grid(modeltype~siteloc) 

Answer 3

мне очень нравится с помощью interaction по agstudy - я бы, наверное попробовать это первый. Но если сохранить вещи без изменений, то:

4 фактора могут быть оборудованы огранкой и 2 осями.Тогда есть 2 метрики х и у: один из вариантов является пузырь диаграмма с обе метрики отличить по цвету или форме, или обоих (добавленной джиттера, чтобы сделать формы меньше, перекрывание):

testm = melt(test, id=c('siteloc', 'modeltype', 'mth', 'yr')) 

# by color 
ggplot(testm, aes(x=siteloc, y=modeltype, size=value, colour=variable)) + 
    geom_point(shape=21, position="jitter") + 
    facet_grid(mth~yr) + 
    scale_size_area(max_size=40) + 
    scale_shape(solid=FALSE) + 
    theme_bw() 

#by shape 
testm$shape = as.factor(with(testm, ifelse(variable=='x', 21, 25))) 

ggplot(testm, aes(x=siteloc, y=modeltype, size=value, shape=shape)) + 
    geom_point(position="jitter") + 
    facet_grid(mth~yr) + 
    scale_size_area(max_size=40) + 
    scale_shape(solid=FALSE) + 
    theme_bw() 

# by shape and color 
ggplot(testm, aes(x=siteloc, y=modeltype, size=value, colour=variable, shape=shape)) + 
    geom_point(position="jitter") + 
    facet_grid(mth~yr) + 
    scale_size_area(max_size=40) + 
    scale_shape(solid=FALSE) + 
    theme_bw() 

UPDATE:

Это попытка на основе 1 комментарий от Dominik, чтобы показать, если (х, у) выше или ниже линии 1: 1 и насколько велико отношение х/у или у/х - синий треугольник не является, если х/у> 1, красный круг в противном случае (нет необходимости в melt в данном случае):

test$shape = as.factor(with(test, ifelse(x/y>1, 25, 21))) 
test$ratio = with(test, ifelse(x/y>1, x/y, y/x)) 

ggplot(test, aes(x=siteloc, y=modeltype, size=ratio, colour=shape, shape=shape)) + 
    geom_point() + 
    facet_grid(mth~yr) + 
    scale_size_area(max_size=40) + 
    scale_shape(solid=FALSE) + 
    theme_bw()