13,500 の観測データ (ここでは学生) を含むデータセットにジー モデルを当てはめています。生徒は 52 の異なる学校にグループ分けされます。学生が学校内に入れ子になっている (ICC が低い) という証拠があることはわかっているため、分散共分散行列でこの入れ子効果を調整する必要があります。私が計画しているのは、交換可能な var-cov 構造を持つ gee モデルを最初に適合させることです。その上で、ロバスト分散推定器とも呼ばれる Huber-White Sandwich 推定器を実行します。堅牢な分散推定器用に独自のコードを作成しましたが、完全に機能します。私のジーステートメントは機能せず、以下のエラーが表示されます:
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")
残念ながら、上記のコードは再現性がなく、問題が何であるかを理解するのが難しいかもしれません.
問題を解決する方法を教えていただければ幸いです。