Pokouším se sestavit test hypotézy pro rozptyl dvou populací. Moje statistika testu je $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Chci to implementovat jako funkci v R, takže jsem vypočítal: a=function(B,A){ var(B)/var(A) }
Teď chci vypočítat p-hodnotu pomocí této testovací statistiky v R. Jak to mohu udělat?
Komentáře
Odpověď
Protože vypočítáváte F-statistiku, můžete použít funkci pf.
Jako příklad si vezměme jako příklad následující data.
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 Jak poznamenal gung, neuvědomil jsem si, že se snažíte udělat Hartleyho Test Fmax. V takovém případě nepotřebujete ani vypočítat statistiku F, ani odchylky, stačí použít var.test.
var.test(x, y)$p.value [1] 0.0001897506 Jak poznamenal gung, v balíčku car je také test Levene (leveneTest).
car.