確率変数 ID の PS、TH、HW の 3 つの変数の固定効果を調べたいと思います。また、固定効果ごとに各被験者 (ID) のランダム勾配の推定値を知りたいと思っています。一部の被験者は複数年のデータを持っていたので、ID 内で年をネストする必要があります。これまでのところ、予備的な分析のために次のコードがあります。
proc glimmix data=dat method=laplace;
class year;
class ID;
model use = PS TH HW / dist=binomial link=logit solution;
random intercept / subject = year(ID);
random PS / subject = year(ID) solution;
random TH / subject = year(ID) solution;
random HW / subject = year(ID) solution;
run;
このコードを実行すると、固定効果ごとに year(ID) のランダムな勾配推定値が得られます。
Solution for Random Effects
Std Err
Effect Subject Estimate Pred
PS year (ID) 2 111 1.1714 0.7257
TH year (ID) 2 111 4.5081 1.2953
HW year (ID) 2 111 -1.3109 0.4196
PS year (ID) 4 114 0.5999 0.5608
TH year (ID) 4 114 4.7292 1.0664
HW year (ID) 4 114 -1.8526 0.3289
それは良いことですが、 IDだけのランダムな勾配の推定値も取得したいと思います。理想的には、モデルは1 年以上の ID を考慮し、年 (ID)とID のランダムな勾配推定値を提供します。
質問: SAS に ID のみのソリューションを出力するように指示するにはどうすればよいですか?
それを行うために次のコードを実行しようとしたところ、ランダム効果の出力に奇妙な空白が表示され、パラメーターの推定値が上記のコードで取得した出力とは異なりました。以下のコードは複雑すぎるようです。
proc glimmix data=dat method=laplace;
class year;
class ID;
model use = PS TH HW / dist=binomial link=logit solution;
random intercept / subject = ID;
random PS / subject = ID solution;
random TH / subject = ID solution;
random HW / subject = ID solution;
random intercept / subject = year(ID);
random PS / subject = year(ID) solution;
random TH / subject = year(ID) solution;
random HW / subject = year(ID) solution;
run;
Example of part of the "Solutions for Random Effects" table
Effect Subject Estimate Pred
PS ID 111 0.5587 1.4640
TH ID 111 4.2063 2.1891
HW ID 111 -0.9158 0.6462
PS Year(ID) 2 111 0.4944 1.4721
PS 0 1.8627
TH Year(ID) 2 111 0.3569 1.7628
TH 0 1.8037
HW Year(ID) 2 111 -0.3687 0.5800
HW 0 0.5639
PS ID 114 0.2531 1.3999
TH ID 114 4.4162 2.0840
HW ID 114 -1.2981 0.6914
PS Year(ID) 4 114 0.2240 1.4027
*注意 - 件名なしでリストされている各固定効果内のすべての「pred」値は同じです。