1

我在 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;

任何帮助,将不胜感激

4

0 回答 0