Guest User

Untitled

a guest
Feb 23rd, 2018
101
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.58 KB | None | 0 0
  1. if(!require(coin)){install.packages("coin")}
  2.  
  3. set.seed(169)
  4.  
  5. n = 16
  6.  
  7. N = n + n
  8.  
  9. A = runif(n,1,10)
  10. B = A + 5 + rnorm(n,0,4)
  11.  
  12. Group = factor(c(rep("A", length(A)), rep("B", length(B))))
  13. Y = c(A,B)
  14.  
  15. boxplot(Y ~ Group)
  16.  
  17. library(coin)
  18.  
  19. MW = wilcox_test(Y ~ Group)
  20. MW
  21.  
  22. Z = MW@statistic@teststatistic
  23. names(Z) = "Z"
  24.  
  25. Z
  26.  
  27. ### Z
  28. ### -4.258848
  29.  
  30. MWa = wilcox.test(Y ~ Group)
  31.  
  32. MWa
  33.  
  34. Za = qnorm(MWa$p.value)
  35.  
  36. Za
  37.  
  38. ### [1] -4.584451
  39.  
  40. r = abs(Z)/sqrt(N)
  41. names(r) = "r"
  42.  
  43. r
  44.  
  45. ### r
  46. ### 0.752865
  47.  
  48. ra = abs(Za)/sqrt(N)
  49. names(ra) = "ra"
  50.  
  51. ra
  52.  
  53. ### ra
  54. ### 0.8104241
Add Comment
Please, Sign In to add comment