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