我在 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;
任何帮助,将不胜感激