<!DOCTYPE html>
<html>
<head>
<title></title>
</head>
<body><div>Em Ter 7 jun. 2016, às 20:26, Andre Oliveira escreveu:<br></div>
<blockquote type="cite"><div style="color:rgb(0, 0, 0);background-color:rgb(255, 255, 255);font-family:bookman old style,new york,times,serif;font-size:16px;"><div dir="ltr">boa noite, alguém poderia ajudar descobrir onde está o erro? <br></div>
<div dir="ltr"> </div>
<div>remove(list=ls())<br></div>
<div style="font-family:Arial;">diurno=c(68,36,40,30)<br></div>
<div style="font-family:Arial;">noturno=c(46,26,55,22)<br></div>
<div style="font-family:Arial;">quadra= c("A", "B", "G","D")<br></div>
<div style="font-family:Arial;">dados=data.frame(diurno,noturno,quadra)<br></div>
<div style="font-family:Arial;">fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, data=dados)<br></div>
<div style="font-family:Arial;"><b>Error in x[good, , drop = FALSE] : <br></b></div>
<div dir="ltr"><b> (subscript) subscrito lógico muito longo</b><br></div>
</div>
</blockquote><div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Andre, conforme está escrito na documentação do glm, a resposta pode ser fornecida como duas colunas (sucessos e fracassos) <i>apenas para regressões binomiais e quasi-binomiais</i>; não para Poisson.<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Supondo que seus dados sejam o resumo de 68 + 46 = 114 observações para a quadra A, 36 + 26 = 62 observações para a quadra B, e por aí em diante; e que você queira calcular uma razão de prevalência / risco / outra proporção entre as quadras, o jeito mais óbvio de fornecer seus dados para a função é em formato longo:<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">dados <- data.frame(horariodiurno = c(1, 0, 1, 0, 1, 0, 1, 0), quadra = c("A", "A", "B", "B", "C", "C", "D", "D"))<br></div>
<div style="font-family:Arial;">dados <- dados[rep(1:nrow(dados), c(68, 46, 36, 26, 40, 55, 30, 22)), ]<br></div>
<div style="font-family:Arial;">fit = glm(horariodiurno ~ quadra, family = quasipoisson, data = dados)<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Repare também que, se você quer trabalhar com razão de proporções, a distribuição Poisson em princípio não se aplica, daí a necessidade de você especificar uma quasipoisson ou então um estimador de variância consistente com heterocedasticidade (pacote sandwich).<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Por favor, use o código com cuidado porque não estou com meu computador aqui para testá-lo.<br></div>
<div style="font-family:Arial;"> </div>
<div style="font-family:Arial;">Atenciosamente,<br></div>
<div style="font-family:Arial;"> </div>
<div id="sig4487139"><div id="signature" class="signature"><a title="Currículo Lattes" href="http://lattes.cnpq.br/9234772336296638">Leonardo Ferreira Fontenelle</a><br></div>
</div>
</body>
</html>