Plotting два Пуассон процессов на одном участке

Question

У меня есть два пуассоновский процессы:

n <- 100 
x <- seq(0, 10, length = 1000) 
y1 <- cumsum(rpois(1000, 1/n)) 
y2 <- -cumsum(rpois(1000, 1/n)) 

Я хотел бы построить их на одном графике и ожидают, что y1 лежит выше оси х и y2 лежит ниже оси х. Я пробовал использовать следующий код:

plot(x, y1) 
par(new = TRUE) 
plot(x, y2, col = "red", 
axes = FALSE, 
xlab = '', ylab = '', 
xlim = c(0, 10), ylim = c(min(y2), max(y1))) 

но это не сработало. Может кто-нибудь, пожалуйста, скажите мне, как это исправить? (Я работаю с R для моего кода)

Большого спасибо заранее

Answer 1

Как насчет

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

Я хотел бы предложить, пытаясь избежать несколько вызовов участка с par(new=TRUE). Это обычно очень грязно. Здесь мы используем lines(), чтобы добавить к существующему сюжету. Единственный улов в том, что пределы x и y не будут меняться в зависимости от новых данных, поэтому мы используем ylim в первом вызове plot(), чтобы установить диапазон, соответствующий всем данным.

Answer 2

Или, если вы не хотите, чтобы беспокоиться о пределах (как MrFlick упоминается) или число строк, вы можете также поток вверх ваши данные и используя melt и ggplot

df <- data.frame(x, y1, y2) 
library(reshape2) 
library(ggplot2) 
mdf <- melt(df, "x") 
ggplot(mdf, aes(x, value, color = variable)) + 
    geom_line()