Hello stata community,

I want to graph the total number of jobs available yearly based on SOC ( Standard Occupation Characteristics ). Here, I've given the sample of data which provides the year and SOC 2 variable which says which industry the person is in ( soc2=11 if // Management, Business, Science, and Arts, soc2=13 if Business Operations Specialists, soc2=23 if legal, soc2=53 if Transportation and Material Moving and so on ). Would like to produce a graph like the following one where the x variable will be a year and the y variable is the total number of jobs available in each industry ( when soc2= 52 , soc2= 49 , soc2=11, and soc2 =13 or any number. )

Could anyone kindly give me some suggestions on that?

input double cpsidp float soc2 int(occ90 year) byte month
20081000084401 . 999 2010 1
20081000084401 . 999 2010 1
20081000084401 . 999 2010 1
20081000085001 . 999 2010 1
20081000085001 . 999 2010 1
20081000085001 . 999 2010 1
20081000085002 53 889 2011 1
20081000085002 53 889 2011 1
20081000085002 53 889 2011 1
20081000088401 . 999 2011 1
20081000088401 . 999 2011 1
20081000088401 . 999 2011 1
20081000089801 . 999 2011 1
20081000089801 . 999 2012 1
20081000089801 . 999 2012 1
20081000089802 49 518 2012 1
20081000089802 49 518 2012 1
20081000089802 49 518 2012 1
20081000097001 35 436 2012 1
20081000097001 35 436 2013 1
20081000097001 35 436 2013 1
20081000097002 45 498 2013 1
20081000097002 45 498 2013 1
20081000097002 45 498 2013 1
20081000097101 39 458 2014 1
20081000097101 39 458 2014 1
20081000097101 39 458 2014 1
20081000098201 37 453 2014 1
20081000098201 37 453 2014 1
20081000098201 37 453 2014 1
20081000098202 43 313 2014 1
20081000098202 43 313 2014 1