Próbuję skonstruować test hipotezy dla wariancji dwóch populacji. Moja statystyka testowa to $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Chcę zaimplementować to jako funkcję w R, więc obliczyłem: a=function(B,A){ var(B)/var(A) }
Teraz chcę obliczyć wartość p przy użyciu tej statystyki testowej w R. Jak mogę to zrobić?
Komentarze
Odpowiedź
Ponieważ obliczasz statystykę F, możesz użyć funkcji pf.
Jako przykład weźmy następujące dane jako przykład.
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 zauważył Gung, nie zdawałem sobie sprawy, że próbujesz zrobić Hartleya Test Fmax. W takim przypadku nie musisz nawet obliczać statystyki F ani nawet wariancji, możesz po prostu użyć var.test.
var.test(x, y)$p.value [1] 0.0001897506 Jak zauważył Gung, istnieje test Levenea (leveneTest) w pakiecie car.
car.