glm Binomial com somas de sucessos e fracassos diferentes

Olá pessoal, Quero fazer um modelo que descreve a probabilidade da presença de uma espécie de peixe em um lago. Tenho algumas caraterísticas dos lagos (variáveis independentes) e relatos (entrevistas) da presença, ou não, do peixe em cada lago (variável dependente). O problema é que não tenho um número fixo de relatos por lago. Ou seja, para alguns lagos, dez pessoas sabiam se tinha a espécie, ou não. Mas para outros lagos somente duas pessoas sabiam. Assim, no modelo glm (family=binomial) entro com o número de sucessos (presença da espécies) e número de fracassos (ausência da espécie), mas a soma de sucessos e fracassos não é fixa. Isto é um problema para a análise? Existe como eu especificar pesos a priori? Muito Obrigado. #Dados# com = c(0,5,3,6,1,9,0,0,4,1,0,0,4,0,2,0,0,0,2) sem = c(9,4,8,3,8,2,9,9,8,10,1,2,0,2,0,2,4,2,0) n.tot = c(340.52,344.07,147.43,411.64,264.75,263.24,1087.00,1028.32,697.57, 428.56,414.46,538.18,519.69,590.99,472.54,381.82,465.55,528.43,334.95) D.L70 = c(4.44,3.92,6.02,3.80,3.23,4.56,8.53,9.56,2.29,3.61,14.99,5.85, 3.96,15.30,3.96,1.01,1.32,1.79,4.51) GLM = glm(cbind(com,sem)~n.tot+D.L70,family = binomial) summary(GLM) Diego PJ

Não há necessidade do size dessas binomiais serem o mesmo. Pode passar pesos com glm(..., weights=). Tome cuidado com superdispersão nessa sua análise. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Diego Pujoni
-
Walmes Zeviani