General questions regarding the plot

Question

Using R , I drew a shaded plot similar to this one . I want to do the following 4 things in R

• Add a legend as shown in the link.
• Replace the x-axis name with the Greek delta li>
• To indicate different intersection points on the plot. for example, with x = 0.75, several curves correspond to the x axis, I want to put a value of 0.75 next to this point.
• If you see curves, they are not smooth. How to make them smooth? Even great stories are much smoother.

How to do it?

Here is the plot.

The following code is used to draw a graph.

 plot(NA,xlim=c(0,1),ylim=c(0,1),xlab="delta",ylab="K", xaxs="i",yaxs="i") # Empty plot a1 <- curve((x+x^7-x^2-x^4)/(1+xx^3-x^4), from=0, n=450000, add = TRUE) # First curve a2 <- curve((x^2+x^3-xx^5)/(x+x^2), to=0.9, n=45000, add = TRUE) a3 <- curve((x+x^7-x^2-x^4)/(1+xx^2-x^3-x^4+x^7),from=0, n=45000, add = TRUE) a4 <- curve((x+x^8-x^3-x^5)/(x+x^8-x^3-x^5+1-x^4),from=0, to=0.9, n=45000, add = TRUE) a5 <- curve((x+x^8-x^3-x^5)/(1+xx^5-x^4),from=0, n=45000, add = TRUE) a6 <- curve((x+x^2-x^4-1)/(xx^4), to=0.84, n=45000, add = TRUE) a7 <- curve((x+x^6-x^3-x^4)/(1+xx^3-x^4), from=0.83 ,to=1, n=45000, add = TRUE) a8 <- curve((1+x^7-x^2-x^4)/(1+x^3-xx^4), from=0.819, n=45000, add = TRUE) a9 <- curve((x)/(1+x), n=45000,from=0.819, to =1, add = TRUE) names(a1) <- c('xA1','yA1') names(a2) <- c('xA2','yA2') names(a3) <- c('xA3','yA3') names(a4) <- c('xA4','yA4') names(a5) <- c('xA5','yA5') names(a6) <- c('xA6','yA6') names(a7) <- c('xA7','yA7') names(a8) <- c('xA8','yA8') names(a9) <- c('xA9','yA9') with(as.list(c(a1,a2,a3,a4,a5,a6,a7,a8,a9)),{ idA <- yA3 >=0 idB <- yA2 >=0 & yA2 <= yA4 idC <- yA4 >= yA2 idD <- yA5 >=0 idE <- yA6 >=0 & yA6 <= yA7 idF <- yA7 <= yA6 idG <- yA8 >=0 & yA8 <= yA9 idH <- xA9 >= xA8 & xA9 >0.8 idI <- xA1 >=0 & xA1 <= 0.755 idJ <- xA3 >=0 & xA3 <= 0.755 polygon(x = c(xA3[idA],xA2[idB],rev(xA4[idC])), y = c(yA3[idA],yA2[idB],rev(yA4[idC])), density=20, angle=90, border=NULL) polygon(x = c(xA5[idD],1,1,0), y = c(yA5[idD],0,1,1), density=20, angle=0, border=NULL) polygon(x = c(xA6,xA7), y = c(yA6,yA7), density=20, angle=45, border=NULL) polygon(x = c(rev(xA8[idG]),xA9[idH],1), y = c(rev(yA8[idG]),yA9[idH],0), density=20, angle=135, border=NULL) polygon(x = c(xA1[idI],rev(xA3[idJ])), y = c(yA1[idI],rev(yA3[idJ])), col="black", border=NULL) }) 

Answer 1

 layout(matrix(c(1,2),nrow=1), width=c(4,1)) #Divide your plotting region in two inequal part par(mar=c(5,4,4,0)) #Get rid of the margin on the right side plot(NA,xlim=c(0,1),ylim=c(0,1), xlab=expression(delta),ylab="K", xaxs="i",yaxs="i") # Here your delta a1 <- curve((x+x^7-x^2-x^4)/(1+xx^3-x^4), from=0, n=450000, add = TRUE) ... par(mar=c(5,0,4,2)) #No margin on the left side plot(c(0,1),type="n", axes=F, xlab="", ylab="") #Empty plot legend("top",legend=c("1","2","3","4","5"), density=c(20,20,20,20,NA), angle=c(90,0,45,135,NA), col=c(NA,NA,NA,NA,"black"), bty="n", cex=1.5) 

Regarding the point you want to mark, use the text (or mtext ) function to make it "programmaticaly" or the locator to do it interactively.

Change Alternatively (as I said in the comments), this will also work to put your legend outside the plot area and probably easier:

 par(mar=c(5,4,4,8)) plot(NA,xlim=c(0,1),ylim=c(0,1), xlab=expression(delta),ylab="K", xaxs="i",yaxs="i") # Here your delta a1 <- curve((x+x^7-x^2-x^4)/(1+xx^3-x^4), from=0, n=450000, add = TRUE) ... legend(1,1,legend=c("1","2","3","4","5"), density=c(20,20,20,20,NA), angle=c(90,0,45,135,NA), col=c(NA,NA,NA,NA,"black"), bty="n", cex=1.5, xpd=TRUE) 

Answer 2

The traditional graphics system provides the legend () function to add a legend or key to the plot. A legend is usually drawn inside the plot area, and is located relative to the user's coordinates.

The function has many arguments, here we need to use the angle and density arguments to distinguish between hashed areas.

 legend(0.5, 0.8, paste("region", 1:5), density=c(20,20,20,20,0), angle=c(90,0,45,135,0))