[statnet_help] Question about Model Stuck in MCMLE

[Carter T. Butts]

If you are using distance as a predictor, you are positing that tie 
probability falls exponentially (approximately) in distance; this is 
probably too sharp an SIF. Try using the log distance instead.

-Carter

On 4/9/23 6:30 AM, 陈博川 wrote:

> Statnet Developer:

> Good morning. Thank you so much for your patience to look into my 

> question. I am encountering some questions beyond my capability so I 

> have no choice but to turn to you for help.

>

> I am encountering the problem of MCMLE stuck. My original code looked like

>

> gw1_help <- ergm(helpnet ~ edges + nodefactor("work") + 

> nodefactor("income") + nodefactor("edu") + nodefactor("religious") + 

> nodefactor("gender") + nodecov("age") + edgecov(kin_sup,"kinweight") + 

> gwidegree(.1, T) + gwesp(.1, T) + gwdsp(.1, T), control = 

> control.ergm(MCMC.samplesize = 1e+5, MCMC.burnin = 1e+6, MCMC.interval 

> = 1000, seed = 567), eval.loglik = T, verbose = T)

>

> which worked qiuite well. However, once I add 'edgecov(sup_loca2_net, 

> "distance") ' in, which looked like:

>

> gw1_help <- ergm(helpnet ~ edges + nodefactor("work") + 

> nodefactor("income") + nodefactor("edu") + nodefactor("religious") + 

> nodefactor("gender") + nodecov("age") + edgecov(sup_loca2_net, 

> "distance") + edgecov(kin_sup,"kinweight") + gwidegree(.1, T) + 

> gwesp(.1, T) + gwdsp(.1, T), control = control.ergm(MCMC.samplesize = 

> 1e+5, MCMC.burnin = 1e+6, MCMC.interval = 1000, seed = 567), 

> eval.loglik = T, verbose = T)

>

> the model stuck and said "*_MCMLE estimation stuck. There may be 

> excessive correlation between model terms, suggesting a poor model for 

> the observed data. If target.stats are specified, try increasing SAN 

> _**param*eter."

>

> I don't know what happened and I desperately want to know how to solve 

> this problem, Thanks for your time and patience! Hope to hear from you 

> soon!

>

>

> _______________________________________________

> statnet_help mailing list

> statnet_help at u.washington.edu

> https://urldefense.com/v3/__http://mailman13.u.washington.edu/mailman/listinfo/statnet_help__;!!CzAuKJ42GuquVTTmVmPViYEvSg!MXKxthAW4qsCYVUIoNakj72_N0d2pl2p1VQqNWKwlCEVOKwOY_y_R8rqy2SuSzxRcojQkfOTzeaCjimCM1NbTpS-EGs$  

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman13.u.washington.edu/pipermail/statnet_help/attachments/20230409/ddf0d365/attachment.html>