Gere amostras de dados da regressão de Poisson

13

Eu queria saber como você geraria dados de uma equação de regressão de Poisson em R? Estou meio confuso sobre como abordar o problema.

Portanto, se eu assumir que temos dois preditores e que são distribuídos . E a interceptação é 0 e ambos os coeficientes são iguais a 1. Então, minha estimativa é simplesmente: $X_1$ $X_2$ $N(0,1)$

registro (Y) = 0 0 + 1 \cdot X_{1} + 1 \cdot X_{2}

$\log(Y) = 0+ 1\cdot X_1 + 1\cdot X_2$

Mas uma vez calculado o log (Y) - como faço para gerar contagens de poisson com base nisso? Qual é o parâmetro rate para a distribuição Poisson?

Se alguém pudesse escrever um breve script R que gere amostras de regressão de Poisson, isso seria incrível!

r regression poisson-distribution simulation

— Michael
fonte

24

$Y$ $\log$ $x$ $Y \sim P(\mu)$ $E(Y) = V(Y) = \mu$ $\log(\mu) = \beta_0 + \beta_1 x$

>   #sample size
> n <- 10
>   #regression coefficients
> beta0 <- 1
> beta1 <- 0.2
>   #generate covariate values
> x <- runif(n=n, min=0, max=1.5)
>   #compute mu's
> mu <- exp(beta0 + beta1 * x)
>   #generate Y-values
> y <- rpois(n=n, lambda=mu)
>   #data set
> data <- data.frame(y=y, x=x)
> data
   y         x
1  4 1.2575652
2  3 0.9213477
3  3 0.8093336
4  4 0.6234518
5  4 0.8801471
6  8 1.2961688
7  2 0.1676094
8  2 1.1278965
9  1 1.1642033
10 4 0.2830910

— ocram
fonte

3

Se você deseja gerar um conjunto de dados que se encaixe perfeitamente no modelo, você pode fazer algo assim em R:

# y <- exp(B0 + B1 * x1 + B2 * x2)

set.seed(1234)

B0 <-  1.2                # intercept
B1 <-  1.5                # slope for x1
B2 <- -0.5                # slope for x2

y <- rpois(100, 6.5)

x2 <- seq(-0.5, 0.5,,length(y))
x1 <- (log(y) - B0 - B2 * x2) / B1

my.model <- glm(y ~ x1 + x2, family = poisson(link = log))
summary(my.model)

Que retorna:

Call:
glm(formula = y ~ x1 + x2, family = poisson(link = log))

Deviance Residuals: 
       Min          1Q      Median          3Q         Max  
-2.581e-08  -1.490e-08   0.000e+00   0.000e+00   4.215e-08  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept)  1.20000    0.08386  14.309  < 2e-16 ***
x1           1.50000    0.16839   8.908  < 2e-16 ***
x2          -0.50000    0.14957  -3.343 0.000829 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 8.8619e+01  on 99  degrees of freedom
Residual deviance: 1.1102e-14  on 97  degrees of freedom
AIC: 362.47

Number of Fisher Scoring iterations: 3

— Mark Miller
fonte