I have a data set that is formatted according to mlogit's standards using mlogit.data command in Rstudio. Here is the first three rows of the dataset.
Trip SevereEarthquake Night Age Mode 1.NTG 1 0 0 18 FALSE 1.TGNV 1 0 0 18 FALSE 1.TGV 1 0 0 18 TRUE
After some variable creation which I do not include here to make the question short, I ran a multinomial logit model using the code below and it worked fine.
mlogit(Mode ~ SE.TGV + SE.TGNV + AGE.NTG, data, reflevel = "NTG")
However when I define nests in the code (as below) to do a nested logit model I get an error:
mlogit(Mode ~ SE.TGV + SE.TGNV + AGE.NTG, data, reflevel = "NTG", nests = list(notrip = "NTG", trip = c("TGV","TGNV"))) Error in solve.default(crossprod(attr(x, "gradi")[, !fixed])) : system is computationally singular: reciprocal condition number = 8.87901e-37
What is causing this problem? How can I solve this?