It's the first time I'm posting here, and I apologize in advance if the same question has been answered previously.
I am working with data of proportions, and I want to know how that proportion changes over time. For example, below are random telephone survey for what middle schooler in a particular state do with their lunch at different years. I want to know what would be the optimal statistical method to analyze the change in proportion of each option across time.
State A
2000 total: 200 Own lunch:120 Buy lunch:45 School lunch:73 Other:24
2001 total: 220 Own lunch:89 Buy lunch:39 School lunch:85 Other:17
2002 total: 270 Own lunch:70 Buy lunch:52 School lunch:121 Other:32
2003 total: 354 Own lunch:45 Buy lunch:71 School lunch:224 Other:37
2004 total: 241 Own lunch:32 Buy lunch:37 School lunch:168 Other:29
2005 total: 530 Own lunch:51 Buy lunch:94 School lunch:403 Other:14
and so on....
On the basic level, I know I can compare, say, proportion of students who bring their own lunch between 2000 and 2001 with fisher's exact or chi-square. But I want to know if there is a way to analyze changes in all options (as they are not mutually exclusive, students can pick multiple options so one option going up doesn't necessarily mean other options have to go down) and across quantitative time (instead of just two categorical time points). In addition, I have this data for many states, so it would so be nice if I can incorporate that into the same model to see if the change in proportion across time differs by states.
Please help me identify the optimal test/analysis method, ideally something I can run on R.