Jessaie de construire un test dhypothèse pour la variance de deux populations. Ma statistique de test est $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Je veux implémenter cela en tant que fonction dans R donc jai calculé: a=function(B,A){ var(B)/var(A) }
Maintenant, je veux calculer la valeur p en utilisant cette statistique de test dans R. Comment puis-je faire cela?
Commentaires
Réponse
Puisque vous calculez la statistique F, vous pouvez utiliser la fonction pf.
À titre dexemple, prenons les données suivantes comme exemple.
set.seed(123) x <- rnorm(50, mean = 0, sd = 2) y <- rnorm(30, mean = 1, sd = 1) f <- var(x)/var(y) # make sure to double for two-sided pf(f, df1=49, df2=29, lower.tail=F)*2 [1] 0.0001897506 Comme la noté Gung, je nai pas réalisé que vous essayez de faire Hartley « s Test Fmax. Dans ce cas, vous navez même pas besoin de calculer la statistique F ni même les variances, vous pouvez simplement utiliser var.test.
var.test(x, y)$p.value [1] 0.0001897506 Également comme indiqué par gung, il y a le test de Levene (leveneTest) dans le paquet car.
car.