Hi everyone, I'm getting an error code (no specific one, just following note: Stopped by error!) when trying to run following Proc Genmod procedure with an unstructured working correlation:
PROC Genmod DATA= table descending;
Class ID Year (ref='1');
Model TreatmentType=Year/ Dist=Binomial Link=Logit;
Repeated Subject=VID/ Type=UN corrw;
The working correlation matrix in the output looks like this (looks like it stopped running):
Col1 Col2 Col3 Col4 Col5 Col6
Row1 2.9295 -1.0000 0.0000 0.0000 0.0000 0.0000
Row2 -1.0000 0.3414 0.0000 0.0000 0.0000 0.0000
Row3 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Row4 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Row5 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Row6 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
However, when I change the working correlation to exchangeable the procedure will run.
How does this come?
Thanks in advance!
The unstructured matrix requires the estimation of most correlation parameters - 15 in your case - as opposed to much simpler structures, like exchangeable, which only requires one to be estimated. As with any iterative model fitting algorithm, increasing the number of parameters to estimate, either model parameters or correlation parameters, increases the likelihood of fitting problems. This is particularly true for binary response models where data sparseness can make these problems more likely. That said, the more modern PROC GEE is the recommended procedure to fit GEE models rather than the older PROC GENMOD. It is possible, though by no means certain, that you'll have better luck there.
I use MIXED and GLIMMIX much more than I do GENMOD. But.... I think you should make VID a class variable (add it to the CLASS statement). Also, if you have more than one observation for each combination of year and VID, you will have a problem.
Good news: We've extended SAS Hackathon registration until Sept. 12, so you still have time to be part of our biggest event yet – our five-year anniversary!
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.