Graph with day/month pair on x-axis

[Marine Jouvin]

Hi,

I'm trying to plot a variable over the date, defined as a day/month pair. The x-axis would therefore plot values from January 1 to December 31. I'd like three lines to be plotted on the graph, one for each different year (2020, 2021, 2022).

I have a "date" variable in numeric format (dmy, %td); but I'm having trouble isolating day/month couples in order to use it as x-values. I already have variables isolating month (month) and day (day). In my current dataset, I don't have a value for each day/month pair over the year, so when I'm using the following code, I don't have values from 0 to 365 (which I would like) :

egen couple=group(month day)

Any idea?

Many thanks for any help,

[Andrew Musau]

You want year to be your panel identifier. The following has only 1 panel but is illustrative:

webuse sp500, clear
gen year= year(date)
xtset year date
set scheme s1mono
format date %tdDm
xtline volume, byopts(note("")) xtitle("")

[Marine Jouvin]

Many thanks Andrew!

I've tried to reproduce your graph with my data. The "year" variable takes any value between 0,1,2, and the "date" variable is a float variable (format %tdDm).

By changing nothing in the code you've provided, except the variable I'm plotting, I've obtained the following graphs:



Which makes me think there's something wrong with my date variable...

Also, is there a way of displaying these three lines on a single graph ?

[Andrew Musau]

Also, is there a way of displaying these three lines on a single graph ?

You do this using the -overlay- option of xtline.

I've tried to reproduce your graph with my data. The "year" variable takes any value between 0,1,2, and the "date" variable is a float variable (format %tdDm).
By changing nothing in the code you've provided, except the variable I'm plotting, I've obtained the following graphs:

Show us a data example using the dataex command so that we can look at what is going on.

[Marine Jouvin]

Show us a data example using the dataex command so that we can look at what is going on.

Here is a data example for my variable "date" :

* Example generated by -dataex-. To install: ssc install dataex
clear
input float date
22189
22190
22191
22192
22193
22194
22195
22196
22197
22198
22199
22200
22201
22202
22203
22204
22205
22206
22207
22208
22209
22210
22211
22212
22213
22214
22215
22216
22217
22218
22219
22220
22221
22222
22223
22224
22225
22226
22227
22228
22229
22230
22231
22232
22233
22234
22235
22236
22237
22238
22239
22240
22241
22242
22243
22244
22245
22246
22247
22248
22249
22250
22251
22252
22253
22254
22255
22256
22257
22258
22259
22260
22261
22262
22263
22264
22265
22266
22267
22268
22269
22270
22271
22272
22273
22274
22275
22276
22277
22278
22279
22280
22281
22282
22283
22284
22285
22286
22287
22288
end
format %tdDm date

[Maarten Buis]

// get some example data
clear
import delimited using "https://raw.githubusercontent.com/CSSEGISandData/COVID-19/master/csse_covid_19_data/csse_covid_19_time_series/time_series_covid19_confirmed_global.csv", clear
drop v4 lat
keep if countryregion == "Belgium"
drop countryregion provincestate
xpose , clear
gen date = _n + mdy(1,21,2020)
format date %td
drop if year(date) >= 2023
tsset date
gen newcases = S7.v1

// here the real work starts
gen daymonth = mdy(month(date),day(date), 2020)
gen year = year(date)
format daymonth %tdDDMon

separate newcases, by(year) veryshortlabel
twoway line newcases???? daymonth

[Maarten Buis]

Also see this talk by Nick Cox a couple of years ago: https://www.stata.com/meeting/5nasug/cox.zip for a variety of similar ideas

[Andrew Musau]

Here is a data example for my variable "date" :

The date variable alone is not sufficient. We need the variable that you are plotting as well. If you are not able to make progress with Marteen's suggestions, provide a new data example. You can collapse by month so that we get the full set of observations.

gen year= year(date)
gen month= month(date)
*CHANGE TO THE NAME OF YOUR VARIABLE (HERE CALLED "VAR")
collapse VAR, by(year month)
dataex year month VAR

[Marine Jouvin]

Thank you both for your help!

I tried to reproduce the graph suggested by Marteen and ended up with something quite close to that expected. My data cover a period from October to April, so I'd like to be able to start the x-axis at October and end it at April, to avoid the gaps currently appearing on the graph between these two periods :



Here are the code that I’m currently using and the description of all the variables:

    gen daymonth = mdy(month(date_trans),day(date_trans), 2020)
    collapse (sum) kilos, by(year daymonth)
    separate kilos, by(year) veryshortlabel
    format daymonth %tdDDMon
    twoway line kilos_0 kilos_1 kilos_2 daymonth

* Example generated by -dataex-. To install: ssc install dataex
clear
input float(year daymonth) double kilos long(kilos_0 kilos_1 kilos_2)
0 21915  71661  71661 . .
0 21916 385205 385205 . .
0 21917 643502 643502 . .
0 21918 643509 643509 . .
0 21919 676360 676360 . .
0 21920 483264 483264 . .
0 21921 433171 433171 . .
0 21922 451833 451833 . .
0 21923 518764 518764 . .
0 21924 563800 563800 . .
0 21925 541312 541312 . .
0 21926 591733 591733 . .
0 21927 472663 472663 . .
0 21928 277247 277247 . .
0 21929 417866 417866 . .
0 21930 605605 605605 . .
0 21931 491667 491667 . .
0 21932 455746 455746 . .
0 21933 519541 519541 . .
0 21934 537373 537373 . .
0 21935 309939 309939 . .
0 21936 280627 280627 . .
0 21937 262320 262320 . .
0 21938 317909 317909 . .
0 21939 299978 299978 . .
0 21940 382633 382633 . .
0 21941 572021 572021 . .
0 21942 308382 308382 . .
0 21943 268690 268690 . .
0 21944 395303 395303 . .
0 21945 206577 206577 . .
0 21946 151037 151037 . .
0 21947 255222 255222 . .
0 21948  80864  80864 . .
0 21949  21935  21935 . .
0 21950  17483  17483 . .
0 21951   7812   7812 . .
0 21952  27121  27121 . .
0 21953   8382   8382 . .
0 21954  17019  17019 . .
0 21955   6675   6675 . .
0 21956   5969   5969 . .
0 21957   8100   8100 . .
0 21958   7308   7308 . .
0 21960  11820  11820 . .
0 21961  18388  18388 . .
0 21962  29949  29949 . .
0 21963  28561  28561 . .
0 21964  14797  14797 . .
0 21965  37500  37500 . .
0 22002  38110  38110 . .
0 22003  39812  39812 . .
0 22004 119093 119093 . .
0 22189 217441 217441 . .
0 22190 214775 214775 . .
0 22191 156388 156388 . .
0 22192 141009 141009 . .
0 22193 215526 215526 . .
0 22194 129711 129711 . .
0 22195 169852 169852 . .
0 22196 204643 204643 . .
0 22197 281729 281729 . .
0 22198 412143 412143 . .
0 22199 324678 324678 . .
0 22200 396663 396663 . .
0 22201 460965 460965 . .
0 22202 398581 398581 . .
0 22203 597762 597762 . .
0 22204 556337 556337 . .
0 22205 479245 479245 . .
0 22206 624982 624982 . .
0 22207 670106 670106 . .
0 22208 555994 555994 . .
0 22209 411386 411386 . .
0 22210 440721 440721 . .
0 22211 541893 541893 . .
0 22212 450873 450873 . .
0 22213 400985 400985 . .
0 22214 500747 500747 . .
0 22215 406071 406071 . .
0 22216 731186 731186 . .
0 22217 461157 461157 . .
0 22218 330291 330291 . .
0 22219 219186 219186 . .
0 22220 343361 343361 . .
0 22221 374581 374581 . .
0 22222 463049 463049 . .
0 22223 370848 370848 . .
0 22224 562520 562520 . .
0 22225 530294 530294 . .
0 22226 449209 449209 . .
0 22227 528547 528547 . .
0 22228 533763 533763 . .
0 22229 477947 477947 . .
0 22230 418616 418616 . .
0 22231 369379 369379 . .
0 22232 487571 487571 . .
0 22233 501288 501288 . .
0 22234 429535 429535 . .
0 22235 374690 374690 . .
end
format %tdDDMon daymonth
label values year year
label def year 0 "2020-2021", modify

[Nick Cox]

See Section 4 in https://journals.sagepub.com/doi/pdf...867X1501500215 on setting up non-calendar years.

[Marine Jouvin]

Thanks Nick !

I'm having trouble fully understanding the way to get my NCYs. I have daily data, and want my graph to display values between October 1st and March 31, so I guess I should use the following code :

generate ncy = cond(ddate >= mdy(10, 1, year(ddate)), year, year - 1)

I was wondering what the "ddate" variable should look like? Should it be day/month, or day/month/year combination?

Again, many thanks for the help.

[Nick Cox]

ddate must be a daily date variable, based on day, month and year. Otherwise it is not comparable with any other daily date variable or expression evaluating to a daily date.

[Marine Jouvin]

Thanks !

I think that I might be missing something here. Following this code :

generate ncy = cond(ddate >= mdy(10, 1, year(ddate)), year, year - 1)

I now have the following time-related variables and the one that I want to plot (kilos) :

* Example generated by -dataex-. To install: ssc install dataex
clear
input float(ncy date month day year kilos)
-1 22201 10 13 2020  900
-1 22274 12 25 2020  793
-1 22297  1 17 2021  317
-1 22211 10 23 2020  196
-1 22233 11 14 2020  227
-1 22265 12 16 2020  966
-1 22201 10 13 2020  852
-1 22249 11 30 2020 1095
-1 22268 12 19 2020 1005
-1 22246 11 27 2020  520
-1 22269 12 20 2020  508
-1 22274 12 25 2020 1197
-1 22249 11 30 2020  923
-1 22216 10 28 2020  723
-1 22289  1  9 2021  931
-1 22235 11 16 2020  795
-1 22265 12 16 2020 1056
-1 22255 12  6 2020 1170
-1 22299  1 19 2021  324
-1 22268 12 19 2020  384
-1 22249 11 30 2020 1104
-1 22255 12  6 2020  644
-1 22269 12 20 2020  512
-1 22293  1 13 2021  369
-1 22300  1 20 2021  435
-1 22217 10 29 2020  495
-1 22237 11 18 2020  345
-1 22248 11 29 2020  737
-1 22222 11  3 2020  111
-1 22301  1 21 2021  169
-1 22255 12  6 2020  506
-1 22265 12 16 2020  600
-1 22243 11 24 2020  280
-1 22201 10 13 2020  979
-1 22268 12 19 2020  441
-1 22286  1  6 2021  303
-1 22254 12  5 2020  781
-1 22227 11  8 2020   49
-1 22292  1 12 2021  221
-1 22208 10 20 2020  182
-1 22246 11 27 2020 1022
-1 22201 10 13 2020  659
-1 22283  1  3 2021  102
-1 22240 11 21 2020  975
-1 22245 11 26 2020  356
-1 22298  1 18 2021  206
-1 22246 11 27 2020 1088
-1 22206 10 18 2020   78
-1 22274 12 25 2020  930
-1 22254 12  5 2020  897
-1 22265 12 16 2020  360
-1 22293  1 13 2021  484
-1 22216 10 28 2020  240
-1 22301  1 21 2021  333
-1 22240 11 21 2020  804
-1 22294  1 14 2021   67
-1 22257 12  8 2020  420
-1 22201 10 13 2020  563
-1 22279 12 30 2020  110
-1 22244 11 25 2020   54
-1 22246 11 27 2020 1008
-1 22268 12 19 2020  620
-1 22226 11  7 2020  178
-1 22227 11  8 2020   72
-1 22292  1 12 2021  162
-1 22255 12  6 2020  816
-1 22284  1  4 2021  988
-1 22274 12 25 2020  549
-1 22240 11 21 2020 1095
-1 22255 12  6 2020  938
-1 22201 10 13 2020 1000
-1 22293  1 13 2021  578
-1 22212 10 24 2020  487
-1 22302  1 22 2021  314
-1 22301  1 21 2021  278
-1 22246 11 27 2020  949
-1 22255 12  6 2020  923
-1 22274 12 25 2020  670
-1 22228 11  9 2020  230
-1 22210 10 22 2020  278
-1 22254 12  5 2020  852
-1 22245 11 26 2020  109
-1 22201 10 13 2020  150
-1 22283  1  3 2021  302
-1 22303  1 23 2021  194
-1 22246 11 27 2020 1184
-1 22212 10 24 2020   56
-1 22300  1 20 2021  196
-1 22284  1  4 2021  259
-1 22201 10 13 2020  543
-1 22285  1  5 2021  326
-1 22240 11 21 2020  962
-1 22297  1 17 2021  133
-1 22227 11  8 2020  362
-1 22287  1  7 2021  388
-1 22240 11 21 2020 1125
-1 22223 11  4 2020  210
-1 22204 10 16 2020  594
-1 22274 12 25 2020  621
-1 22240 11 21 2020  876
end
format %td date

I'm not sure about the way I'm supposed to introduce the NCY variable in my graph, and whether it is possible to introduce it in my twoway line plot.

[Andrew Musau]

I do not get a discontinuous axis with:

* Example generated by -dataex-. To install: ssc install dataex
clear
input float(ncy date month day year kilos)
-1 22201 10 13 2020 1000
-1 22201 10 13 2020  563
-1 22201 10 13 2020  900
-1 22201 10 13 2020  150
-1 22201 10 13 2020  979
-1 22201 10 13 2020  659
-1 22201 10 13 2020  543
-1 22201 10 13 2020  852
-1 22204 10 16 2020  594
-1 22206 10 18 2020   78
-1 22208 10 20 2020  182
-1 22210 10 22 2020  278
-1 22211 10 23 2020  196
-1 22212 10 24 2020  487
-1 22212 10 24 2020   56
-1 22216 10 28 2020  723
-1 22216 10 28 2020  240
-1 22217 10 29 2020  495
-1 22222 11  3 2020  111
-1 22223 11  4 2020  210
-1 22226 11  7 2020  178
-1 22227 11  8 2020  362
-1 22227 11  8 2020   72
-1 22227 11  8 2020   49
-1 22228 11  9 2020  230
-1 22233 11 14 2020  227
-1 22235 11 16 2020  795
-1 22237 11 18 2020  345
-1 22240 11 21 2020  962
-1 22240 11 21 2020 1095
-1 22240 11 21 2020  876
-1 22240 11 21 2020  975
-1 22240 11 21 2020 1125
-1 22240 11 21 2020  804
-1 22243 11 24 2020  280
-1 22244 11 25 2020   54
-1 22245 11 26 2020  109
-1 22245 11 26 2020  356
-1 22246 11 27 2020 1088
-1 22246 11 27 2020 1022
-1 22246 11 27 2020 1008
-1 22246 11 27 2020  520
-1 22246 11 27 2020 1184
-1 22246 11 27 2020  949
-1 22248 11 29 2020  737
-1 22249 11 30 2020 1095
-1 22249 11 30 2020  923
-1 22249 11 30 2020 1104
-1 22254 12  5 2020  852
-1 22254 12  5 2020  897
-1 22254 12  5 2020  781
-1 22255 12  6 2020  923
-1 22255 12  6 2020  644
-1 22255 12  6 2020  506
-1 22255 12  6 2020  816
-1 22255 12  6 2020  938
-1 22255 12  6 2020 1170
-1 22257 12  8 2020  420
-1 22265 12 16 2020  966
-1 22265 12 16 2020  360
-1 22265 12 16 2020  600
-1 22265 12 16 2020 1056
-1 22268 12 19 2020  441
-1 22268 12 19 2020  620
-1 22268 12 19 2020  384
-1 22268 12 19 2020 1005
-1 22269 12 20 2020  512
-1 22269 12 20 2020  508
-1 22274 12 25 2020 1197
-1 22274 12 25 2020  549
-1 22274 12 25 2020  793
-1 22274 12 25 2020  930
-1 22274 12 25 2020  621
-1 22274 12 25 2020  670
-1 22279 12 30 2020  110
-1 22283  1  3 2021  102
-1 22283  1  3 2021  302
-1 22284  1  4 2021  259
-1 22284  1  4 2021  988
-1 22285  1  5 2021  326
-1 22286  1  6 2021  303
-1 22287  1  7 2021  388
-1 22289  1  9 2021  931
-1 22292  1 12 2021  162
-1 22292  1 12 2021  221
-1 22293  1 13 2021  369
-1 22293  1 13 2021  484
-1 22293  1 13 2021  578
-1 22294  1 14 2021   67
-1 22297  1 17 2021  317
-1 22297  1 17 2021  133
-1 22298  1 18 2021  206
-1 22299  1 19 2021  324
-1 22300  1 20 2021  196
-1 22300  1 20 2021  435
-1 22301  1 21 2021  278
-1 22301  1 21 2021  333
-1 22301  1 21 2021  169
-1 22302  1 22 2021  314
-1 22303  1 23 2021  194
end
format %td date

collapse (sum) kilos, by(date year month)
set scheme s1mono
xtset year date
xtline kilos, overlay xlab(,format(%tdDm)) xtitle("") leg(col(1) pos(1) ring(0))

How do you define a year? October to April?? I think you need to collapse your data and give a data example with at least 2 years, possibly using the -count()- option of dataex. In the above, the axis will keep on going as I am using the date variable to plot, so you need to map the date to day of year.

collapse (sum) kilos, by(date year month)
sort date
dataex, count(300)

[Marine Jouvin]

Thanks Andrew!

Here is my code and the graph I get:

 generate month =month(date)
    generate day =day(date)

    collapse (sum) kilos, by(date year month)
    set scheme s1mono
    xtset year date
    xtline kilos, overlay xlab(,format(%tdDm)) xtitle("") leg(col(1) pos(1) ring(0))

Here is how my variables look like:

* Example generated by -dataex-. To install: ssc install dataex
clear
input float(date year month) double kilos
22189 0 10 217441
22190 0 10 214775
22191 0 10 156388
22192 0 10 141009
22193 0 10 215526
22194 0 10 129711
22195 0 10 169852
22196 0 10 204643
22197 0 10 281729
22198 0 10 412143
22199 0 10 324678
22200 0 10 396663
22201 0 10 460965
22202 0 10 398581
22203 0 10 597762
22204 0 10 556337
22205 0 10 479245
22206 0 10 624982
22207 0 10 670106
22208 0 10 555994
22209 0 10 411386
22210 0 10 440721
22211 0 10 541893
22212 0 10 450873
22213 0 10 400985
22214 0 10 500747
22215 0 10 406071
22216 0 10 731186
22217 0 10 461157
22218 0 10 330291
22219 0 10 219186
22220 0 11 343361
22221 0 11 374581
22222 0 11 463049
22223 0 11 370848
22224 0 11 562520
22225 0 11 530294
22226 0 11 449209
22227 0 11 528547
22228 0 11 533763
22229 0 11 477947
22230 0 11 418616
22231 0 11 369379
22232 0 11 487571
22233 0 11 501288
22234 0 11 429535
22235 0 11 374690
22236 0 11 419362
22237 0 11 455256
22238 0 11 457892
22239 0 11 464516
22240 0 11 431156
22241 0 11 590224
22242 0 11 502984
22243 0 11 703208
22244 0 11 695350
22245 0 11 818355
22246 0 11 822102
22247 0 11 752457
22248 0 11 544013
22249 0 11 825585
22250 0 12 444274
22251 0 12 613058
22252 0 12 513498
22253 0 12 419795
22254 0 12 432261
22255 0 12 211717
22256 0 12 432832
22257 0 12 238591
22258 0 12 345336
22259 0 12 359076
22260 0 12 275059
22261 0 12 258255
22262 0 12 161090
22263 0 12 232518
22264 0 12 284216
22265 0 12 130504
22266 0 12 189958
22267 0 12 195392
22268 0 12 173444
22269 0 12 166438
22270 0 12 212410
22271 0 12 298384
22272 0 12 234968
22273 0 12 224480
22274 0 12 125155
22275 0 12 346394
22276 0 12 352458
22277 0 12 262452
22278 0 12 400913
22279 0 12 395405
22280 0 12 374545
22281 0  1  71661
22282 0  1 385205
22283 0  1 643502
22284 0  1 643509
22285 0  1 676360
22286 0  1 483264
22287 0  1 433171
22288 0  1 451833
22289 0  1 518764
22290 0  1 563800
22291 0  1 541312
22292 0  1 591733
22293 0  1 472663
22294 0  1 277247
22295 0  1 417866
22296 0  1 605605
22297 0  1 491667
22298 0  1 455746
22299 0  1 519541
22300 0  1 537373
22301 0  1 309939
22302 0  1 280627
22303 0  1 262320
22304 0  1 317909
22305 0  1 299978
22306 0  1 382633
22307 0  1 572021
22308 0  1 308382
22309 0  1 268690
22310 0  1 395303
22311 0  1 206577
22312 0  2 151037
22313 0  2 255222
22314 0  2  80864
22315 0  2  21935
22316 0  2  17483
22317 0  2   7812
22318 0  2  27121
22319 0  2   8382
22320 0  2  17019
22321 0  2   6675
22322 0  2   5969
22323 0  2   8100
22324 0  2   7308
22326 0  2  11820
22327 0  2  18388
22328 0  2  29949
22329 0  2  28561
22330 0  2  14797
22331 0  2  37500
22367 0  3  38110
22368 0  3  39812
22369 0  3 119093
22550 1  9  30572
22551 1  9  98512
22552 1  9 182601
22553 1  9 138817
22554 1 10 674764
22555 1 10 558345
22556 1 10 634476
22557 1 10 682916
22558 1 10 789773
22559 1 10 744648
22560 1 10 664479
22561 1 10 622090
22562 1 10 547615
22563 1 10 595351
22564 1 10 702597
22565 1 10 584840
22566 1 10 554814
22567 1 10 622070
22568 1 10 634559
22569 1 10 696918
22570 1 10 723898
22571 1 10 777443
22572 1 10 734257
22573 1 10 651609
22574 1 10 512716
22575 1 10 506459
22576 1 10 553743
22577 1 10 315231
22578 1 10 366495
22579 1 10 410387
22580 1 10 286126
22581 1 10 270082
22582 1 10 474542
22583 1 10 236578
22584 1 10 195127
22585 1 11 388184
22586 1 11 469820
22587 1 11 312197
22588 1 11 217251
22589 1 11 251247
22590 1 11 370173
22591 1 11 324401
22592 1 11 345120
22593 1 11 192244
22594 1 11 357838
22595 1 11 296727
22596 1 11 369587
22597 1 11 209974
22598 1 11 232645
22599 1 11 389248
22600 1 11 297196
22601 1 11 381549
22602 1 11 259946
22603 1 11 448258
22604 1 11 374545
22605 1 11 167131
22606 1 11 148915
22607 1 11 125862
22608 1 11 121988
22609 1 11 308078
22610 1 11 296548
22611 1 11 190741
22612 1 11 228828
22613 1 11 157088
22614 1 11 145906
22615 1 12 125843
22616 1 12 142378
22617 1 12 130781
22618 1 12  98331
22619 1 12 156848
22620 1 12  80604
22621 1 12 155936
22622 1 12 144203
22623 1 12 119197
22624 1 12 136930
22625 1 12 118505
22626 1 12 147002
22627 1 12 152242
22628 1 12 184725
22629 1 12 202855
22630 1 12 249227
22631 1 12 282639
22632 1 12 147041
22633 1 12 194177
22634 1 12 225637
22635 1 12 204281
22636 1 12 222460
22637 1 12 275612
22638 1 12 324776
22639 1 12 115118
22640 1 12 298012
22641 1 12 281119
22642 1 12 324554
22643 1 12 376365
22644 1 12 542017
22645 1 12 647953
22646 1  1 182062
22647 1  1 440316
22648 1  1 590459
22649 1  1 420291
22650 1  1 577376
22651 1  1 453879
22652 1  1 395481
22653 1  1 497996
22654 1  1 459708
22655 1  1 434478
22656 1  1 568668
22657 1  1 484064
22658 1  1 366658
22659 1  1 280233
22660 1  1 393004
22661 1  1 374622
22662 1  1 548439
22663 1  1 489415
22664 1  1 499840
22665 1  1 477615
22666 1  1 455946
22667 1  1 333159
22668 1  1 304074
22669 1  1 546346
22670 1  1 336759
22671 1  1 337529
22672 1  1 222791
22673 1  1 296321
22674 1  1 136709
22675 1  1 324511
22676 1  1 182665
22677 1  2 377622
22678 1  2 171317
22679 1  2 158570
22680 1  2 209780
22681 1  2 198458
22682 1  2 109427
22683 1  2 325986
22684 1  2 310593
22685 1  2 221924
22686 1  2 304629
22687 1  2 147061
22688 1  2 165274
22689 1  2 118849
22690 1  2 230544
22691 1  2 239814
22692 1  2 160647
22693 1  2 148698
22694 1  2 230646
22695 1  2 188020
22696 1  2 372197
22697 1  2 155793
22698 1  2 225878
22699 1  2 221937
22700 1  2  97887
22701 1  2  78709
22702 1  2 143795
22703 1  2 149986
22704 1  2 161023
end
format %td date
label values year year
label def year 0 "2020-2021", modify
label def year 1 "2021-2022", modify