Участка Логнормальной плотность вероятности в R

Question

Я пытаюсь создать сюжет для логарифмически нормальной плотности вероятности в R, с 3-мя различными способами входа и журнал стандартов отклонения. Я пробовал следующее, но мой график настолько уродлив и выглядит совсем неважно.

x<- seq(0,10,length = 100) 
a <- dlnorm(x, meanlog = 0, sdlog = 1, log = FALSE) 
b <- dlnorm(x, meanlog = 0, sdlog = 1.5, log = FALSE) 
g <- dlnorm(x, meanlog = 1.5, sdlog = 0.2, log = FALSE) 
plot(x,a, lty=5, col="blue", lwd=3) 
lines(x,b, lty=2, col = "red") 
lines(x,g, lty=4, col = "green") 

Я даже пытался добавить легенду в правом верхнем углу для каждого среднелогарифмическое и стандартное отклонение журнала, но это не будет работать со мной. Мне было интересно, может ли кто-нибудь помочь мне в этом.

Right top of the graph

Answer 1

Там действительно ничего плохого в вашем коде. Вы только что забыли:

• type = "l" in plot;
• установить хорошее ylim провести все линии.

Вот простое решение с matplot:

matplot(x, cbind(a,b,g), type = "l", ylab = "density", main = "log-normal", 
     col = 1:3, lty = 1:3) 

Чтобы добавить легенду, используйте

legend("topright", 
     legend = c("mu = 0, sd = 1", "mu = 0, sd = 1.5", "mu = 1.5, sd = 0.2"), 
     col = 1:3, 
     lty = 1:3) 

Вы также можете прочитать ?plotmath для добавления выражений. Попробуйте изменить аргумент legend выше:

legend = c(expression(ln(y) %~% N(0,1)), 
      expression(ln(y) %~% N(0,1.5)), 
      expression(ln(y) %~% N(1.5,0.2)))