Advertisement
Guest User

Untitled

a guest
Feb 22nd, 2017
56
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.50 KB | None | 0 0
  1. Generalized linear mixed model fit by maximum likelihood (Laplace Approximation) ['glmerMod']
  2. Family: binomial ( logit )
  3. Formula: survival ~ tree + (1 | tree/rep)
  4. Data: NPV01Datacorr.data
  5. Control: glmerControl(optimizer = "bobyqa")
  6.  
  7. ...
  8.  
  9. Fixed effects:
  10. Estimate Std. Error z value Pr(>|z|)
  11. (Intercept) -0.09531 0.17837 -0.534 0.59311
  12. treeBC3F3 0.84765 0.26019 3.258 0.00112 **
  13. treeD54 -0.40113 0.24606 -1.630 0.10305
  14. treeD58 -0.76691 0.26314 -2.914 0.00356 **
  15. treeEllis1 -0.38426 0.25090 -1.532 0.12563
  16. treeQing 0.53357 0.24338 2.192 0.02836 *
  17. ---
  18. Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
  19.  
  20. levels(NPV01Datacorr.data$tree)
  21. [1] "Ash" "BC3F3" "D54" "D58" "Ellis1" "Qing"
  22.  
  23. contrastmatrix <- rbind(c(-0.5, 0, 0.5, 0.5, -0.5, 0), c(0, 0, 0.5, 0.5, -1, 0))
  24. rownames(contrastmatrix) <- c("Transgenic - Wildtype", "Transgenic - Isogenic")
  25. summary(glht(NPVcorrbinary.glmer, linfct=mcp(tree=contrastmatrix)))
  26.  
  27. Simultaneous Tests for General Linear Hypotheses
  28.  
  29. Multiple Comparisons of Means: User-defined Contrasts
  30.  
  31. Linear Hypotheses:
  32. Estimate Std. Error z value Pr(>|z|)
  33. Transgenic - Wildtype == 0 -0.3919 0.1797 -2.181 0.0459 *
  34. Transgenic - Isogenic == 0 -0.1998 0.2183 -0.915 0.4866
  35.  
  36. (D54 estimate(0.5)) + (D58 estimate(0.5)) + (Intercept estimate(-0.5)) + (Ellis1 estimate(-0.5))
  37. (-0.40113(0.5)) + (-0.76691(0.5)) + (-0.9531(-0.5)) + (-0.38426(-0.5))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement