Não tem como usar a média das duas observações e ao mesmo tempo usar a distribuição de Poisson. Você pdoeria tentar descobrir qual é a distribuição de probabilidade da média de duas variáveis aleatórias Poisson i.i.d, mas acho que essa não é a melhor saída.
Para mim, a melhor saída seria usar um modelo de medidas repetidas, e considerar um efeito aleatório de indivíduo para cada dia.
Para isso, considerando que você tem um banco de dados da seguinte forma:
> dados <- data.frame(
+ id = 1:100,
+ explicativa = runif(100, 0, 20)
+ )
>
> dados$r1 <- rpois(100, dados$explicativa)
> dados$r2 <- rpois(100, dados$explicativa)
>
> head(dados)
id explicativa r1 r2
1 1 9.082513 16 14
2 2 17.741123 14 29
3 3 10.819865 13 12
4 4 18.527938 22 25
5 5 4.828392 6 7
6 6 13.986794 14 15
r1
e r2
são as frequências observadas no dia 1 e no dia 2.
Transforme seus dados no formato tidy/long
:
library(tidyr);library(dplyr)
dados <- dados %>% gather(dia, resposta, starts_with("r"))
Depois ajustaria um modelo da seguinte forma:
library(lme4)
modelo <- glmer(resposta ~ explicativa + (0 + dia | id), data = dados, family = poisson)
Assim, você estará considerando que para cada indivíduo existe uma variação aleatória relacionada ao dia em que ele está sendo medido.