Skip to content

Instantly share code, notes, and snippets.

@drizopoulos
Created August 27, 2018 08:10
fm2 <- update(fm1, zi_random = ~ 1 | id)
fm2
#>
#> Call:
#> mixed_model(fixed = y ~ sex * time, random = ~1 | id, data = DF,
#> family = zi.poisson(), zi_fixed = ~sex, zi_random = ~1 |
#> id)
#>
#>
#> Model:
#> family: zero-inflated poisson
#> link: log
#>
#> Random effects covariance matrix:
#> StdDev Corr
#> (Intercept) 0.7391
#> zi_(Intercept) 0.7157 -0.6334
#>
#> Fixed effects:
#> (Intercept) sexfemale time sexfemale:time
#> 1.5902225072 -0.0117290754 -0.0029500093 0.0001046371
#>
#> Zero-part coefficients:
#> (Intercept) sexfemale
#> -1.1201182 0.4416912
#>
#> log-Lik: -2214.306
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment