7

我有一个长格式的联合调查数据。前几行如下所示：

 ID alt choice size tar length brand flavor gender age yr_smoke num_smoke job
 1  1   1     no    1   1      1     1      1      2  35       10        20   1
 2  1   2     no    1   1      1     1      1      2  35       10        20   1
 3  1   3     no    1   1      1     1      1      2  35       10        20   1
 4  1   4     no    1   1      1     1      1      2  35       10        20   1
 5  1   5     no    1   1      1     1      1      2  35       10        20   1
 6  1   6     no    1   1      1     1      1      2  35       10        20   1

我使用 mlogit.data 作为：

data_mlogit_ct1_test2 <- mlogit.data(data_mlogit_ct1_test1,choice="choice",
shape="long",alt.var="alt")

转换后的数据的前几行如下所示：

   ID alt choice size tar length brand flavor gender age yr_smoke num_smoke job
   1.1  1   1  FALSE    1   1      1     1      1      2  35       10        20   1
   1.2  1   2  FALSE    1   1      1     1      1      2  35       10        20   1
   1.3  1   3  FALSE    1   1      1     1      1      2  35       10        20   1
   1.4  1   4  FALSE    1   1      1     1      1      2  35       10        20   1
   1.5  1   5  FALSE    1   1      1     1      1      2  35       10        20   1
   1.6  1   6  FALSE    1   1      1     1      1      2  35       10        20   1

现在我正在尝试将 mlogit 拟合为：

 fit_mlogit_ct11 <- mlogit(choice~size+tar+length+brand+flavor, 
                           data_mlogit_ct1_test1,shape="long",chid.var="ID",
                           alt.var="alt",method="bfgs",heterosc=TRUE,tol=10)

我收到了这个错误：

solve.default(crossprod(attr(x, "gradi")[, !fixed])) 中的错误：Lapack 例程 dgesv：系统完全是奇异的

然后我做了：

  fit_mlogit_ct11 <- mlogit(choice~size+tar+length+brand+flavor, data_mlogit_ct1_test2)

现在我得到了：

solve.default(H, g[!fixed]) 中的错误：系统在计算上是奇异的：倒数条件数 = 3.4767e-18

有人能帮帮我吗？

4

1 回答 1

11

我认为问题在于您没有正确定义您的个人特定变量。根据包手册，当您使用 mlgit 功能构建模型时，您可以使用：| 将替代特定变量与单个特定变量分开。例如：

fit_mlogit_ct11 <- mlogit(choice~1|size+tar+length+brand+flavor, 
                           data_mlogit_ct1_test1,shape="long",chid.var="ID",
                           alt.var="alt",method="bfgs",heterosc=TRUE,tol=10)

这只是说您的所有变量都是特定于个人的，并且您只接受替代特定的截距。但实际上正如上面已经提到的那样，您的数据看起来很奇怪，因为它们看起来相同，因为您只有单个特定变量。这最终会导致奇点（您的变量与 1 相互关联）。

于 2013-01-15T11:05:46.293 回答