2014-12-03 26 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