0

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

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

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

data example;
   input price cat time x2 x3;
   cards; 
   5000 1 1 5.4  50
   6000 1 2 6    45
   3000 1 3 7     60
   4000 2 1 5    50
   4500 2 2 5.4  75
   4786 3 1 6    33
   6500 3 2 5.8  36
   1010 3 3 4    41
;;;;
run;


proc glm data=example PLOTS(UNPACK)=DIAGNOSTIC;
    class cat;
   model price= cat time x2 x3/ noint solution;
run;

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

ありがとう

4

1 に答える 1

0

出力をファイルに追加する必要があります: (意図的にスクリーン印刷を無効にしましたが、必要に応じて自由に有効にしてください。)

proc glm data=example noprint;
    class cat;
    model price= cat time x2 x3  / noint solution ;
   **output out= from_glm COVRATIO = Cov ;**
run; quit;

結果:

price cat time x2   x3 COV 
5000  1   1    5.4  50 597.2565 
6000  1   2    6    45 8.312725 
3000  1   3    7    60 0.0493 
....

編集：出力ステートメントを更新しました。

キーワードの詳細については、https: //support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_glm_sect020.htm を参照してください。

うまくいけば、これがあなたが求めているものです。

于 2017-09-04T08:16:04.583 に答える