在R中修改混合效应公式

我有一个公式:

my_formula <- a ~ b + (1|d) + e

e

我对固定效应公式所做的是:

drop_es <- function(f) {
    e_idxs <- grep(pattern = 'e',
                   x = all.vars(f[[3]]))

    stats::formula(stats::drop.terms(termobj = stats::terms(f),
                                     dropx = e_idxs,
                                     keep.response = TRUE))
}

这只适用于固定效果公式:

drop_es(f = a ~ b + e)

drop_es(f = a ~ b + (1|d) + e)

a~b+1 | d

a ~ b|d

如何删除 e