Time-series plots:Full dataset:Truncated dataset:
Basic statistics:Full dataset (only dropped first three days):. tabstat merit, s(n mean sd p50 p25 p75 min max) format(%9.1f)
variable | N mean sd p50 p25 p75 min max
-------------+--------------------------------------------------------------------------------
merit | 564.0 745.6 443.5 641.0 534.0 777.5 295.0 4500.0
----------------------------------------------------------------------------------------------
Applied formulas in previous weeks, potential outliers are days have intra-day merits beyond 169 or 1143.
. di 777.5-534
243.5
. di 243.5*1.5
365.25
. di 777.5+365.25
1142.75
. di 534-365.25
168.75
There are
49 outliers in full dataset, in total.
. count if (merit >= 1143 | merit <= 169) & merit != .
49
Those days are:
. list id merit date if (merit >= 1143 | merit <= 169) & merit != .
+-------------------------+
| id merit date |
|-------------------------|
1. | 3 4500 26jan2018 |
2. | 4 3490 27jan2018 |
3. | 5 3190 28jan2018 |
4. | 6 3800 29jan2018 |
5. | 7 4193 30jan2018 |
|-------------------------|
6. | 8 2821 31jan2018 |
7. | 9 2546 01feb2018 |
8. | 10 2569 02feb2018 |
9. | 11 1868 03feb2018 |
10. | 12 2182 04feb2018 |
|-------------------------|
11. | 13 2078 05feb2018 |
12. | 14 2310 06feb2018 |
13. | 15 2142 07feb2018 |
14. | 16 2143 08feb2018 |
15. | 17 1449 09feb2018 |
|-------------------------|
16. | 18 1748 10feb2018 |
17. | 19 1443 11feb2018 |
18. | 20 1332 12feb2018 |
19. | 21 1580 13feb2018 |
20. | 22 2514 14feb2018 |
|-------------------------|
21. | 23 1992 15feb2018 |
22. | 24 1416 16feb2018 |
23. | 25 1618 17feb2018 |
24. | 26 1293 18feb2018 |
25. | 27 1404 19feb2018 |
|-------------------------|
26. | 28 1170 20feb2018 |
27. | 29 1268 21feb2018 |
28. | 30 1280 22feb2018 |
30. | 32 1410 24feb2018 |
31. | 33 1187 25feb2018 |
|-------------------------|
32. | 34 1392 26feb2018 |
33. | 35 1327 27feb2018 |
35. | 37 1335 01mar2018 |
36. | 38 1706 02mar2018 |
39. | 41 1246 05mar2018 |
|-------------------------|
46. | 48 1355 12mar2018 |
48. | 50 1160 14mar2018 |
54. | 56 1324 20mar2018 |
55. | 57 1229 21mar2018 |
66. | 68 1258 01apr2018 |
|-------------------------|
67. | 69 1147 02apr2018 |
234. | 236 2464 16sep2018 |
235. | 237 1863 17sep2018 |
236. | 238 1295 18sep2018 |
237. | 239 1271 19sep2018 |
|-------------------------|
349. | 351 1162 09jan2019 |
426. | 428 1250 27mar2019 |
473. | 475 1151 13may2019 |
502. | 504 1188 11jun2019 |
+-------------------------+
Only
four of them occured in 2019, on 09/1/2019, 27/3/2019, 13/5/2019, and 11/6/2019, at 1162, 1250, 1151, and 1188 merits circulated in total, respectively.
Truncated dataset (first 25 days dropped):.
. tabstat merit, s(n mean sd p50 p25 p75 min max) format(%9.1f)
variable | N mean sd p50 p25 p75 min max
-------------+--------------------------------------------------------------------------------
merit | 540.0 674.6 226.5 628.0 530.0 757.0 295.0 2464.0
----------------------------------------------------------------------------------------------
Applied same formulas I used in earlier analyses, potential outliers are days have intra-day merits beyond
190 or
1098.
. di 757-530
227
. di 227*1.5
340.5
. di 757+340.5
1097.5
. di 530-340.5
189.5
There are
29 outliers in total, only
five of them occured in 2019, on 09/1/2019, 14/01/2019, 27/3/2019, 13/5/2019, and 11/6/2019, at 1162, 1128, 1250, 1151, and 1188, respectively.
. count if (merit >= 1098 | merit <= 190) & merit != .
29
List of those 29 outliers in truncated dataset
. list id merit date if (merit >= 1098 | merit <= 190) & merit != .
+-------------------------+
| id merit date |
|-------------------------|
1. | 27 1404 19feb2018 |
2. | 28 1170 20feb2018 |
3. | 29 1268 21feb2018 |
4. | 30 1280 22feb2018 |
6. | 32 1410 24feb2018 |
|-------------------------|
7. | 33 1187 25feb2018 |
8. | 34 1392 26feb2018 |
9. | 35 1327 27feb2018 |
11. | 37 1335 01mar2018 |
12. | 38 1706 02mar2018 |
|-------------------------|
15. | 41 1246 05mar2018 |
17. | 43 1111 07mar2018 |
22. | 48 1355 12mar2018 |
24. | 50 1160 14mar2018 |
25. | 51 1131 15mar2018 |
|-------------------------|
30. | 56 1324 20mar2018 |
31. | 57 1229 21mar2018 |
42. | 68 1258 01apr2018 |
43. | 69 1147 02apr2018 |
127. | 153 1139 25jun2018 |
|-------------------------|
210. | 236 2464 16sep2018 |
211. | 237 1863 17sep2018 |
212. | 238 1295 18sep2018 |
213. | 239 1271 19sep2018 |
325. | 351 1162 09jan2019 |
|-------------------------|
330. | 356 1128 14jan2019 |
402. | 428 1250 27mar2019 |
449. | 475 1151 13may2019 |
478. | 504 1188 11jun2019 |
+-------------------------+