我有离散选择实验(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)选项(扇区)
我试过了,结果相同。我还在上面的R代码中提供了这些参数。 – 2014-12-03 14:08:45