BookmarkSubscribeRSS Feed
raisins
Calcite | Level 5
Hi,

I'm fitting a generalized estimating equation-type model by specifying the empirical model option in proc glimmix. I am using this in order to model spatial autocorrelation among residuals. The model response is 0/1 and is modeled with a binomial distribution. I know that proc genmod outputs Pan's QIC (quasilikelihood version of AIC) but proc genmod does not (to my knowledge) include the spatial modeling options of glimmix. Has anyone found a way to extract or calculate QIC based upon glimmix model output? I could potentially "borrow" the fitted correlation structure from glimmix and use it as a user-defined correlation structure in genmod, and then get QIC from genmod. However, my data set is quite large with >100 spatial clusters and so handling this seems it would be quite tricky. The empirical option is being used because without it my machine runs out of memory in its efforts to populate a huge matrix that includes the correlation among all data points, even the very distant. My data is quite "cheap" and I have discarded many points (it is satellite imagery). So, for now I am simply validating with an independent data set for model comparison purposes. I believe this gets at the purpose of AIC, but you know how some reviewers are with their insistence upon including the latest buzz word they've latched ahold of. And for a few, it seems "information theoretic" is it. Any guidance is appreciated.-Seth
2 REPLIES 2
Dale
Pyrite | Level 9
Before there was support for the QIC statistic in the GENMOD procecdure, SAS Institute offered a macro for computing QIC. You might see if you can modify that macro to produce the QIC statistic when the GLIMMIX procedure is employed to estimate model parameters. See:

http://support.sas.com/kb/26/100.html
http://support.sas.com/kb/26/addl/fusion26100_1_qic.sas.txt

The first link provides a description of the macro and demonstrates use of the macro (under the Results tab). The second link provides the macro code.

HTH
bno
Calcite | Level 5 bno
Calcite | Level 5

Hi,

I have similar problem. I have over 300 clusters, where sizes of clusters varied from 1 (there are many totally independent subjects) to 250 of individuals with total of over 5000 individuals. I have already obtained correlation matrix between indiviuals using R for each of the clusters. I would like to used PROC GENMOD in SAS and use my own correlation structure instead of using say, EXCH for each cluster. There is an option for userdefined matrix in SAS for a single matrix, but how would I specify these several matrices for these clusters? Also, how do I make sure that the SAS picks up the correct correlation between any two individuals from my data in the GENMOD analysis?

Is there any way to specify user defined correlation matrix of all individuals in GENMOD without spefifying clustering variable?

Many thanks in advance.

BNO.

sas-innovate-2024.png

Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.

Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.

 

Register now!

What is ANOVA?

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.

Discussion stats
  • 2 replies
  • 1827 views
  • 0 likes
  • 3 in conversation