我正在估计 R 中的广义估计方程。我有一个多项式（序数）结果，所以我一直在尝试使用这个包multgee，因为据我所知，包喜欢geepack或gee不允许估计多项式结果。
但是，我遇到了一些问题。文档看起来不错，但特别是，它似乎需要模型中的初始（起始）值。如果我尝试在没有它的情况下运行模型，我会得到一条请求起始值的行。这是模型：
formula <- PAINAD_recode~Age + Gender + group_2part + Cornell + stim_intensity_scale
fitmod <- ordLORgee(formula,data=data,bstart=c(1,0,1,0,1),id=data$subject,
repeated=data$trial)
我只是为那里的起始值添加了一些 1 和 0。但是，当我输入起始值（甚至是合理的值）时，它声称：
起始值和参数向量的长度不同
我认为有五个预测变量，我需要五个起始值。我找不到有关此特定矩阵的更多信息。有人对此有任何想法吗？这里的结果有五个级别（序数），重复分量有 20 个级别。任何建议，将不胜感激。