plot(x,y1,type="l",ylim=c(0,0.13),col="red",main=expression(sigma^2))

运行一下下面的代码，里面有很多格式的。

xlabel <- expression(paste(Delta, italic(s), sep = ""))
ylabel <- expression(alpha[1] * " in (kg)"^2)
plotname <- expression(sin * (beta))
plot(rnorm(50), rnorm(50), xlab = xlabel, ylab = ylabel, main = plotname, xlim = c(-pi, pi), ylim = c(-3, 3), axes = FALSE)
axis(1, at = c(-pi, -pi/2, 0, pi/2, pi), labels = expression(-pi, -pi/2, 0, pi/2, pi))
axis(2)
box()
text(-pi/2, -2, expression(hat(beta) == (X^t * X)^{-1} * X^t * y))
text(pi/2, 2, expression(paste(frac(1, sigma*sqrt(2*pi)), exp*(frac(-(x-mu)^2, 2*sigma^2)), sep = "")), cex = 1.5)

我试了，可是等于3却弄不出来，奇怪

我也试了，但对于我举的例子来说，还是找不到明确的线索，老兄费心了！！

啊！
再次尝试，成功了！
qoiqpwqr 版主，真强！