Advertisement
Guest User

Untitled

a guest
Jul 1st, 2015
574
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.30 KB | None | 0 0
  1. dat<- read.csv(file="C:/data.csv",head=TRUE,sep=",")
  2. library(metafor)
  3.  
  4. dat<- escalc(measure="ROM", m1i = mean_t, m2i = mean_c, sd1i = sd_t, sd2i = sd_c, n1i = n, n2i = n, data = dat)
  5. metaa<- rma(yi, vi, method="DL", data=dat, mods=cbind(ppm))
  6.  
  7. wi <- 1/sqrt(dat$vi)
  8. size <- 0.5 + 3.0 * (wi - min(wi))/(max(wi) - min(wi))
  9. plot(dat$ppm, exp(dat$yi), pch=19, cex=size, xlab="PPM", ylab="Reaction rate",las=1, bty="l", log="y")
  10.  
  11. preds <- predict(metaa, newmods=c(540:740), transf=exp)
  12. lines(540:740, preds$pred)
  13. lines(540:740, preds$ci.lb, lty="dashed")
  14. lines(540:740, preds$ci.ub, lty="dashed")
  15.  
  16. number mean_c mean_t sd_c sd_t n ppm
  17. 1 1 36.85 48.61 28.40 24.54 20 700
  18. 2 2 31.36 29.01 16.83 21.04 20 700
  19. 3 8 29.00 35.00 3.03 3.03 4 700
  20. 4 9 26.12 41.05 7.50 4.14 12 700
  21. 5 13 38.20 34.90 9.68 11.23 15 550
  22. 6 14 38.20 36.30 9.68 12.01 15 550
  23. 7 15 21.00 55.20 12.07 8.05 20 550
  24. 8 16 62.00 62.00 9.80 9.80 6 700
  25. 9 17 53.00 53.00 7.35 9.80 6 700
  26. 10 18 76.00 63.00 7.35 17.15 6 700
  27. 11 23 258.00 249.00 101.19 199.22 10 700
  28. 12 24 12.00 23.75 6.48 6.48 8 560
  29. 13 25 11.25 20.63 6.48 6.48 8 560
  30. 14 26 17.63 25.75 6.48 6.48 8 560
  31. 15 27 16.38 19.00 6.48 6.48 8 560
  32. 16 46 360.00 360.00 200.92 259.81 12 600
  33. 17 47 170.00 234.00 90.07 62.35 12 600
  34. 18 48 228.00 284.00 38.11 131.64 12 600
  35. 19 49 260.00 340.00 263.27 443.41 12 600
  36. 20 50 75.00 147.00 65.82 110.85 12 600
  37. 21 51 138.00 240.00 110.85 242.49 12 600
  38. 22 52 94.00 157.00 110.85 138.56 12 600
  39. 23 82 154.00 154.00 90.07 31.18 12 540
  40. 24 83 156.00 329.00 110.85 76.21 12 540
  41. 25 84 163.00 293.00 100.46 45.03 12 540
  42. 26 94 376.00 418.00 148.63 132.82 10 740
  43. 27 95 29.00 36.00 41.11 82.22 10 740
  44. 28 96 188.00 403.00 117.00 94.87 10 740
  45. 29 97 121.30 207.80 34.47 43.64 10 700
  46. 30 98 278.30 146.20 82.54 25.93 10 700
  47. 31 120 212.00 226.00 153.36 169.79 30 700
  48. 32 121 568.00 663.00 83.14 121.24 12 550
  49. 33 122 677.00 648.00 131.64 173.21 12 550
  50. 34 123 279.00 449.00 117.00 154.95 10 730
  51. 35 124 266.00 352.00 139.14 211.87 10 730
  52. 36 125 51.66 53.94 52.81 40.49 8 700
  53. 37 126 44.81 44.19 66.90 61.03 8 700
  54. 38 127 14.56 21.10 26.76 17.60 8 700
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement