-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy patheffect-size.R
117 lines (77 loc) · 18.5 KB
/
effect-size.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
## initialization ##
cols <- c('d0.01', 'd0.02', 'd0.04', 'd0.08', 'd0.16')
rows <- c('u25%','u50%','u75%','u100%')
effect_size <- matrix(0.0, nrow=length(rows), ncol=length(cols), dimnames=list(rows, cols))
effect_size1 <- matrix(0.0, nrow=length(rows), ncol=length(cols), dimnames=list(rows, cols))
statistical_diff <- matrix(0.0, nrow=length(rows), ncol=length(cols), dimnames=list(rows, cols))
x <- matrix(0.0, nrow=length(rows), ncol=100) # flat sample
x1 <- matrix(0.0, nrow=length(rows), ncol=100) # rand sample
y <- matrix(0.0, nrow=length(cols)*length(rows), ncol=100) # HDR-aware sample
# *********** u25% *****************
x[1,] <- c(1503, 1503, 2004, 1503, 1503, 1503, 2004, 1002, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 1503, 2004, 2004, 1503, 1503, 2004, 1503, 1503, 1503, 2004, 1002, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 1503, 2004, 2004, 1503, 1503, 2004, 1503, 1503, 1503, 2004, 1002, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 1503, 2004, 2004, 1503, 1503, 2004, 1503, 1503, 1503, 2004, 1002, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 1503, 2004, 2004, 1503, 1503, 2004, 1503, 1503, 1503, 2004, 1002, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 1503, 2004, 2004)
x1[1,] <- c(2004, 1503, 2004, 1503, 1503, 2004, 2004, 2004, 2004, 1503, 1503, 1503, 2004, 1503, 1503, 2004, 2004, 2004, 1503, 2004, 2004, 2004, 1503, 1503, 2004, 2004, 2004, 1503, 1503, 2004, 2004, 2004, 2004, 1503, 1503, 1503, 1002, 2004, 1503, 2004, 2004, 2004, 1503, 1503, 1503, 2004, 2004, 2004, 1002, 1503, 2004, 1503, 2004, 1503, 2004, 1503, 1503, 2004, 1503, 1503, 2004, 1503, 1503, 1503, 1002, 1503, 2004, 2004, 2004, 1503, 2004, 2004, 2004, 2004, 1503, 1503, 1503, 2004, 1503, 2004, 2004, 1503, 2004, 2004, 2004, 2004, 1503, 1503, 1503, 2004, 1503, 2004, 2004, 2004, 2004, 1002, 2004, 1503, 1503, 2004)
# dist 0.01
y[1,] <- c(1999, 2499, 2499, 1499, 1499, 1999, 1999, 2499, 1999, 1499, 1499, 1999, 1999, 1499, 1999, 1999, 1999, 2499, 1499, 2499, 1499, 2499, 1499, 1499, 1999, 1999, 1999, 2499, 1999, 1499, 1999, 1499, 2499, 2499, 2499, 1499, 1999, 1999, 1999, 1999, 2499, 1499, 2499, 2499, 1499, 2499, 2499, 1499, 1499, 2499, 1999, 1499, 1999, 1499, 1499, 1499, 1999, 1499, 1499, 2499, 1999, 1999, 1499, 1499, 1999, 1499, 2499, 1499, 1499, 2499, 1499, 1499, 2499, 1999, 2499, 2499, 2499, 2499, 1999, 1999, 1999, 1999, 2499, 1499, 1999, 2499, 2499, 1499, 2499, 1499, 1999, 1999, 1999, 1999, 2499, 1999, 1499, 1499, 2499, 1499)
# dist 0.02
y[2,] <- c(3499, 999, 1999, 1499, 3999, 2499, 2499, 1999, 1999, 499, 5499, 499, 1999, 2999, 1999, 5499, 2499, 999, 5499, 1999, 3999, 3499, 999, 499, 5499, 1499, 999, 2999, 999, 1999, 1499, 2499, 3499, 999, 1999, 999, 2999, 1999, 1499, 2499, 5499, 1499, 2999, 999, 2499, 999, 1999, 1999, 999, 499, 999, 499, 1999, 1499, 499, 999, 2999, 999, 999, 1999, 2499, 999, 1999, 499, 499, 999, 1999, 999, 499, 1999, 2999, 2499, 1999, 1999, 3999, 499, 999, 999, 999, 999, 1999, 1499, 3499, 2999, 2999, 2999, 1999, 2499, 3999, 1499, 999, 2499, 999, 2999, 3999, 999, 3499, 1499, 999, 1499)
# dist 0.04
y[3,] <- c(999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 1499, 1499, 999, 1499, 999, 999, 1499, 1499, 999, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 1499, 1499, 999, 999, 999, 999, 999, 999, 999, 1499, 1499, 1499, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 1499, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 1499, 1499, 1499, 999, 1499, 999, 1499, 999, 1499, 999, 999, 999, 1499)
# dist 0.08
y[4,] <- c(499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499)
# dist 0.16
y[5,] <- c(499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499, 499)
# *********** u50% *****************
x[2,] <- c(8116, 10680, 7992, 10794, 7981, 10788, 10568, 10673, 8020, 10727, 7967, 8176, 8041, 10800, 10773, 10659, 8067, 10692, 10735, 10692, 8116, 10680, 7992, 10794, 7981, 10788, 10568, 10673, 8020, 10727, 7967, 8176, 8041, 10800, 10773, 10659, 8067, 10692, 10735, 10692, 8116, 10680, 7992, 10794, 7981, 10788, 10568, 10673, 8020, 10727, 7967, 8176, 8041, 10800, 10773, 10659, 8067, 10692, 10735, 10692, 8116, 10680, 7992, 10794, 7981, 10788, 10568, 10673, 8020, 10727, 7967, 8176, 8041, 10800, 10773, 10659, 8067, 10692, 10735, 10692, 8116, 10680, 7992, 10794, 7981, 10788, 10568, 10673, 8020, 10727, 7967, 8176, 8041, 10800, 10773, 10659, 8067, 10692, 10735, 10692)
x1[2,] <- c(10568, 7967, 10659, 10568, 10568, 10568, 10568, 10568, 10680, 8116, 10692, 7981, 10659, 10800, 10568, 10727, 10568, 10568, 10568, 10788, 10773, 10692, 7967, 7981, 10794, 7992, 7981, 7992, 10735, 10568, 8067, 8067, 10794, 8176, 10692, 8020, 10568, 10568, 10773, 10692, 10800, 10680, 10673, 10568, 10659, 10568, 10568, 8116, 10773, 10568, 10794, 10735, 8176, 8176, 10673, 10568, 7992, 10692, 8067, 10568, 10568, 10680, 10735, 8116, 10673, 8041, 10568, 8067, 10568, 10727, 7967, 10692, 10788, 10788, 7981, 10568, 10568, 8116, 8041, 10568, 10794, 10568, 10568, 8041, 10800, 10788, 10800, 8020, 10659, 10692, 10692, 10680, 10568, 7992, 10568, 10773, 8020, 10727, 10735, 8041)
# dist 0.01
y[6,] <- c(1499, 5999, 1499, 2999, 4499, 2999, 2999, 2999, 1499, 1499, 2999, 2999, 1499, 4499, 2999, 1499, 1499, 2999, 1499, 1499, 2999, 1499, 2999, 1499, 2999, 2999, 1499, 1499, 2999, 2999, 2999, 1499, 2999, 1499, 2999, 2999, 1499, 2999, 1499, 2999, 2999, 2999, 2999, 4499, 2999, 1499, 1499, 1499, 2999, 2999, 2999, 1499, 1499, 1499, 2999, 1499, 5999, 2999, 2999, 1499, 2999, 2999, 2999, 1499, 4499, 1499, 2999, 4499, 1499, 2999, 1499, 1499, 2999, 2999, 2999, 5999, 1499, 5999, 1499, 2999, 1499, 1499, 2999, 2999, 1499, 2999, 1499, 2999, 2999, 1499, 2999, 2999, 1499, 5999, 2999, 2999, 1499, 2999, 2999, 1499)
# dist 0.02
y[7,] <- c(1999, 2999, 2999, 1999, 1999, 2499, 3999, 1999, 1999, 1999, 2499, 1999, 1999, 2999, 1999, 1999, 1999, 1999, 1999, 1999, 2999, 2999, 1999, 2999, 1999, 2999, 1999, 1999, 2499, 1999, 2499, 3999, 3999, 1999, 2499, 1999, 1999, 1999, 3999, 2999, 1999, 3999, 2999, 2999, 2999, 2499, 1999, 2499, 2999, 1999, 1999, 1999, 2999, 1999, 2999, 2999, 1999, 2999, 1999, 3999, 1999, 2499, 1999, 1999, 2999, 1999, 2999, 2999, 1999, 2499, 2499, 1999, 2999, 1999, 1999, 2999, 1999, 1999, 2499, 1999, 3999, 3999, 3999, 2999, 1999, 1999, 2499, 2999, 3999, 2999, 2499, 2499, 1999, 1999, 1999, 1999, 2999, 2499, 1999, 1999)
# dist 0.04
y[8,] <- c(1499, 1499, 1499, 1499, 1499, 1999, 1999, 1999, 1999, 1499, 1499, 1499, 1999, 1999, 1499, 999, 1499, 1499, 1999, 1999, 1499, 1499, 1499, 1499, 1999, 1999, 999, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1999, 999, 1499, 1499, 1499, 1499, 1499, 1999, 1499, 999, 1499, 1499, 1999, 1499, 1999, 999, 999, 1999, 1999, 1499, 1499, 1999, 1499, 1499, 1999, 1499, 1499, 1999, 1499, 1999, 1499, 1999, 1499, 999, 999, 1999, 1999, 1999, 1499, 1499, 1499, 1999, 1499, 1499, 1499, 999, 1999, 1499, 1499, 1499, 1499, 1499, 999, 1499, 1999, 1499, 1499, 1999, 1499, 1499, 1499, 1999, 1999)
# dist 0.08
y[9,] <- c(999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 1499, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999)
# dist 0.16
y[10,] <- c(999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999, 999)
# *********** u75% *****************
x[3,] <- c(25866, 26577, 26456, 25780, 26196, 25873, 25842, 26523, 26538, 34420, 34892, 35364, 26144, 35067, 34563, 26465, 26427, 35032, 26534, 25987, 26717, 26041, 26641, 25772, 26584, 26695, 25828, 25939, 26515, 35168, 34832, 35231, 26005, 34779, 35265, 26188, 26254, 34515, 26390, 26070, 25865, 26082, 26505, 26346, 26705, 26417, 25828, 26356, 26243, 35148, 35036, 35136, 25802, 34710, 35138, 25918, 26811, 35112, 26038, 25889, 25886, 26672, 26756, 26060, 26866, 26148, 26299, 26092, 26601, 35023, 34487, 34616, 25973, 34730, 34566, 26395, 26762, 34656, 26154, 26031, 26309, 26854, 26740, 26437, 26685, 26310, 25960, 26611, 26344, 34718, 34253, 35204, 26388, 34496, 35070, 26452, 26311, 34655, 26503, 26705)
x1[3,] <- c(26196, 25939, 35168, 26611, 35168, 35112, 34718, 26762, 25886, 35168, 35168, 26705, 26196, 35138, 26601, 26196, 26717, 26196, 35168, 25918, 26740, 25987, 26515, 35168, 26196, 34496, 26344, 26196, 26092, 26310, 25973, 26854, 35168, 35168, 35168, 34566, 26060, 26196, 26577, 26866, 35032, 26452, 26437, 35168, 35070, 35168, 26584, 26641, 26534, 35168, 34487, 26311, 26503, 35204, 26196, 35168, 26456, 34655, 26685, 35168, 35023, 26196, 26154, 35168, 26672, 34656, 26695, 26395, 25772, 26756, 26196, 25866, 26196, 35168, 34710, 26811, 25828, 35168, 26041, 26299, 35168, 35168, 34253, 26196, 26148, 25889, 25802, 25780, 26031, 26309, 26196, 26388, 34616, 35168, 26038, 35168, 25960, 35168, 34730, 35136)
# dist 0.01
y[11,] <- c(1999, 1999, 1999, 3999, 3999, 1999, 3999, 3999, 2999, 3999, 7999, 3999, 2999, 2999, 1999, 3999, 3999, 3999, 3999, 3999, 2999, 3999, 3999, 2999, 1999, 1999, 2999, 3999, 3999, 5999, 3999, 1999, 3999, 1999, 3999, 1999, 1999, 3999, 2999, 1999, 3999, 7999, 2999, 7999, 3999, 3999, 3999, 3999, 1999, 3999, 1999, 1999, 2999, 3999, 3999, 3999, 3999, 5999, 3999, 3999, 3999, 3999, 3999, 2999, 1999, 1999, 1999, 3999, 1999, 7999, 1999, 1999, 2999, 3999, 2999, 2999, 3999, 5999, 2999, 3999, 2999, 1999, 3999, 2999, 7999, 2999, 3999, 3999, 3999, 5999, 5999, 2999, 2999, 3999, 3999, 1999, 2999, 3999, 1999, 1999)
# dist 0.02
y[12,] <- c(2999, 2499, 2499, 2999, 3999, 3499, 2499, 1999, 2499, 2999, 2499, 2999, 2499, 2499, 1999, 1999, 2499, 2999, 2499, 2499, 1999, 2499, 3499, 2499, 2499, 2999, 2499, 2499, 2999, 2999, 2999, 2499, 2999, 3999, 2999, 2499, 2999, 2499, 2499, 3499, 2499, 2499, 2999, 3499, 2999, 3999, 2499, 3499, 2499, 2499, 3999, 2999, 3999, 2499, 2999, 2499, 2999, 3499, 2499, 3499, 3999, 3499, 3499, 2499, 2499, 2499, 2499, 3499, 3999, 3999, 3999, 2999, 2999, 2499, 3499, 2499, 2499, 3999, 3499, 1999, 3499, 2499, 2499, 3499, 2999, 2999, 2499, 3499, 2999, 3499, 3499, 2999, 3499, 2499, 3499, 3499, 2499, 2999, 2999, 2499)
# dist 0.04
y[13,] <- c(1999, 1499, 1999, 1999, 1999, 1499, 1999, 1999, 1499, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1499, 1999, 1499, 1999, 1999, 1499, 1999, 1999, 1999, 1499, 1999, 1499, 1499, 1999, 1999, 1999, 1999, 1499, 1999, 1999, 1499, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1999, 1499, 1499, 1499, 1499, 1999, 1499, 1999, 1999, 1999, 1499, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1499, 1999)
# dist 0.08
y[14,] <- c(1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499)
# dist 0.16
y[15,] <- c(1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499, 1499)
# *********** u100% *****************
x[4,] <- c(220237, 216083, 300061, 299479, 222083, 291465, 220216, 296045, 220969, 225058, 227726, 227675, 220610, 220440, 302262, 226944, 218520, 298324, 229519, 225646, 221691, 217040, 299381, 300678, 221821, 290827, 222422, 295638, 222340, 223678, 229641, 225303, 221672, 220800, 306597, 225755, 221548, 302612, 229313, 224138, 221382, 216780, 296777, 300041, 220981, 291743, 222426, 296575, 221976, 226529, 228766, 226353, 223259, 219185, 305071, 223953, 222177, 302358, 228824, 228437, 222424, 218353, 298574, 298298, 222895, 292290, 217959, 298071, 221893, 226569, 226954, 228737, 223004, 220671, 306708, 226234, 222480, 299818, 229601, 225144, 219159, 218095, 297682, 299211, 223328, 291187, 219280, 297989, 223075, 226777, 229377, 227589, 220660, 219752, 304209, 226479, 221955, 302114, 225347, 225797)
x1[4,] <- c(221893, 300061, 302262, 295638, 297989, 218353, 218520, 300678, 216083, 225646, 229377, 298574, 217040, 222422, 296575, 229641, 291187, 220237, 226954, 300061, 299818, 296575, 225347, 300061, 300061, 296575, 221691, 221672, 296575, 223004, 292290, 300061, 220671, 221548, 300061, 298071, 296575, 223328, 229313, 302114, 302114, 302114, 296575, 226569, 216780, 290827, 296575, 300041, 299381, 291743, 298324, 224138, 226777, 217959, 222895, 226234, 298298, 300061, 306708, 221382, 226944, 218095, 296575, 296575, 229601, 300061, 306597, 223075, 219159, 302114, 296575, 302114, 220981, 222340, 228737, 299211, 229519, 296777, 225303, 300061, 220800, 223678, 300061, 302114, 220440, 225144, 225755, 302612, 219280, 221821, 225797, 296575, 296575, 300061, 297682, 222480, 302114, 222424, 296575, 222426)
# dist 0.01
y[16,] <- c(1999, 7999, 1999, 2999, 7999, 2999, 7999, 5999, 5999, 5999, 5999, 1999, 1999, 5999, 3999, 2999, 5999, 2999, 1999, 1999, 7999, 8999, 3999, 3999, 1999, 7999, 2999, 5999, 5999, 5999, 2999, 5999, 1999, 1999, 5999, 1999, 5999, 1999, 2999, 6999, 5999, 3999, 5999, 7999, 2999, 5999, 5999, 3999, 1999, 7999, 6999, 2999, 5999, 3999, 1999, 2999, 5999, 1999, 5999, 7999, 6999, 5999, 1999, 5999, 7999, 5999, 5999, 5999, 3999, 1999, 2999, 3999, 1999, 7999, 2999, 8999, 3999, 8999, 8999, 1999, 5999, 6999, 5999, 3999, 2999, 5999, 5999, 5999, 3999, 8999, 3999, 5999, 2999, 1999, 1999, 2999, 3999, 3999, 3999, 6999)
# dist 0.02
y[17,] <- c(2999, 4999, 3499, 4499, 2999, 2999, 2999, 2999, 4999, 2999, 2999, 2999, 2999, 2999, 3499, 3499, 4499, 2499, 2999, 2499, 2999, 4499, 3499, 3499, 2999, 3499, 4999, 2499, 2999, 3499, 4999, 3499, 4499, 4499, 2999, 3499, 2999, 2999, 2999, 3499, 2999, 3499, 4999, 3499, 3499, 3499, 3499, 4999, 2999, 2999, 3499, 4999, 4999, 4499, 2999, 2999, 4499, 3499, 4499, 2999, 4499, 4499, 2999, 3499, 4499, 4999, 2999, 3499, 2999, 2999, 2999, 2999, 4499, 3499, 2999, 2999, 2499, 2999, 4999, 2999, 3499, 2999, 4499, 4499, 4499, 2999, 2999, 3499, 2999, 3499, 2999, 3499, 2999, 2999, 2999, 3499, 2999, 2499, 2999, 2999)
# dist 0.04
y[18,] <- c(1999, 1999, 2499, 2499, 2499, 2499, 2499, 2499, 2999, 2499, 2499, 2499, 2499, 2499, 2499, 1999, 2499, 1999, 2499, 2499, 2499, 2499, 2499, 2499, 1999, 2499, 2499, 2499, 2499, 2499, 2499, 2499, 2499, 1999, 2999, 2499, 2499, 2499, 2499, 2499, 1999, 2499, 2499, 2499, 2499, 2499, 1999, 1999, 2499, 2499, 1999, 1999, 2499, 2499, 1999, 1999, 2499, 1999, 1999, 2999, 2499, 2999, 1999, 1999, 1999, 2499, 2499, 2499, 2499, 2499, 2499, 2499, 2499, 1999, 2499, 2499, 1999, 2499, 2499, 2499, 1999, 2499, 2499, 2499, 2499, 1999, 1999, 2499, 2499, 2499, 2499, 1999, 1999, 2499, 2999, 2499, 2499, 2499, 2499, 2499)
# dist 0.08
y[19,] <- c(1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999)
# dist 0.16
y[20,] <- c(1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999, 1999)
# *********** A12 effect-size + wilcox statistical-diff computation **************
for (i in 1:length(rows)) {
for (j in 1:length(cols)) {
k <- (i-1)*length(cols) + j
m <- length(x[i,])
n <- length(y[k,])
r <- sum(rank(c(x[i,],y[k,]))[seq_along(x[i,])])
r1 <- sum(rank(c(x1[i,],y[k,]))[seq_along(x1[i,])])
effect_size[i, j] <- (r/m - (m + 1)/2)/n
effect_size1[i, j] <- (r1/m - (m + 1)/2)/n
statistical_diff[i, j] <- wilcox.test(x[i,], y[k,])$p.value
}
}
effect_size
effect_size1
statistical_diff