我正在使用ggplot2叠加几个密度图来创建图形。假设我有每年12月的数据点(在这种情况下是2004-2012年),我想为每个12月份+年份绘制一个密度函数并将它们叠加在一起。我想通过将其中一个月的密度线变成虚线来突出显示其中一个密度线,而其他所有密度线都是实线。以下是一个可复制的示例:
#Create vector of data for December
set.seed(12345)
dec_emas = as.matrix(rnorm(496, 122, 250))
#create indicators for Dec04 ... Dec11, then attach to data frame w/ estimates
declab = c('Dec04', 'Dec05', 'Dec06', 'Dec07', 'Dec08', 'Dec09', 'Dec10', 'Dec11')
declabs = rep(declab, 62)
rownames(dec_emas) = declabs
colnames(dec_emas) = 'EMA'
#add in factor ID for the 8 levels
dec04 = as.numeric(rownames(dec_emas) == 'Dec04')
dec05 = as.numeric(rownames(dec_emas) == 'Dec05')
dec06 = as.numeric(rownames(dec_emas) == 'Dec06')
dec07 = as.numeric(rownames(dec_emas) == 'Dec07')
dec08 = as.numeric(rownames(dec_emas) == 'Dec08')
dec09 = as.numeric(rownames(dec_emas) == 'Dec09')
dec10 = as.numeric(rownames(dec_emas) == 'Dec10')
dec11 = as.numeric(rownames(dec_emas) == 'Dec11')
dec_ema = as.numeric(dec_emas)
dec_ema_df = data.frame(dec_ema, dec04,dec05,dec06,dec07,dec08,dec09,dec10,dec11)
dec_ema_df$new[dec_ema_df$dec04==1]<-4
dec_ema_df$new[dec_ema_df$dec05==1]<-5
dec_ema_df$new[dec_ema_df$dec06==1]<-6
dec_ema_df$new[dec_ema_df$dec07==1]<-7
dec_ema_df$new[dec_ema_df$dec08==1]<-8
dec_ema_df$new[dec_ema_df$dec09==1]<-9
dec_ema_df$new[dec_ema_df$dec10==1]<-10
dec_ema_df$new[dec_ema_df$dec11==1]<-11
#using GGPLOT2...
library(ggplot2)
library('RColorBrewer')
#use different colours for each december
cbbPalette <- brewer.pal(8, "Accent")
#Plot overlaid densities and highlight December 08 with a dashed line.
dec.dens1 = ggplot(dec_ema_df, aes(x=dec_ema, group = new)) +
geom_density(aes(colour=factor(new), linetype = factor(new==8)), size = 1, alpha = .60, kernel='epanechnikov') +
xlab(NULL) + ylab(NULL) +
scale_linetype_manual(values = c('solid', 'dashed'), guide = FALSE) +
scale_colour_manual(values = cbbPalette,
labels= c('Dec `04','Dec `05','Dec `06','Dec `07','Dec `08','Dec `09','Dec `10','Dec `11'),
name= '')
print(dec.dens1)
目前,我能够生成这个。不幸的是,右侧的图例显示了一个实线蓝色线条,而不是虚线。如何使图例与仅在2008年12月的虚线密度线相对应(例如,所有其他图例键都是实线)?我尝试使用legend.key,但这仅适用于轮廓和背景,而不适用于图例中的实际对象。