2014-12-03 81 views
4

我有离散选择实验(DCE)的数据,查看来自不同行业的个人的招聘偏好。我已经格式化为长格式。我想使用mlogit进行建模。我已经导出的数据,并可以使用asclogit命令成功运行模型在Stata,但我无法得到它在R.mlogit:错误的值需要TRUE/FALSE

这里运行是第25行数据的快照:

> data[1:25,] 
    userid chid item sector outcome cul fit ind led prj rel 
1 11275 211275 2  1  1 0 1 0 1 1 1 
2 11275 211275 2  2  0 1 0 0 0 0 0 
3 11275 211275 2  0  0 0 0 1 1 0 1 
4 11275 311275 3  0  1 1 1 0 0 0 1 
5 11275 311275 3  2  0 0 1 0 0 0 1 
6 11275 311275 3  1  0 0 1 0 0 0 0 
7 11275 411275 4  0  0 1 0 1 1 0 0 
8 11275 411275 4  2  1 0 1 1 1 1 0 
9 11275 411275 4  1  0 0 1 0 1 0 0 
10 11275 511275 5  1  1 1 0 1 0 1 1 
11 11275 511275 5  2  0 0 0 1 1 0 0 
12 11275 511275 5  0  0 0 0 1 1 1 0 
13 11275 611275 6  0  0 0 1 1 0 0 1 
14 11275 611275 6  1  1 1 1 1 0 0 1 
15 11275 611275 6  2  0 1 1 1 0 1 0 
16 11275 711275 7  1  0 0 0 0 0 1 0 
17 11275 711275 7  0  0 1 0 0 1 1 0 
18 11275 711275 7  2  1 1 0 0 1 1 1 
19 11275 811275 8  0  1 0 1 0 0 1 1 
20 11275 811275 8  1  0 1 0 1 1 1 1 
21 11275 811275 8  2  0 0 0 0 0 1 1 
22 11275 911275 9  0  0 1 1 0 0 1 0 
23 11275 911275 9  2  1 1 1 1 1 0 1 
24 11275 911275 9  1  0 1 0 1 1 0 0 
25 11275 1011275 10  0  0 0 0 0 0 0 0 

userid和chid是因子变量,其余都是数字。变量: 用户ID是唯一的应答者ID CHID每申请人独特的选择组ID 项被选择集ID(它们被整个受访重复) 扇区是替代品(3个不同的扇区) 结果是替代由申请人在给定的选择的选择集 cul-rel是二元因子变量,根据实验设计,备选方案在不同备选方案中有所不同。

这里是我的mlogit语法:

mlogit(outcome~cul+fit+ind+led+prj+rel,shape="long", 
     data=data,id.var=userid,chid.var="chid", 
     choice=outcome,alt.var="sector") 

这里是我的错误:

Error in if (abs(x - oldx) < ftol) { : 
    missing value where TRUE/FALSE needed 

我已确认没有丢失数据,而且每个选择集都只有1选择替代品。 当模型在Stata中使用完全相同的数据集运行良好时,为什么会出现此错误的任何想法?我可能在某处误读了mlogit语法。如果有帮助,我的Stata语法如下: asclogit结果适用于指定适用的led prj,case(chid)选项(扇区)

回答

-2

您可能需要使用mlogit.data()来塑造数据。有一个例子在?mlogit。希望有所帮助。

+0

我试过了,结果相同。我还在上面的R代码中提供了这些参数。 – 2014-12-03 14:08:45

5

回答我自己的问题,就像我想出来的一样。

当选择集中没有任何选项被选中时,m mlogit无法处理。 R还需要正确排序的数据,选择集中的每个选项都必须连续排列。由于一些数据管理,我没有这样做。有趣的是,Stata可以处理这两种情况,所以这就是我的Stata命令起作用的原因。

另外,对于那些有兴趣的人,Stata的asclogit和R的mlogit给出了完全相同的结果。当这种情况发生时总是很好

+0

对不起,再次挖掘它,但我不知道该做什么时,有不同的选择集。即使替代选项在设置中不可用,R是否需要每个备选选项设置一行? Stata处理这个问题吗? – gannawag 2017-02-21 20:21:46