I see some values in September, so are you sure that your calender runs from 1Oct - 31Mar? Otherwise, you will still have gaps with some dates in September as there are no observations between Apr-Aug. I will drop these in my illustration below. Below, I use labmask from the Stata Journal.
Res.:
Code:
net install gr0034, from(http://www.stata-journal.com/software/sj8-2)
Code:
* 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 gen doy= doy(date) gen time= doy if inrange(month(date), 10, 12) bys year (date): replace time= cond(doy==1, time[_n-1]+1, time[_n-1]+ doy - doy[_n-1]) if missing(time) drop if missing(time) bys year (time): gen time1= time[1] replace time = time-time1+1 gen mask= string(date, "%tdDm") labmask time, values(mask) xtset year time set scheme s1color xtline kilos, overlay xlab(1(30)181,val) xtitle("") leg( order(1 "2000-2001" 2 "2001-2002") col(1) pos(1) ring(0))
Comment