Calcular la distribución beta en la programación R: funciones dbeta(), pbeta(), qbeta() y rbeta()

La distribución beta en lenguaje R se define como una propiedad que representa los posibles valores de probabilidad. Este artículo es una ilustración de las funciones dbeta, pbeta, qbeta y rbeta de Beta Distribution.

Función dbeta()

Se define como función de densidad beta y se utiliza para crear un valor de densidad beta correspondiente al vector de cuantiles.

Sintaxis:
dbeta(vec, forma1, forma2)

Parámetro:
vec: Vector a utilizar
forma1, forma2: densidad beta de los valores de entrada

Devuelve: valores de densidad beta para un vector de cuantiles

Ejemplo :

# R program to illustrate the use of 
# dbeta() function
  
# Creating a vector
x_beta <- seq(0, 1.5, by = 0.025 )  
  
# Apply beta function          
y_beta <- dbeta(x_beta, shape1 = 2, shape2 = 4.5)  
  
# Plot beta values
plot(y_beta)                                       

Producción:

Función pbeta()

Se utiliza para crear una función de distribución acumulativa de la distribución beta.

Sintaxis:
pbeta(vec, forma1, forma2)

Parámetro:
vec: Vector a utilizar
forma1, forma2: densidad beta de los valores de entrada

Ejemplo:

# Specify x-values for pbeta function
x_pbeta <- seq(0, 1, by = 0.025)      
  
# Apply pbeta() function
y_pbeta <- pbeta(x_pbeta, shape1 = 1, shape2 = 4)  
  
# Plot pbeta values
plot(y_pbeta)

Producción:

Función qbeta()

Se conoce como función cuantil beta y se utiliza para devolver valores cuantiles de la función.

Sintaxis:
qbeta(vec, forma1, forma2)

Parámetros:
vec: Vector a utilizar
forma1, forma2: densidad beta de los valores de entrada

Ejemplo:

    
# Specify x-values for qbeta() function
x_qbeta <- seq(0, 1, by = 0.025)
  
# Apply qbeta() function
y_qbeta <- qbeta(x_qbeta, shape1 = 1, shape2 = 4)  
  
# Plot qbeta() values
plot(y_qbeta) 

Producción:

Función rbeta()

Se define como un generador de números aleatorios que se utiliza para establecer la semilla y especificar el tamaño de la muestra.

Sintaxis:
rbeta(N, forma1, forma2 )

Parámetros:
vec: Vector a utilizar
forma1, forma2: densidad beta de los valores de entrada

Ejemplo:

# Set seed for reproducibility
set.seed(13579)
  
# Specify sample size
N <- 10000  
  
# Draw N beta distributed values
y_rbeta <- rbeta(N, shape1 = 1, shape2 = 5)   
y_rbeta
  
# Plot of randomly drawn beta density
plot(density(y_rbeta), 
     main = "beta Distribution in R")

Producción: