Negative Binomial model is crashing
Posted: Wed Feb 20, 2019 4:46 pm
Hi it's me again. I need to run different models for different types of outcomes of social support, and now I tried a Poisson model and it went fine, but the Neg binomial is constantly crashing.
This is the formula.
#######################################################
m8 = log(AltSupportCount) ~ 1 + TieIsFamily + TieSameSex + TieSameAge + EgoMarital + EgoSex +
EgoTMig + AltSex + AltIsRoma + AltNatRomanian + AltResidence + NetPropFemCen + NetPropFamCen +
NetPropIsRomaCen + NetPropNatRomanianCen + NetPropResidCen + EgoTMig:TieSameSex + EgoSex:AltIsRoma +
EgoTMig:AltIsRoma + (1 | EgoIdentifier)
(VarCompModel <- runMLwiN(Formula = m8, D = "Negbinom", data = RomaAllData,
estoptions = list(EstM = 1)))
This is the error on the console
"error while obeying batch file /var/folders/3n/zvw9lr5s1c1dnl0jkrc8y2xh0000gn/T//Rtmp1LfSTP/macrofile_9f41447b133.txt at line number 195:
MCMC 0 500 1 5.8 50 10 C2498] C2499] 1 1 1 1 1 8
MCMC Error 0010 : Wrong parameters for Gamma distribution generator
Execution completed"
Best,
This is the formula.
#######################################################
m8 = log(AltSupportCount) ~ 1 + TieIsFamily + TieSameSex + TieSameAge + EgoMarital + EgoSex +
EgoTMig + AltSex + AltIsRoma + AltNatRomanian + AltResidence + NetPropFemCen + NetPropFamCen +
NetPropIsRomaCen + NetPropNatRomanianCen + NetPropResidCen + EgoTMig:TieSameSex + EgoSex:AltIsRoma +
EgoTMig:AltIsRoma + (1 | EgoIdentifier)
(VarCompModel <- runMLwiN(Formula = m8, D = "Negbinom", data = RomaAllData,
estoptions = list(EstM = 1)))
This is the error on the console
"error while obeying batch file /var/folders/3n/zvw9lr5s1c1dnl0jkrc8y2xh0000gn/T//Rtmp1LfSTP/macrofile_9f41447b133.txt at line number 195:
MCMC 0 500 1 5.8 50 10 C2498] C2499] 1 1 1 1 1 8
MCMC Error 0010 : Wrong parameters for Gamma distribution generator
Execution completed"
Best,