r - 使用 gee() : NA/NaN/Inf in foreign function call (arg 3) 拟合模型时出错

Question

我正在一个包含 13,500 个观察(这里是学生)的数据集上拟合一个 Gee 模型。学生分为 52 所不同的学校。我知道有证据表明学生嵌套在学校内(低 ICC)，因此我应该在方差协方差矩阵中调整这种嵌套效应。我打算做的是首先拟合具有可交换 var-cov 结构的 Gee 模型。然后，最重要的是，我将运行 Huber-White Sandwich 估计器，也称为稳健方差估计器。我为强大的方差估计器编写了自己的代码，它运行良好。我的 gee 语句不起作用并给出以下错误:

NA/NaN/Inf in foreign function call (arg 3)

这是我的代码:

STMath.OneYr.C1 = gee(postCSTMath1Yr ~ TRT1Yr + preCSTMath + preCSTENG + 
post1YrGradeRef + ELLBaseLine + GENDER + ECODIS + ETHNICITY.F + 
as.factor(FailedInd1Yr), data = UCI.clone[UCI.clone$COHORT0809 == "C1",], 
id =  post1YrSchIID, corstr = "exchangeable") 

不幸的是，上面的代码对于你们来说是不可重现的，也许很难弄清楚问题是什么。

如果您能帮我解决问题，我将不胜感激。

Answer 1

好的，这个问题已经很老了，但我最终来到了这里，所以这最终可能会对某人有所帮助。

基本上，这个错误是因为与其他库不同，id 引起的。参数被视为数字向量。

确实，gee功能正在类型转换 id作为一个双重，我真的不明白。以下是隐含的行(函数的 l. 119-120):

if (!(is.double(id))) 
  id <- as.double(id)

如果您的 id column 是一个字符，只需将其转换为一个因子，或者使用一些函数(如 dplyr::min_rank )将其转换为数字变量。

这应该可以解决问题。