Jeg prøver å konstruere en hypotesetest for variansen til to populasjoner. Teststatistikken min er $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Jeg vil implementere dette som en funksjon i R, så jeg beregnet: a=function(B,A){ var(B)/var(A) }
Nå vil jeg beregne p-verdien ved hjelp av denne teststatistikken i R. Hvordan kan jeg gjøre dette?
Kommentarer
Svar
Siden du beregner, kan F-statistikken bruke pf -funksjonen.
La oss ta et eksempel på følgende 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 Som gung bemerket, klarte jeg ikke å innse at du prøver å gjøre Hartleys Fmax-test. I så fall trenger du ikke engang å beregne F-statistikken eller til og med avvikene, du kan bare bruke var.test.
var.test(x, y)$p.value [1] 0.0001897506 Også som bemerket av gung, er det Levene test (leveneTest) i car -pakken.
car-pakken.