Par3 <- c(4.27, 2.90, 0.0868, 0.0224, -0.0287)
xx1 <- seq(-4.27, 4.27, by=0.001)
yy1 <- EPE(P=Par3, x=xx1, simpver=NULL)
yy2 <- -EPE(P=Par3, x=xx1, simpver=NULL)
dev.new()
plot(xx1, yy1, asp=1, type="l", col=4, cex.lab=1.5, cex.axis=1.5,
xlim=c(-5, 5), ylim=c(-5, 5),
xlab=expression(italic(x)), ylab=expression(italic(y)))
lines(xx1, yy2, col=2)
graphics.off()
Run the code above in your browser using DataLab