[statnet_help] Question about Model Stuck in MCMLE
陈博川
chenbc20 at mails.tsinghua.edu.cn
Sun Apr 9 06:30:58 PDT 2023
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 parameter."
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!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman13.u.washington.edu/pipermail/statnet_help/attachments/20230409/9dfa97a3/attachment.html>
More information about the statnet_help
mailing list