这是一个非常具体的问题。 我从R（3.2.5）Ubuntu系统的geepack安装了geeglm模型。这是我的命令
fga <- geeglm(response~covariates, data=dataframe, id=pid,
family=binomial, scale.fix=T, corstr='exchangeable')
数据不均衡。不同的pid组具有1到16个测量值。我现在想要的是提取最大组的zcor或工作相关矩阵（具有16个测量的组）。我知道这可能会失效，结果将是矢量而不是矩阵，但显示成对相关。然而，经过几个小时的搜索，我找不到办法做到这一点。我也知道如果我使用gee包中的gee，你可以得到这个矩阵。但是我们怎样才能通过使用geepack获得这个work.correlation矩阵？
欣赏任何输入！