如何从gamlss提取多分类混合效应模型的k-1组回归系数？

gamlss针对多分类MN族（如MN4）确实会生成k-1组回归系数，默认输出仅展示基准组外第一组的结果，其余系数需要通过指定参数或直接访问模型对象提取：

• 提取固定效应系数：使用coef()函数，指定what = "mu"并开启print.all = TRUE，即可输出所有非基准水平的系数组：

  # 先将模型赋值给变量
  model_obj <- gamlss::gamlss(
    formula = y ~ wk + gamlss::random(factor(trt)),
    sigma.formula = ~1,
    family = gamlss.dist::MN4,
    data = schizo
  )
  # 提取所有mu系数
  coef(model_obj, what = "mu", print.all = TRUE)
  

• 直接访问模型内部存储：多分类模型的mu系数以列表形式存储在model_obj$mu.coefficients中，每个列表元素对应一组非基准水平的系数。

• 提取随机效应：使用ranef()函数查看分组随机效应的估计值：

  ranef(model_obj)
  

注：MN4族默认将因变量的第一个水平作为基准组，所有系数均为该组相对于基准组的效应估计。

内容的提问来源于stack exchange，提问作者bison2178