# script_gva_wei
# Autor
# 08/04/2020
# Função para gerar valores com distribuição weibull 
#  com parâmetro lambda e alfa
rm(list = ls())

# Semente
set.seed(23)

# gerador de variável aleatória
# tamanho da amostra
n <- 100

# valores dos parâmetros
alfa <- 2
lambda <- 1

# gerador
y <- rep(0,n)
for(i in 1:n){
  u <- runif(1)
  y[i] <- lambda*((-log(1-u))^(1/alfa))
}

#função
wei <- function(alfa,lambda,n){
  y <- rep(0,n)
  for(i in 1:n){
    u <- runif(1)
    y[i] <- lambda*((-log(1-u))^(1/alfa))
  }
  return(y)
}

#valores da distribuição Weibull
x <- wei(2,1,100)

# Para verificar se a amostra gerada segue a distribuição 
# Função de distribuição empírica
plot(ecdf(x))
# Função de distribuição Weibull 
curve(pweibull(x,2,1), add = TRUE, type = 'l', col = 'red')

# ou
# Histograma
hist(x, freq = FALSE, main = '', xlab = 'Valores de y',
     ylab = 'Densidade')
# Função densidade Weibull
curve(dweibull(x,2,1), add = TRUE, type = 'l', 
      col = 'red')
# Legenda
legend(1.5,0.75,c('Weibull'),col='red',lty=1, bty = 'n')