I tried to use substitute in legend as follows: pval <- 0.04 plot(0) legend(1,0.5,substitute(hat(theta) == p, list(p = pval))) For some reason the legend is repeated 3 times. Any suggestions or is this a bug?