Boa tarde pessoal, desejo fazer um teste de hipoteses bem simples, porem nao consigo visualizar como deveria fazer isso no R, uma vez que os testes que ja fiz sempre foram utilizando um modelo (lm) e a funcao linearHypothesis. Nesse caso tenho os dados da seguinte forma:
cov V1 V2 V3 V4 V5 V6 45 65 81 55 31 0 41 80 73 51 36 0 42 79 74 40 45 0 43 77 76 62 * 0
Onde as colunas V1...V6 sao niveis de agua adicionados a uma planta, e cada linha eh uma caixa onde o experimento foi realizado, nesse caso 4 caixas para cada nivel de agua. Os valores sao o numero de semestes que germinaram. Quero testar a hipotese de que o nivel de agua nao influencia na germinacao das sementes. A principio pensei em usar anova, mas nao tenho claro como devo escrever meu modelo. Alguem poderia me dar uma dica? Abracos e bom dia! Daniela -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "Not everything that can be counted counts; and not everything that counts can be counted.” George Gallup
Daniela, Ao que parece é algo simples mesmo. Escreva os seus dados no formato longo, duas colunas (água, sementes), coloque num data.frame, delcare que os níveis de água são qualitativos (factor()) (caso queira ajustar um modelo de regressão, que por sinal é mais informativo no geral, use água como numérica), use a função lm() m0 <- lm(sementes~agua, data=dados) depois aplique contrastes (contrast::contrast(), gmodels::fit.contrast()) ou comparações múltiplas (stats::TukeyHSD(), agricole::HSD.test(), ScottKnott::SK(), ...). Caso o número de sementes inicial colocadas para germinar for conhecido, pode-se atribuir a distribuição binomial para a resposta "número de sementes germinadas para um número de submetidas a germinação". Então, use a função glm(..., family=binomial). Exemplos desse tipo de análise você encontra na página do Professor Paulo Justiniano (http://leg.ufpr.br/~paulojus/embrapa/Rembrapa/). À 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 (2)
-
Daniela Recchia -
Walmes Zeviani