BookmarkSubscribeRSS Feed
Jerem
Calcite | Level 5

Many thanks to both of you!

I've been extremely busy over the last days so I couldn't work that much on the model. I've thus need to have a closer look at this later.

Now the model do not converge....I tried different dependent variable and this is always the same problem: "did not converge".

Is it a problem in the modelling or the data?

Best,

Jérémy

SteveDenham
Jade | Level 19

"Did not converge" isn't the same as "Will not noverge".  We can address the first with the NLOPTIONS statement, and by adjusting the convergence criterion in the PROC GLIMMIX statement.  In order to be of assistance, though, I need to see code and output for that code, especially the iteration history

Steve Denham

Jerem
Calcite | Level 5

Hi Steve,

I included the code as well as the SAS output. Is that all you need?

I'm really grateful for your help here!

proc glimmix data=candidates;

nloptions maxiter=100;

  class commune ;

  model  binouverture (event='1') = EurostatA EurostatC sizelog ENP commune

      / solution dist=binary link=logit ddfm=kenwardrogers ; /* Addresses individual commune levels, and adds the Kenward-Rogers adjustment for repeated measures designs.  This also applies the averaged effects of the continuous covariates across all communes */

  random intercept / subject= commune solution residual; /* Converts to a repeated measures design.  Even though this is a RANDOM statement, the residual option means treat the values within a commune as correlated.  */

lsmeans commune/ilink cl; /* Outputs the estimates for each level of commune, on both the logit and original scale */

run;

Jérémy

Le Système SAS

The GLIMMIX Procedure

Data SetWORK.CANDIDATES
Response Variablebinouverture
Response DistributionBinary
Link FunctionLogit
Variance FunctionDefault
Variance Matrix Blocked ByCommune
Estimation TechniqueResidual PL
Degrees of Freedom MethodKenward-Roger
Fixed Effects SE AdjustmentPrasad-Rao-Kackar-Harville-Jeske

Commune252Aiseau-Presles Amay Amel Andenne Anderlues Anhée Ans Anthisnes Antoing Arlon Assesse Ath Attert Aubange Aubel Awans Aywaille Baelen Bassange Bastogne Beaumont Beauraing Beauvechain Beloeil Berloz Bernissart Bertogne Bertrix Beyne-Heusay Binche Bièvre Blegny Bouillon Boussu Braine-l'Alleud Braine-le-Château Braine-le-Comte Braives Brugelette Brunehaut Burdinne Burg-Reuland Büllingen Bütgenbach Celles Cerfontaine Chapelle-Lez-Herlaimont Charleroi Chastre Chaudfontaine Chaumont-Gistoux Chimay Chiny Chièvres Châtelet Ciney Colfontaine Comblain-au-Pont Comines-Warneton Courcelles Court-St-Etienne Couvin Crisnée Dalhem Dinant Dison Doische Donceel Dour Durbuy Ecaussinnes Eghezée Ellezelles Enghien Engis Erezée Erquelinnes Esneux Estaimpuis Estinnes Etalle Eupen Faimes Farciennes Fauvillers Fernelmont Ferrières Fexhe-le-Haut-Clocher Flemalle Fleurus Flobecq Floreffe Florennes Florenville Fléron Fontaine-l'Evêque Fosses-la-Ville Frameries Frasnes-lez-Anvaing Froidchapelle Gedinne Geer Gembloux Genappe Gerpinnes Gesves Gouvy Grace-Hollogne Grez-Doiceau Habay Ham-sur-Heure-Nalinnes Hamois Hannut Hastière Havelange Hensies Herbeumont Herstal Herve Honnelles Hotton Houffalize Houyet Huy Hélecine Héron Incourt Ittre Jalhay Jemeppe-sur-Sambre Jodoigne Juprelle Jurbise Kelmis La Bruyère La Hulpe La Louvière La Roche-en-Ardenne Lasne Le Roeulx Lens Les-Bons-Villers Lessines Leuze-en-Hainaut Libramont-Chevigny Lierneux Limbourg Lincent Liège Lobbes Lontzen Malmedy Manage Manhay Marche-en-Famenne Marchin Martelange Meix-devant-Virton Merbes-le-Château Messancy Mettet Modave Momignies Mons Mont-St-Guibert Mont-de-l'Enclus Montigny-le-Tilleul Morlanwelz Mouscron Musson Namur Nandrin Nassogne Neufchâteau Neupré Nivelles Ohey Olne Onhaye Oreye Orp-Jauche Ottignies-LLN Ouffet Oupeye Paliseul Pecq Pepinster Perwez Philippeville Plombières Pont-à-Celles Profondeville Péruwelz Quaregnon Quiévrain Quévy Raeren Ramilies Rebecq Remicourt Rixensart Rochefort Rouvroy Rumes Saint Nicolas Saint-Georges-sur-Meuse Saint-Ghislain Saint-Hubert Sainte-Ode Sambreville Sankt Vith Seneffe Seraing Silly Sivry-Rance Soignies Sombreffe Somme-Leuze Soumagne Spa Sprimont Stavelot Stoumont Tellin Theux Thimister-Clermont Thuin Tinlot Tintigny Tournai Trooz Tubize Vaux-sur-Sûre Verlaine Verviers Vielsalm Villers-la-Ville Villers-le-Bouillet Viroinval Virton Visé Vresse-sur-Semois Walcourt Walhain Wanze Waremme Wasseiges Waterloo Wavre Welkenraedt Wellin Yvoir

Number of Observations Read804
Number of Observations Used803

10476
21327

R-side Cov. Parameters1
Columns in X257
Columns in Z per Subject0
Subjects (Blocks in V)252
Max Obs per Subject14

Optimization TechniqueNone
Parameters0
Lower Boundaries0
Upper Boundaries0
Fixed EffectsProfiled
Residual VarianceProfiled
Starting FromData

0001.340781E1542.00000000.
1001.340781E1541.13896683.
2001.340781E1542.00000000.
3001.340781E1542.00000000.
4001.340781E1540.91900698.
5003460.72436132.00000000.
6001.340781E1542.00000000.
7003728.57513632.00000000.
8001.340781E1542.00000000.
9003993.86542062.00000000.
10004127.23543912.00000000.
11001.340781E1542.00000000.
12004394.11295972.00000000.
13004525.23404232.00000000.
14004658.35009272.00000000.
15004791.71562.00000000.
16004924.19137272.00000000.
17005059.80102772.00000000.
18005191.39611152.00000000.
19005324.18239922.00000000.
20001.340781E1542.00000000.
21001.340781E1542.00000000.
22001.340781E1542.00000000.
23001.340781E1542.00000000.
24001.340781E1542.00000000.
25001.340781E1540.99810273.
26001.340781E1540.81393176.
27001.340781E1542.00000000.
28003148.70888112.00000000.
29003281.26246662.00000000.
30001.340781E1542.00000000.
31003547.66602212.00000000.
32003677.98464132.00000000.
33003814.07183042.00000000.
34001.340781E1542.00000000.
35001.340781E1540.14868120.
36004210.50722542.00000000.
37004346.00036962.00000000.
38001.340781E1542.00000000.
39001.340781E1540.09425822.
40001.340781E1540.09072392.
41004877.31974272.00000000.
42005009.61826142.00000000.
43001.340781E1542.00000000.
44001.340781E1540.06120650.
45005409.79126142.00000000.
46001.340781E1542.00000000.
47001.340781E1542.00000000.
480016935.2892472.00000000.
49001.340781E1542.00000000.
50001.340781E1542.00000000.
51001.340781E1540.99810273.
52001.340781E1540.81393176.
53001.340781E1542.00000000.
54003148.70888112.00000000.
55003281.26246662.00000000.
56001.340781E1542.00000000.
57003547.66602212.00000000.
58003677.98464132.00000000.
59003814.07183042.00000000.
60001.340781E1542.00000000.
61001.340781E1540.14868120.
62004210.50722542.00000000.
63004346.00036962.00000000.
64001.340781E1542.00000000.
65001.340781E1540.09425822.
66001.340781E1540.09072392.
67004877.31974272.00000000.
68005009.61826142.00000000.
69001.340781E1542.00000000.
70001.340781E1540.06120650.
71005409.79126142.00000000.
72001.340781E1542.00000000.
73001.340781E1542.00000000.
740016935.2892472.00000000.
75001.340781E1542.00000000.
76001.340781E1542.00000000.
77001.340781E1540.99810273.
78001.340781E1540.81393176.
79001.340781E1542.00000000.
80003148.70888112.00000000.
81003281.26246662.00000000.
82001.340781E1542.00000000.
83003547.66602212.00000000.
84003677.98464132.00000000.
85003814.07183042.00000000.
86001.340781E1542.00000000.
87001.340781E1540.14868120.
88004210.50722542.00000000.
89004346.00036962.00000000.
90001.340781E1542.00000000.
91001.340781E1540.09425822.
92001.340781E1540.09072392.
93004877.31974272.00000000.
94005009.61826142.00000000.
95001.340781E1542.00000000.
96001.340781E1540.06120650.
97005409.79126142.00000000.
98001.340781E1542.00000000.
99001.340781E1542.00000000.

Did not converge.

Residual (VC)1.0328.

SteveDenham
Jade | Level 19

I may have spoken too hastily.  This has all the symptoms of quasi-separation--for at least some of the communes, all of the values are either 1 or zero, and all of the continuous covariates are identical.  I think you will either have to relinquish the hope of finding marginal estimates (i.e., consider commune to be a random effect), or find a way to collapse the communes into some larger categorization that eliminates the quasi-separation.  There is a major problem in trying to estimate 258 parameters from only 803 observations.  A rule of thumb for getting decent estimates is at least 5, and preferably 10, observations per parameter.  I think you may be able to build a respectable model if you could collapse the commune variable into 70 or 80 "super-commune" levels.

Steve Denham

SAS Innovate 2025: Register Now

Registration is now open for SAS Innovate 2025 , our biggest and most exciting global event of the year! Join us in Orlando, FL, May 6-9.
Sign up by Dec. 31 to get the 2024 rate of just $495.
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
  • 18 replies
  • 5286 views
  • 0 likes
  • 3 in conversation