Estou tentando construir um teste de hipótese para a variância de duas populações. Minha estatística de teste é $ F = $$ {\ sigma_b ^ 2} $ / $ {\ sigma_A ^ 2} $. Quero implementar isso como uma função em R, então calculei: a=function(B,A){ var(B)/var(A) }
Agora, quero calcular o valor p usando essa estatística de teste em R. Como posso fazer isso?
Comentários
Resposta
Já que você está calculando a estatística F, pode usar a função pf.
Como exemplo, vejamos os seguintes dados como exemplo.
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 Como Gung observou, não percebi que você está tentando fazer Hartley “s Teste Fmax. Nesse caso, você nem mesmo precisa calcular a estatística F ou mesmo as variâncias, você pode apenas usar var.test.
var.test(x, y)$p.value [1] 0.0001897506 Também conforme observado por gung, há o teste de Levene (leveneTest) no pacote car.
car.