Как построить функцию плотности вероятности в R?

Question

Я пытаюсь построить функцию плотности вероятности с ограничением x от -1 до 0 и от 0 до 1, поэтому я делаю два графика:

x1 = seq(-1, 0, 0.01)
x2 = seq(0, 1, 0.01)

eq1 = function(x) {(1+x)^2}
eq2 = function(x) {(1+x)^3}

plot(x1, eq1, col="red")
par(new = TRUE)
plot(x2, eq2, type = "l", col = "green")

Однако я получаю следующую ошибку:

Ошибка в xy.coords (x, y, xlabel, ylabel, log): длины 'x' и 'y' отличается.

Я не уверен, что случилось.

Answer 1

Как отмечалось в комментариях, вторым аргументом plot() (т.е. y) должен быть вектор:

x1 = seq(-1, 0, 0.01)
x2 = seq(0, 1, 0.01)

eq1 = function(x) {(1+x)^2}
eq2 = function(x) {(1+x)^3}

plot(x1, eq1(x1), col="red")
par(new = TRUE)
plot(x2, eq2(x2), type = "l", col = "green")