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)

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.,

Livre de vírus. www.avast.com.

_______________________________________________
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.