我有下面的数据集
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
我的因变量是估值，我想建立一个线性混合效应模型，其中school_id是我的随机效应。但我希望该模型因主题和分数水平而有所不同。我想用python建立这个。
如果有比线性混合效应更好的模型，谁能建议我也应该在python中这样做吗？