Julia 如何在广义混合模型中加入权重参数

如何将观察的权重添加到混合模型中？ 我想我可以将Freq列添加到wt参数中，但显然不行

using RDatasets MixedModels
titanic = RDatasets.dataset("datasets", "Titanic")
titanic.surv_flg = titanic.Survived .== "Yes";

这是：

MixedModels.fit(GeneralizedLinearMixedModel, @formula(surv_flg ~ 1 + Age + Sex + (1 | Class)), titanic, Bernoulli(), nAGQ = 2, fast = true)

但事实并非如此

MixedModels.fit(GeneralizedLinearMixedModel, @formula(surv_flg ~ 1 + Age * Sex + (1 | Class)), titanic, wt = Freq, Bernoulli(), nAGQ = 2, fast = true)

---

我在另一个论坛上发现了这一点

参数应该是

wts

而不是

wt

因此，它应该是：

---

MixedModels.fit（广义线性Armixedmodel，@formula（surv_flg~1+年龄*性别+（1 |等级）），泰坦尼克号，wts=Freq，Bernoulli（），nAGQ=2）

我在另一个论坛上发现了这一点

参数应该是

wts

而不是

wt

因此，它应该是：

MixedModels.fit（广义线性Armixedmodel，@formula（surv_flg~1+年龄*性别+（1 |等级）），泰坦尼克号，wts=Freq，Bernoulli（），nAGQ=2）