Ajuste os modelos individualmente: diurno=c(68,36,40,30) noturno=c(46,26,55,22) quadra= c("A", "B", "G","D") dados=data.frame(diurno,noturno,quadra) fit1 = glm(diurno ~ factor(quadra), family = poisson, data=dados) fit2 = glm(noturno ~ factor(quadra), family = poisson, data=dados) Isto deve resolver o teu problema. 2016-06-07 20:26 GMT-03:00 Andre Oliveira <andreolsouza@yahoo.com.br>:
boa noite, alguém poderia ajudar descobrir onde está o erro?
remove(list=ls()) diurno=c(68,36,40,30) noturno=c(46,26,55,22) quadra= c("A", "B", "G","D") dados=data.frame(diurno,noturno,quadra) fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, data=dados)
*Error in x[good, , drop = FALSE] : * * (subscript) subscrito lógico muito longo*
obrigado
André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
_______________________________________________ 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.
-- Marcus Nunes http://marcusnunes.me/