I am getting an error when I try and combine using expression with do.call and plot.
expression
do.call
plot
x <- 1:10 y <- x^1.5
do.call evaluates the parameters before running the function; try wrapping the expression in quote:
quote
p <- list(xlab=quote(expression(paste("Concentration (",mu,"M)")))) do.call("plot", c(y~x, p))