Prueba F

A continuación se comparte el código en R para la prueba F:

#Cargar los datos

# x=datos con mayor varianza

#Información de cada muestra
var(x)
nx=length(x)
var(y)
ny=length(y)

#Elaboración de la prueba F - manual
alfa<-0.05

#F tabulado
Ftab=qf((1-(alfa/2)),nx-1,ny-1)
Ftab

#F calculado
Fcal=var(x)/var(y)           
Fcal

#p-value
p_value<-2*(pf(var(y)/var(x),nx-1,ny-1))
p_value

#IC para la varianza
IC=c((1/qf((1-(alfa/2)),ny-1,nx-1))*(var(x)/var(y)),(1/qf((alfa/2),ny-1,nx-1))*(var(x)/var(y)))
IC

#Prueba F - codigo de R
var.test(x,y, conf.level = 1-alfa)