求助：中介变量及因变量都是分类变量的中介效应分析，其中因变量为无序三分类变量

以下是我的数据格式。其中，a是因变量，为无序三分类变量；b是中介变量，为二分类变量；c是自变量，为连续变量；d是控制变量。请问我该如何实现中介效应分析呢？


<

a	b	c	d
1	0	3	1
1	0	4.8	0
1	1	2.95	0
1	0	4.8	0
0	1	3.7	0
0	0	3	0
1	0	1.9	0
1	1	3.3	0
0	1	2.75	1
1	1	2.6	0
0	0	4.6	0
0	0	2.8	0
0	0	2.1	0
1	0	2.15	1
1	0	4.3	0
0	1	3	0
1	0	3.45	0
2	0	2.1	1
1	1	2.9	1
1	0	2.9	0
0	0	4.8	0
1	0	2.15	0
0	1	2.5	1
0	0	3.9	0
0	0	4.5	0
1	0	4	0
1	1	2.3	1
1	1	2.3	0
0	0	4.12	0
0	1	3.54	0
1	0	2.5	0
1	0	2.1	1
1	0	2.5	0
1	0	3.75	1
1	0	2.5	0
0	0	3	0
1	1	2.02	0
1	0	3.75	0
1	0	2.31	1
1	0	2.55	0

这是一个模拟数据；变量a实际不会只有1个取值2，敬请提供解决思路~