Como traçar uma distribuição weibull em r

Para traçar a função de densidade de probabilidade para uma distribuição Weibull em R, podemos usar as seguintes funções:

• dweibull(x, shape, scale = 1) para criar a função de densidade de probabilidade.
• curve(function, from = NULL, to = NULL) para traçar a função de densidade de probabilidade.

Para traçar a função de densidade de probabilidade, precisamos especificar o valor do parâmetro de forma e escala na função dweibull , bem como os valores de e para na função curve() .

Por exemplo, o código a seguir ilustra como traçar uma função de densidade de probabilidade para uma distribuição Weibull com parâmetros forma = 2 e escala = 1 onde o eixo x do gráfico vai de 0 a 4:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4)

Podemos adicionar um título, alterar o rótulo do eixo Y, aumentar a largura da linha e até mesmo alterar a cor da linha para tornar o gráfico mais esteticamente agradável:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, 
    main = 'Weibull Distribution (shape = 2, scale = 1)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

Também podemos adicionar várias curvas ao gráfico para comparar distribuições Weibull com diferentes parâmetros de forma e escala:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

Podemos adicionar uma legenda ao gráfico usando a função legend() , que assume a seguinte sintaxe:

legenda(x, y=NULL, legenda, preenchimento, col, bg, lty, cex)

Ouro:

• x, y: as coordenadas xey usadas para posicionar a legenda
• legenda: o texto a ser colocado na legenda
• fill: preenche a cor dentro da legenda
• col: a lista de cores a serem usadas nas linhas dentro da legenda
• bg: a cor de fundo da legenda
• lty: estilo de linha
• cex: tamanho do texto na legenda

Em nosso exemplo usaremos a seguinte sintaxe para criar uma legenda:

 #create density plots
curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

#add legend
legend(2, .7, legend=c("shape=2, scale=1", "shape=1.5, scale=1"),
       col=c("red", "blue"), lty=1, cex=1.2)