[R-br] Resposta binária
Elias T Krainski
eliaskrainski em yahoo.com.br
Terça Agosto 28 14:44:45 -03 2018
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 em 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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180828/8b5b95a3/attachment.html>
Mais detalhes sobre a lista de discussão R-br