ajuda teste hipótese duas variáveis

OI PESSOAL sou uma aprendiz bem iniciante do R, preciso de ajuda com um teste de hipótese, alguém pode me ajudar? bom, é o seguinte: tenho um conj de dados com duas variáveis, quero fazer um t.test para verificar se as médias dos salarios (variável numerica) de funcionários com diferente grau de instrução (2ºgrau e superior, variável categorica) são diferentes ou não. tentei fazer o seguinte: meu arquivo é da apostila do prof. justiniano o velho arquivo 'milsa'.(bem de iniciante) # teste de hipótese: H0:salario 2grau = superior# formulei essa hipotese # H1:salario 2grau diferente superior milsa$instrucao <- factor(milsa$instrucao, label = c("2oGrau", "Superior"), lev = 1:2, ord = T) var.test(milsa$salario~milsa$instrucao) alpha <- 0.05 # especifica valor da signific^ancia if (var.test$p.value > alpha) diff.test <- t.test(milsa$salario ~ milsa$instrucao,alternative = "two.sided", var.equal = TRUE, conf.level = 1-alpha) else diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative = "two.sided", var.equal = FALSE,conf.level = 1-alpha) diff.test não deu certo, o que estou esquecendo? não acho o erro. muito obrigada... Atenciosamente, Carla Eloize Carducci Engº Agronoma (UFGD) Msc. Ciência do Solo (UFLA) Doutoranda em Ciência do Solo (UFLA)

Não deu certo como? O que saiu no R? A variável salário esta adequadamente reconhecida como numérica? 2011/5/24 carla eloize carducci <elocarducci@hotmail.com>
OI PESSOAL sou uma aprendiz bem iniciante do R, preciso de ajuda com um teste de hipótese, alguém pode me ajudar? bom, é o seguinte: tenho um conj de dados com duas variáveis, quero fazer um t.test para verificar se as médias dos salarios (variável numerica) de funcionários com diferente grau de instrução (2ºgrau e superior, variável categorica) são diferentes ou não. tentei fazer o seguinte: meu arquivo é da apostila do prof. justiniano o velho arquivo 'milsa'.(bem de iniciante) # teste de hipótese: H0:salario 2grau = superior# formulei essa hipotese # H1:salario 2grau diferente superior
milsa$instrucao <- factor(milsa$instrucao, label = c("2oGrau", "Superior"), lev = 1:2, ord = T) var.test(milsa$salario~milsa$instrucao) alpha <- 0.05 # especifica valor da signific^ancia if (var.test$p.value > alpha) diff.test <- t.test(milsa$salario ~ milsa$instrucao,alternative = "two.sided", var.equal = TRUE, conf.level = 1-alpha) else diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative = "two.sided", var.equal = FALSE,conf.level = 1-alpha) diff.test não deu certo, o que estou esquecendo? não acho o erro.
muito obrigada...
*Atenciosamente*,
*Carla Eloize Carducci *Engº Agronoma (UFGD) Msc. Ciência do Solo (UFLA) Doutoranda em Ciência do Solo (UFLA)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

Carla, Você esqueceu de atribuir a algum objeto o resultado do var.test() que usa na pergunta lógica. Contruí uma função (super simples por sinal) na disciplina de Estatística Computacional que faz justamente o procedimento que você apresentou. Confira em http://www.leg.ufpr.br/doku.php/disciplinas:ce223-2011-01 o script da aula15.R a função auto.t.test(). # meus dados simulados milsa milsa <- data.frame(instrucao=gl(2,14)) milsa$salario <- rnorm(nrow(milsa)) str(milsa) var.test <- var.test(milsa$salario~milsa$instrucao) # aqui ó <---------------- alpha <- 0.05 # especifica valor da signific^ancia if(var.test$p.value > alpha){ diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative="two.sided", var.equal=TRUE, conf.level=1-alpha) } else { diff.test <- t.test(milsa$salario ~ milsa$instrucao, alternative="two.sided", var.equal=FALSE, conf.level=1-alpha) } diff.test À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (3)
-
carla eloize carducci
-
Daniel C Bezerra
-
Walmes Zeviani