問題タブ [longitudinal]

Proc GLM を使用して基本的な固定効果モデルに適合させており、分散/共分散行列を取得したいと考えています。モデルを proc reg に適合させる場合、これは非常に簡単なことですが、私が適合させているモデルには、クラスの各メンバー (クラスのメンバーが 50 人以上) ごとに個別の勾配があるため、そうしたくありません。それらすべてにダミー変数をコーディングします。

proc glm を使用して適合から分散共分散行列を取得する方法はありますか?

これは、構成されたデータと私のコードの例です。推定値の分散共分散行列を取得したいと思います。

各カテゴリのパラメーター推定値を取得し (これらは本質的に厄介なパラメーターです)、推定時間 x2 および x3 の共分散行列に関心があります。

ありがとう

次の長いデータセットがあります (「dat」と呼ばれる R データフレーム)。1200 人の子供の約 10,000 の観測があり、子供ごとに最大 10 の観測があります。これらの観察結果は、出生から 5 歳までの幼児期のさまざまな時点で収集された子供の体重です。データセットの年齢は日単位 (変数 = agedays) で、体重は変数 wtkg のキログラム単位で与えられます。

データセット

SITAR パッケージを使用して、肥満の母親 (変数 = BMI_group) の間で (変数 = GWG_T3_cat2) (変数 = GWG_T3_cat2) 妊娠 3 期の過度の妊娠中の体重増加に応じて、子供の速度、テンポ、およびサイズを調査したいと考えています。 .

このリンクにあるモデルを実行してみました: https://www.rdocumentation.org/packages/sitar/versions/1.0.9

しかし、私はエラーが発生します:

誰かがこの問題の解決に役立つことができれば、本当に感謝しています。