Yritän rakentaa hypoteesitestin kahden populaation varianssille. Testitilastoni on $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Haluan toteuttaa tämän funktiona R: ssä, joten laskin: a=function(B,A){ var(B)/var(A) }
Haluan nyt laskea p-arvon käyttämällä tätä testitilastoa R: ssä. Kuinka voin tehdä tämän?
Kommentit
Vastaa
Koska lasket F-tilastoa, voit käyttää pf -toimintoa.
Otetaan esimerkkinä seuraavat tiedot esimerkkeinä.
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 Kuten gung totesi, en tajunnut, että yrität tehdä Hartleylle Fmax-testi. Tällöin sinun ei tarvitse edes laskea F-tilastoa tai edes variansseja, voit vain käyttää var.test.
var.test(x, y)$p.value [1] 0.0001897506 Kuten myös gung huomautti, car -paketissa on Levene-testi (leveneTest).
car.