Ha duas formas, dados agrupados por tratamentos ou os dados brutos, onde você coloca a resposta (Bernoulli) em cada repetição. Veja a forma de dados agrupados: trat <- factor(letters[1:5]) n <- c(10, 10, 10, 10, 10) ## 10 repeticoes em cada tratamento y <- c(5, 7, 3, 1, 2) result <- glm(cbind(y, n-y)~trat, family=binomial) summary(result) anova(result, test='Ch') ### Caso haja algum y=0 ou y=n, uma opção é usar library(brglm) res2 <- brglm(cbind(y, n-y)~trat, family=binomial) summary(res2) anova(res2) Elias On 28/08/2018 14:31, ana paula coelho madeira via R-br wrote:
Prezados colegas,
Boa tarde! Estou avaliando um experimento com 5 tratamentos e a varável resposta é binária ( conta a presença ou ausência de uma dada característica). Sei (espero não estar enganada) que se trata de uma regressão logística, visto que a variável reposta é binária e a explicativa categórica. Não estou conseguindo tratar esses dados no R. Encontrei um script no livro do Crawley porém não consegui adaptar. Peço indicações. Desde já agradeço
Att.,
<https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Livre de vírus. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.