lines()
or points()
will add to the existing graph, but will not create a new window. So you'd need to do
plot(x,y1,type="l",col="red") lines(x,y2,col="green")
ID : 3026
viewed : 139
96
lines()
or points()
will add to the existing graph, but will not create a new window. So you'd need to do
plot(x,y1,type="l",col="red") lines(x,y2,col="green")
89
You can also use par
and plot on the same graph but different axis. Something as follows:
plot( x, y1, type="l", col="red" ) par(new=TRUE) plot( x, y2, type="l", col="green" )
If you read in detail about par
in R
, you will be able to generate really interesting graphs. Another book to look at is Paul Murrel's R Graphics.
80
When constructing multilayer plots one should consider ggplot
package. The idea is to create a graphical object with basic aesthetics and enhance it incrementally.
ggplot
style requires data to be packed in data.frame
.
# Data generation x <- seq(-2, 2, 0.05) y1 <- pnorm(x) y2 <- pnorm(x,1,1) df <- data.frame(x,y1,y2)
Basic solution:
require(ggplot2) ggplot(df, aes(x)) + # basic graphical object geom_line(aes(y=y1), colour="red") + # first layer geom_line(aes(y=y2), colour="green") # second layer
Here + operator
is used to add extra layers to basic object.
With ggplot
you have access to graphical object on every stage of plotting. Say, usual step-by-step setup can look like this:
g <- ggplot(df, aes(x)) g <- g + geom_line(aes(y=y1), colour="red") g <- g + geom_line(aes(y=y2), colour="green") g
g
produces the plot, and you can see it at every stage (well, after creation of at least one layer). Further enchantments of the plot are also made with created object. For example, we can add labels for axises:
g <- g + ylab("Y") + xlab("X") g
Final g
looks like:
UPDATE (2013-11-08):
As pointed out in comments, ggplot
's philosophy suggests using data in long format. You can refer to this answer in order to see the corresponding code.
69
I think that the answer you are looking for is:
plot(first thing to plot) plot(second thing to plot,add=TRUE)
59
Use the matplot
function:
matplot(x, cbind(y1,y2),type="l",col=c("red","green"),lty=c(1,1))
use this if y1
and y2
are evaluated at the same x
points. It scales the Y-axis to fit whichever is bigger (y1
or y2
), unlike some of the other answers here that will clip y2
if it gets bigger than y1
(ggplot solutions mostly are okay with this).
Alternatively, and if the two lines don't have the same x-coordinates, set the axis limits on the first plot and add:
x1 <- seq(-2, 2, 0.05) x2 <- seq(-3, 3, 0.05) y1 <- pnorm(x1) y2 <- pnorm(x2,1,1) plot(x1,y1,ylim=range(c(y1,y2)),xlim=range(c(x1,x2)), type="l",col="red") lines(x2,y2,col="green")
Am astonished this Q is 4 years old and nobody has mentioned matplot
or x/ylim
...