3

我正在尝试使用该mlogit包来估计多项概率模型，但出现以下错误。

Error in if (is.null(initial.value) || lnl <= initial.value) break : 
      missing value where TRUE/FALSE needed

这些是我使用的语句

dat <- read.table("C:\\Users\\Data.txt", header=TRUE, sep="\t", row.names=NULL)
brand <- mlogit.data(dat, shape="long", choice = "choice", varying = 3:7, sep = " ", alt.levels=1:7, id.var = "id")
brand.mprobit <- mlogit(choice ~ last + price + ad + shelf | pers_1 + pers_2 + spend, probit = TRUE)
summary(brand.mprobit)

我的数据如下所示：

  obs   choice  person  option  last    price   ad  shelf   pers_1  pers_2  spend   ago
1   0   101 1   0   0.38    1   1   1   9   20  4
1   0   101 2   0   0.44    0   0   1   9   20  4
1   0   101 3   0   0.33    0   0   1   9   20  4
1   0   101 4   0   0.64    0   0   1   9   20  4
1   1   101 5   0   0.55    0   0   1   9   20  4
1   0   101 6   1   0.18    0   0   1   9   20  4
2   1   102 1   0   0.99    0   1   1   9   23  5
2   0   102 2   0   0.66    0   0   1   9   23  5
2   0   102 3   0   0.45    1   1   1   9   23  5
2   0   102 4   0   0.56    0   0   1   9   23  5
2   0   102 5   0   0.78    0   0   1   9   23  5
2   0   102 6   1   0.18    0   0   1   9   21  5

最后 4 列仅是个人特定的，而其他列是替代特定的。有人可以帮忙吗？

4

0 回答 0