Traccia due grafici nella stessa trama in R

Vorrei tracciare y1 e y2 nella stessa trama.

x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x, 1, 1)
plot(x, y1, type = "l", col = "red")
plot(x, y2, type = "l", col = "green")

Ma quando lo faccio in questo modo, non vengono tracciati insieme nella stessa trama.

In Matlab si può fare hold on, ma qualcuno sa come farlo in R?

lines()o points()aggiungerà al grafico esistente, ma non creerà una nuova finestra. Quindi dovresti farlo

plot(x,y1,type="l",col="red")
lines(x,y2,col="green")

Puoi anche usare pare tracciare sullo stesso grafico ma assi diversi. Qualcosa come segue:

plot( x, y1, type="l", col="red" )
par(new=TRUE)
plot( x, y2, type="l", col="green" )

Se leggi in dettaglio su parin R, sarai in grado di generare grafici davvero interessanti. Un altro libro da guardare è la R Graphics di Paul Murrel.

Quando si costruiscono grafici multistrato, è necessario considerare il ggplotpacchetto. L'idea è di creare un oggetto grafico con un'estetica di base e migliorarlo in modo incrementale.

ggplotlo stile richiede il confezionamento dei dati data.frame.

# Data generation
x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x,1,1)
df <- data.frame(x,y1,y2)

Soluzione di base:

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

Qui + operatorè usato per aggiungere livelli extra all'oggetto base.

Con ggplotte hai accesso all'oggetto grafico in ogni fase della stampa. Supponiamo che la normale configurazione dettagliata possa assomigliare a questa:

g <- ggplot(df, aes(x))
g <- g + geom_line(aes(y=y1), colour="red")
g <- g + geom_line(aes(y=y2), colour="green")
g

gproduce la trama e puoi vederla in ogni fase (beh, dopo la creazione di almeno un livello). Ulteriori incantesimi della trama sono fatti anche con l'oggetto creato. Ad esempio, possiamo aggiungere etichette per gli assi:

g <- g + ylab("Y") + xlab("X")
g

La finale gè simile a:

AGGIORNAMENTO (2013-11-08):

Come sottolineato nei commenti, ggplotla filosofia suggerisce di utilizzare i dati in formato lungo. Puoi fare riferimento a questa risposta per vedere il codice corrispondente.

Penso che la risposta che stai cercando sia:

plot(first thing to plot)
plot(second thing to plot,add=TRUE)

Usa la matplotfunzione:

matplot(x, cbind(y1,y2),type="l",col=c("red","green"),lty=c(1,1))

usare questo se y1e y2vengono valutati negli stessi xpunti. Ridimensiona l'asse Y per adattarlo a qualsiasi sia più grande ( y1o y2), a differenza di alcune delle altre risposte qui che verranno tagliate y2se diventa più grande di y1(le soluzioni ggplot per lo più vanno bene con questo).

In alternativa, e se le due linee non hanno le stesse coordinate x, impostare i limiti degli assi sul primo grafico e aggiungere:

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")

Stupito che questa Q abbia 4 anni e nessuno ha menzionato matploto x/ylim...

tl; dr: vuoi usare curve(con add=TRUE) o lines.

Non sono d'accordo par(new=TRUE)perché ciò consentirà di stampare due segni di spunta e le etichette degli assi. Per esempio

L'output di plot(sin); par(new=T); plot( function(x) x**2 ).

Guarda come sono incasinate le etichette degli assi verticali! Poiché le gamme sono diverse, dovresti impostarle ylim=c(lowest point between the two functions, highest point between the two functions), il che è meno facile di quello che sto per mostrarti --- e molto meno facile se vuoi aggiungere non solo due curve, ma molte.

Ciò che mi ha sempre confuso riguardo alla trama è la differenza tra curvee lines. (Se non è possibile ricordare che questi sono i nomi dei due importanti comandi plottaggio, basta cantare esso.)

Ecco la grande differenza tra curvee lines.

curvetraccia una funzione, come curve(sin). linespiazzole punti con valori di xey, come: lines( x=0:10, y=sin(0:10) ).

Ed ecco una differenza minore: curvedeve essere chiamato add=TRUEper quello che stai cercando di fare, mentre linesgià presume che stai aggiungendo a una trama esistente.

Ecco il risultato della chiamata plot(0:2); curve(sin).

Dietro le quinte, dai un'occhiata methods(plot). E controlla body( plot.function )[[5]]. Quando chiami plot(sin)R capisci che sinè una funzione (non valori y) e usa il plot.functionmetodo, che finisce per chiamare curve. Quindi curvelo strumento è pensato per gestire le funzioni.

se vuoi dividere il diagramma in due colonne (2 diagrammi uno accanto all'altro), puoi farlo in questo modo:

par(mfrow=c(1,2))

plot(x)

plot(y) 

Link di riferimento

Come descritto da @redmode, è possibile tracciare le due linee nello stesso dispositivo grafico usando ggplot. In quella risposta i dati erano in un formato "ampio". Tuttavia, quando ggplotlo si utilizza è generalmente più conveniente mantenere i dati in un frame di dati in un formato "lungo". Quindi, usando diverse "variabili di raggruppamento" negli aesargomenti di thetics, le proprietà della linea, come il tipo di linea o il colore, varieranno in base alla variabile di raggruppamento e appariranno le legende corrispondenti.

In questo caso, possiamo usare l' colouraessetica, che abbina il colore delle linee ai diversi livelli di una variabile nel set di dati (qui: y1 vs y2). Ma prima dobbiamo fondere i dati dal formato largo a quello lungo, usando ad es. La funzione 'melt' dal reshape2pacchetto. Altri metodi per rimodellare i dati sono descritti qui: Rimodellare data.frame da formato largo a lungo .

library(ggplot2)
library(reshape2)

# original data in a 'wide' format
x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x, 1, 1)
df <- data.frame(x, y1, y2)

# melt the data to a long format
df2 <- melt(data = df, id.vars = "x")

# plot, using the aesthetics argument 'colour'
ggplot(data = df2, aes(x = x, y = value, colour = variable)) + geom_line()

Se stai usando la grafica di base (cioè non la grata / la griglia), puoi imitare la funzione di blocco di MATLAB usando le funzioni punti / linee / poligoni per aggiungere ulteriori dettagli ai tuoi grafici senza iniziare una nuova trama. Nel caso di un layout multiplot, è possibile utilizzare par(mfg=...)per scegliere a quale trama aggiungere elementi.

È possibile utilizzare i punti per l'overplot, ovvero.

plot(x1, y1,col='red')

points(x2,y2,col='blue')

Invece di mantenere i valori da tracciare in un array, archiviarli in una matrice. Per impostazione predefinita, l'intera matrice verrà trattata come un set di dati. Tuttavia, se aggiungi lo stesso numero di modificatori al grafico, ad esempio col (), poiché hai le righe nella matrice, R scoprirà che ogni riga deve essere trattata in modo indipendente. Per esempio:

x = matrix( c(21,50,80,41), nrow=2 )
y = matrix( c(1,2,1,2), nrow=2 )
plot(x, y, col("red","blue")

Questo dovrebbe funzionare a meno che i tuoi set di dati non abbiano dimensioni diverse.

Idiomatic Matlab plot(x1,y1,x2,y2)può essere tradotto in R con ggplot2ad esempio in questo modo:

x1 <- seq(1,10,.2)
df1 <- data.frame(x=x1,y=log(x1),type="Log")
x2 <- seq(1,10)
df2 <- data.frame(x=x2,y=cumsum(1/x2),type="Harmonic")

df <- rbind(df1,df2)

library(ggplot2)
ggplot(df)+geom_line(aes(x,y,colour=type))

^{Ispirato dai grafici a doppia linea di Tingting Zhao con diversa gamma di assi x Usando ggplot2 .}

Potresti usare la ggplotly()funzione dal pacchetto plotly per trasformare uno qualsiasi degli esempi gggplot2 qui in un grafico interattivo, ma penso che questo tipo di trama sia migliore senza ggplot2 :

# call Plotly and enter username and key
library(plotly)
x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x, 1, 1)

plot_ly(x = x) %>%
  add_lines(y = y1, color = I("red"), name = "Red") %>%
  add_lines(y = y2, color = I("green"), name = "Green")

Puoi anche creare la tua trama usando ggvis :

library(ggvis)

x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x,1,1)
df <- data.frame(x, y1, y2)

df %>%
  ggvis(~x, ~y1, stroke := 'red') %>%
  layer_paths() %>%
  layer_paths(data = df, x = ~x, y = ~y2, stroke := 'blue')

Questo creerà il seguente diagramma:

Utilizzando plotly(aggiunta di soluzione plotlycon asse y primario e secondario- Sembra mancare):

library(plotly)     
x  <- seq(-2, 2, 0.05)
y1 <- pnorm(x)
y2 <- pnorm(x, 1, 1)

df=cbind.data.frame(x,y1,y2)

  plot_ly(df) %>%
    add_trace(x=~x,y=~y1,name = 'Line 1',type = 'scatter',mode = 'lines+markers',connectgaps = TRUE) %>%
    add_trace(x=~x,y=~y2,name = 'Line 2',type = 'scatter',mode = 'lines+markers',connectgaps = TRUE,yaxis = "y2") %>%
    layout(title = 'Title',
       xaxis = list(title = "X-axis title"),
       yaxis2 = list(side = 'right', overlaying = "y", title = 'secondary y axis', showgrid = FALSE, zeroline = FALSE))

Screenshot dalla demo funzionante:

possiamo anche usare la libreria reticolare

library(lattice)
x <- seq(-2,2,0.05)
y1 <- pnorm(x)
y2 <- pnorm(x,1,1)
xyplot(y1 + y2 ~ x, ylab = "y1 and y2", type = "l", auto.key = list(points = FALSE,lines = TRUE))

Per colori specifici

xyplot(y1 + y2 ~ x,ylab = "y1 and y2", type = "l", auto.key = list(points = F,lines = T), par.settings = list(superpose.line = list(col = c("red","green"))))