13

我希望将以下内容（是的，我在这里以 TeX 格式编写，只是为了清楚
$P_{M1}(\tilde{z}>z) - P_{M0}(\tilde{z}>z)$
起见）替换绿线的图例条目（cptsdtbehavioralm）和
$P_{M2}(\tilde{z}>z) - P_{M0}(\tilde{z}>z)$
蓝线的图例条目（fullbehavioralmodel）。

这是我生成绘图的代码（我省略了 10,000 次观测数据集和生成函数 Fm0、Fm1 和 Fm2 的转换）：

bmp("bias_plot_v4.bmp", width=540, pointsize=10)
ggplot(data.frame(x=c(0,80)),aes(x) ) +
   stat_function(fun=function(x)((1-Fm1(x)) - (1- Fm0(x))), geom="line", 
aes(colour="cptsdtbehavioralm"), n=1000) +
   stat_function(fun=function(x)((1-Fm2(x)) - (1- Fm0(x))), geom="line", 
aes(colour="fullbehavioralmodel"), n=1000) +
 theme_bw() +
   ylab("Probability") +
 xlab(expression(paste("Security breaches per 3-week spear-phishing campaign ", 
italic( (z) )))) +  
 theme(aspect.ratio=.618) + 
   theme(legend.position=c(0.845,0.8)) +
   theme(legend.key = element_blank()) +
   scale_color_manual(values = c("cptsdtbehavioralm"="green2", 
"fullbehavioralmodel" = "blue"), name="Bias")
dev.off()

每当我尝试用数学表达式替换“cptsdtbehavioralm”和“fullbehavioralmodel”时——例如为了简单起见expression(P[{M1}])——我会收到以下类型的错误：

Error: unexpected '=' in:
"       theme(legend.key = element_blank()) +
   scale_color_manual(values = c(expression(P[{M1}]))="

我已经束手无策了——任何见解或建议都会非常非常受欢迎。

4

1 回答 1

33

这是一个简单的例子，说明如何使用values和labels获得我认为你想要的东西。

请注意，values将数据中的值映射到您希望使用的颜色，同时labels是您希望显示的标签（所以这是您放置表达式的地方）。

ggplot(data = data.frame(x= c(0,5)),aes(x=x)) +
    stat_function(fun=dnorm,aes(colour = 'red')) + 
    stat_function(fun = dexp, aes(colour = 'blue')) + 
    scale_colour_manual(values = c('red' = 'red','blue' = 'blue'),name = '', 
          labels = expression(P[M1](tilde(z)>0),P[M0](tilde(z)>0)))

在此处输入图像描述

于 2013-10-22T02:16:38.333 回答