python混合线性模型_Python中的多元线性混合效应模型

我在玩这个代码，这是一元线性混合效应建模。数据集表示：学生作为s

讲师(d)

部门为部门

服务即服务

在R的lme4包语法中(Bates等人，2015)，实现的模型可以概括为：y ~ 1 + (1|students) + (1|instructor) + (1|dept) + service

其中1表示截距项，(1 | x)表示x的随机效应，x表示固定效应。在

^{pr2}$

s_train = train['s'].values

d_train = train['dcodes'].values

dept_train = train['deptcodes'].values

y_train = train['y'].values

service_train = train['service'].values

n_obs_train = train.shape[0]

s_test = test['s'].values

d_test = test['dcodes'].values

dept_test = test['deptcodes'].values

y_test = test['y'