Két populáció varianciájának hipotézis tesztjét próbálom megalkotni. A tesztstatisztikám $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Ezt egy függvényként szeretném megvalósítani az R-ben, ezért kiszámoltam: a=function(B,A){ var(B)/var(A) }
Most a p-értéket szeretném kiszámítani az R tesztstatisztikájával. Hogyan tudom ezt megtenni?
Megjegyzések
Válasz
Mivel az F-statisztikát számolja, használhatja a pf függvényt.
Példaként vegyük a következő adatokat példaként.
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 Amint azt Gung megjegyezte, nem vettem észre, hogy Hartley-t próbálod megtenni Fmax teszt. Ebben az esetben nem is kell kiszámítania az F statisztikát, sőt a varianciákat sem, csak a következőt használhatja: var.test.
var.test(x, y)$p.value [1] 0.0001897506 Amint azt Gung is megjegyezte, van egy Levene teszt (leveneTest) a car csomagban.
carcsomagban.