我希望用
$P_{M1}(\tilde{z}>z) - P_{M0}(\tilde{z}>z)$
绿线的传奇条目(cptsdtbehavioralm)和
$P_{M2}(\tilde{z}>z) - P_{M0}(\tilde{z}>z)$
蓝线的图例条目(fullbehavioralmodel)替换以下(是的,我在这里用TeX格式写的,只是为了清楚).
这是我生成绘图的代码(我省略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 = …