Distribución t inversa en R

En probabilidad y estadística, la distribución t de Student (o simplemente la distribución t) es cualquier miembro de una familia de distribuciones de probabilidad continuas. Viene a la imagen al estimar la media de una población distribuida normalmente, donde el tamaño de la muestra es pequeño y la desviación estándar (SD) es indefinida y desconocida. La distribución t con gl = n grados de libertad tiene densidad

f(x) = \gamma(\frac{n+1}{2}) / (\sqrt(n pi) \gamma(n/2)) (1 + x^2/n)^-((n+1)/2)

Si una variable aleatoria X tiene una distribución T con ν grados de libertad, entonces Pr (X ≤ x) = P. Suponga que Z e Y son variables aleatorias independientes. Sea Z una variable aleatoria normal estándar (media de 0, la varianza de 1). Y Y una variable chi-cuadrado con m grados de libertad. Ahora tenemos, 

T = \frac{Z}{\sqrt{\frac{Y}{m}}}

La variable T tiene la distribución t con m grados de libertad.

Método 1: Usar el método tinv()

PEIP es un paquete en R que se usa para aplicar muchas traducciones y cálculos de MATLAB en R. Es útil para estadísticas de datos y operaciones matemáticas. Este método es más preciso en comparación con el método qt(). El paquete se puede instalar en el espacio de trabajo usando el siguiente comando: 

install.packages("PEIP")

El método tinv() de este paquete se basa en la función betaincinv. Se utiliza para calcular la distribución t inversa. Devuelve el quintil para la distribución T. Un valor p es una probabilidad de que los resultados de sus datos de muestra se hayan producido por casualidad.

Sintaxis: tinv(p, nu)

Argumentos: 

  • p : valor p
  • nu : grados de libertad

El siguiente fragmento de código ilustra el cálculo de la distribución t inversa considerando la probabilidad aleatoria de un número 0.4 (valor p) con 2 grados de libertad: 

R

# loading the required library
library("PEIP")
  
# defining p-value
p <- 0.4
nu <- 2
print ("Inverse t-distribution")
tinv(p, nu)

Producción:

[1] "Inverse t-distribution" 
[1] -0.2887

Método 2: Usar el método qt()

El método qt() en base R se utiliza para devolver la densidad acumulada de probabilidad inversa de la distribución t de Student, también conocida como distribución T. Es básicamente una función cuantil. 

Sintaxis: qt(p, nu )

Argumentos:

  • p : vector de probabilidades.
  • nu : grados de libertad

El siguiente fragmento de código ilustra el cálculo de la distribución t inversa considerando la probabilidad aleatoria de un pequeño número exponencial 1e-20 (valor p) con 2 grados de libertad. El signo negativo ilustra que está declinando en la naturaleza mucho hacia el -ve NaN: 

R

# defining p-value
p <- 1e-20
nu <- 2
  
print ("Inverse t-distribution")
qt(p, nu)

Producción:

[1] "Inverse t-distribution" 
[1] -7071067812

Se toma un vector de valores x en forma de secuencia, considerando que el número de grados de libertad es 15. Los valores y correspondientes se pueden calcular usando el método qt() en R. El siguiente fragmento de código ilustra el uso de este método sobre un vector de valores: 

R

# creating a vec
x <- seq(.1, .9, by = .1)
  
# Calling qt() Function
y <- qt(x, 15)
print ("inverse t-distribution")
print (y)
  
# plotting data
plot (x, y, type = "o")

Producción:

[1] “distribución t inversa” 

[1] -1,3406056 -0,8662450 -0,5357291 -0,2578853 0,0000000 0,2578853 0,5357291 0,8662450 1,3406056

Publicación traducida automáticamente

Artículo escrito por yippeee25 y traducido por Barcelona Geeks. The original can be accessed here. Licence: CCBY-SA

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *