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
[email protected]
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.
_______________________________________________
R-br mailing list
[email protected]
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.
