Advertisement
Guest User

Untitled

a guest
Mar 18th, 2019
83
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.83 KB | None | 0 0
  1. > zinf.poiss <- glmmadmb(resp ~ var1 + var2 - 1 + (1|var3),
  2. data = exampledata, family = 'poisson', zeroInflation = TRUE)
  3.  
  4. > summary(zinf.poiss)
  5.  
  6. Call:
  7. glmmadmb(formula = resp ~ var1 + var2 - 1 + (1 | var3), data = exampledata,
  8. family = "poisson", zeroInflation = TRUE)
  9.  
  10. AIC: 2124.5
  11.  
  12. Coefficients:
  13. Estimate Std. Error z value Pr(>|z|)
  14. var1D -0.552 0.197 -2.80 0.0051 **
  15. var1J -0.314 0.207 -1.52 0.1295
  16. var22016 -0.516 0.162 -3.19 0.0014 **
  17. var22017 -0.279 0.127 -2.20 0.0281 *
  18. var22018 -0.387 0.132 -2.93 0.0034 **
  19. ---
  20. Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
  21.  
  22. Number of observations: total=1365, var3=12
  23. Random effect variance(s):
  24. Group=var3
  25. Variance StdDev
  26. (Intercept) 0.267 0.5167
  27.  
  28. Zero-inflation: 0.2634 (std. err.: 0.052798 )
  29.  
  30. Log-likelihood: -1055.25
  31. Warning message:
  32. In .local(x, sigma, ...) :
  33. 'sigma' and 'rdig' arguments are present for compatibility only: ignored
  34.  
  35. > zinf.poissref <- as.list(ref_grid(zinf.poiss))
  36. > zinf.poiss.em <- as.emmGrid(zinf.poissref)
  37. > estimates.zinf.poiss <- confint(zinf.poiss.em, adjust = 'none', level = 0.95)
  38.  
  39. > estimates.zinf.poiss
  40. var1 var2 prediction SE df asymp.LCL asymp.UCL
  41. D 2015 -0.5518259 0.1970900 Inf -0.9381152 -0.1655366
  42. J 2015 -0.3136274 0.2068700 Inf -0.7190851 0.0918304
  43. D 2016 -1.0682002 0.2963979 Inf -1.6491294 -0.4872710
  44. J 2016 -0.8300017 0.3561272 Inf -1.5279982 -0.1320052
  45. D 2017 -0.8309125 0.2366763 Inf -1.2947895 -0.3670354
  46. J 2017 -0.5927140 0.2452413 Inf -1.0733781 -0.1120498
  47. D 2018 -0.9386216 0.2334296 Inf -1.3961352 -0.4811081
  48. J 2018 -0.7004231 0.2461742 Inf -1.1829157 -0.2179306
  49.  
  50. Results are given on the log (not the response) scale.
  51. Confidence level used: 0.95
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement