0

離散選択実験のデータを分析していますmlogitが、指定したときに重みが何を使用するのかわかりませんweights

次のコード:

mlogit(formula = RES ~ -1 + V1 + V2, data = data, 
    reflevel = 1, rpar = c(V1 = "n", V2 = "n"), weights = Weight1, correlation = FALSE, 
    halton = NA, panel = TRUE, seed = 1234567890, method = "bfgs")

次の推定値を生成します。

Frequencies of alternatives:
      1       2 
0.22987 0.77013 

bfgs method
19 iterations, 0h:15m:34s 
g'(-H)^-1g = 4.29E-08 
gradient close to zero 

Coefficients :
      Estimate Std. Error t-value  Pr(>|t|)    
V1     0.859789   0.019076  45.072 < 2.2e-16 ***
V2     2.705395   0.039205  69.006 < 2.2e-16 ***
sd.V1  0.483573   0.023502  20.576 < 2.2e-16 ***
sd.V2  3.916796   0.062557  62.612 < 2.2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Log-Likelihood: -9297.9

random coefficients
       Min.    1st Qu.    Median      Mean  3rd Qu. Max.
V1  -Inf 0.53362451 0.8597892 0.8597892 1.185954  Inf
V2  -Inf 0.06355681 2.7053955 2.7053955 5.347234  Inf

ただし、Stata で同じ混合ロジット モデルを実行すると、次のコマンドが実行されます。

mixlogit res [pweight=weight1], group(str) id(id) rand(V1 V2) ln(0)

次の見積もりが得られます。

Mixed logit model                               Number of obs     =     41,154
                                                Wald chi2(2)      =     395.55
Log likelihood = -9089.7906                     Prob > chi2       =     0.0000

------------------------------------------------------------------------------
             |               Robust
         res |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Mean         |
          V1 |   1.207748   .0774815    15.59   0.000     1.055887    1.359608
          V2 |   4.458814   .2356245    18.92   0.000     3.996998    4.920629
-------------+----------------------------------------------------------------
SD           |
          V1 |   1.107036   .0765884    14.45   0.000     .9569252    1.257146
          V2 |   4.444472   .3586858    12.39   0.000     3.741461    5.147483
------------------------------------------------------------------------------

Stata で使用する重み付けスキーム ( pweight、iweight、またはfweight) に関係なく、同様の結果が得られますが、R が提供する結果は得られません。

ただし、いずれかのプログラムで重み付けされていない混合ロジット モデルを実行すると、同じ推定値が得られます。これにより、重み付けが明らかな問題であることがわかりますが、Rが何をしているのかわかりません。

ヘルプ？

4

1 に答える 1

1

これは私にしか見えないかもしれませんが、誰かがこれに遭遇した場合に備えて、問題は重み付けではありませんでした. Rコマンド( R = 10, tol = 10)の停止条件が甘すぎたのが原因です。R大きい値と小さい値を設定するとtol、2 つの推定値は収束します。

于 2016-02-15T17:43:30.233 に答える