Para Poisson com função de ligação log o modelo fica exp(X%*%beta) que no caso fica
* tratamento base: exp(b0+b1*trat)
* tratamento outro: exp(b0+b0i+b1*trat);
Para porção zero que considera Bernoulli com função de ligação logit fica exp(X%*%beta)/(1+exp(X%*%beta)), assim
* tratamento base: exp(b0)/(1+exp(b0))
* tratamento outro: exp(b0+b0i)/(1+exp(b0+b0i));
As contas matriciais que fiz (veja onde tem o objeto eta) estão correspondendo à essas especificações que foram os modelos que eu declarei ao gerar os dados artificiais.
À 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
skype: walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
_______________________________________________ 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.
--
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr@yahoo.com.br
alexandre.santos@cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================