Linear mixed effect model that differs by variable factor level

I have below dataset

student_id  school_id   non_verbal  quantitative    spatial   verbal    subject score_level valuation
69044       198        117.0         116.0          112.000000   105.0  26        1          77.0
57589       119        104.0         96.0            101.554674  89.0   16        1           62.0
65452       195        83.0          111.0           84.000000   97.0   15        1           62.0
26203       107        98.0          109.0           107.000000  102.0  11        1             47.0
25501       17         111.0         90.0            113.000000  118.0  27        1              62.0

My dependent variable is valuation and I want to build a linear mixed effect model where school_id is my random effect. But I want this model to differ by subject and score level. I want to build this in python.

Can anyone please suggest how should I do this in python also if any better model than linear mixed effect?

Any help will be appreciated.