DISTRIBUCIÓN F
Suponiendo que damos los valores de equivalencia a pt= nu, siendo este 4 y pf, tenemos en cuanta una secuancia de x que va desde 0.001 al 5 con longitud de 100. Aplicar la distribución de F
x =seq(0.001, 5, len = 100)
nu = 4
stopifnot(all.equal(2*pt(x,nu) - 1, pf(x^2, 1,nu)),all.equal(2*pt(x,nu, lower=FALSE),pf(x^2, 1,nu, lower=FALSE)))
#aplicar el upper.tails
all.equal(df(x^2, 1, 5), dt(x, 5)/x)
## [1] TRUE
Evaluar la secuencia para que el resultado sea apróximado a 7e-9
p = seq(1/2, .99, length = 50); df <- 10
rel.err = function(x, y) ifelse(x == y, 0, abs(x-y)/mean(abs(c(x,y))))
quantile(rel.err(qf(2*p - 1, df1 = 1, df2 = df), qt(p, df)^2), .90)
## 90%
## 5.159512e-16
2) Encontrar el percentil 95 con grados de libertad entre (5,2)
qf(.95, df1=5, df2=2)
## [1] 19.29641
Conclusión
La distribución de F es usada para casos especificos donde las graficas tengan grados de liberación dadas por el mismo problema, en este caso combinamos dos ejercicios en uno para practicar los comandos.