Guest User

Untitled

a guest
Mar 21st, 2018
72
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.60 KB | None | 0 0
  1. library(metafor)
  2. m1i <- 15.6
  3. m2i <- 12.2
  4. n1i <- 15
  5. n2i <- 20
  6. sd1i <- 3.82
  7. sd2i <- 3.22
  8.  
  9. log(m1i) - log(m2i)
  10. # [1] 0.245835
  11. log(m1i + .5) - log(m2i + .5)
  12. # [1] 0.2372173
  13.  
  14. sd1i^2/(n1i * (m1i)) + sd2i^2/(n2i * (m2i))
  15. # [1] 0.1048541
  16. sd1i^2/(n1i * (m1i + .5)) + sd2i^2/(n2i * (m2i + .5))
  17. # [1] 0.1012445
  18.  
  19.  
  20. lnR <- escalc(measure="ROM", m1i=m1i, m2i=m2i, sd1i=sd1i, sd2i=sd2i, n1i=n1i, n2i=n2i)
  21. lnRk <- escalc(measure="ROM", add=0.5, to="all", m1i=m1i, m2i=m2i, sd1i=sd1i, sd2i=sd2i, n1i=n1i, n2i=n2i)
  22. lnR; lnRk
  23. # yi vi
  24. # 1 0.2458 0.0075
  25. # yi vi
  26. # 1 0.2458 0.0075
  27. all.equal(lnR, lnRk)
  28. # [1] TRUE
Add Comment
Please, Sign In to add comment