[R-br] glm Binomial com somas de sucessos e fracassos diferentes

[Diego Pujoni]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120702/85e8b2b9/attachment.html>