SAS に次のコードがあり、R で同じデータ セットに対してこのモデルを再作成する必要があります。モデルを再作成するための適切なコードを見つけるのに苦労しており、特にデータの分布をガンマ分布に変更しています。
SAS コード
proc glimmix data=combined plot=(residualpanel studentpanel);
class block pen hour daymod;
model mean=hour|pen/solution ddfm=kr dist=gamma;
random _residual_/ subject=steer type=cs;
random block;
random intercept /subject=pen;
lsmeans hour/cl plots=meanplot;
output out=FinalModel pred=p resid=r lcl=lowlimit ucl=uplimit;
run;
どんな助けでも大歓迎です