Încerc să construiesc un test de ipoteză pentru varianța a două populații. Statistica mea de testare este $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Vreau să implementez acest lucru ca o funcție în R, așa că am calculat: a=function(B,A){ var(B)/var(A) }
Acum, vreau să calculez valoarea p folosind această statistică de test în R. Cum pot face acest lucru?
Comentarii
Răspuns
Deoarece calculați statistica F, puteți utiliza funcția pf.
Ca exemplu, să luăm următoarele date ca exemplu.
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 După cum a menționat gung, nu am reușit să realizez că încercați să faceți Hartley Test Fmax. În acest caz, nici măcar nu trebuie să calculați statistica F sau chiar varianțele, puteți folosi doar var.test.
var.test(x, y)$p.value [1] 0.0001897506 De asemenea, așa cum este menționat de gung, există testul lui Levene (leveneTest) în pachetul car.
car.