Programming the statistical procedures from SAS

Glimmix did not converge or covariance matrix doesn't make sense

Reply
New Contributor
Posts: 2

Glimmix did not converge or covariance matrix doesn't make sense

Hello,

 

I have several clinical data set with binary response. For all data sets, subjects are enrolled from different sites and are taking either treatment A or B (parallel design). I am trying to assess the site variation for the response variable. Since there are 2 treatment, I assume that the site variation could be different for the 2 treatments and there could be a correlation. The following is my SAS code:

 

proc glimmix data=dset;
class sitenum trt;
model resp= trt/dist=bin link=logit solution;
random trt/SUBJECT=SITENUM type=UN solution G;
run;

 

For most data sets, Proc Glimmix doesn't converge. For one data set, Glimmix does converge and I got the following covariance matrix. I think here 0.00907 is the site variance of treatment A and 3.54E-18 is the site variance of treatment B. However, the covariance between treatments A and B are even bigger than the both variance, so the correlation will be >1. Do I understand it correct? I don't know why it can happened.

 

Estimated G Matrix
EffectDescription of Actual ArmRowCol1Col2
trtA10.009070.02008
trtB20.020083.54E-18
Super User
Posts: 13,523

Re: Glimmix did not converge or covariance matrix doesn't make sense

Posted in reply to Fang_Merck

I am a little confused by:

For most data sets, Proc Glimmix doesn't converge

 

I would expect that you would combine all of the data to have all of the sites in one data set to answer a question about the impact of "site" on the analysis. Why are these separate data sets analyzed one a time? Are you comparing all pairs of individual sites?

New Contributor
Posts: 2

Re: Glimmix did not converge or covariance matrix doesn't make sense

I meant that I have several independent clinical trial data set for
different drug. That is why I did the analysis separately.


Ask a Question
Discussion stats
  • 2 replies
  • 160 views
  • 1 like
  • 2 in conversation