2017-07-21 99 views
2

以下のデータに対して離散選択モデリングを実行しようとしています。基本的に、30人のお客様には16種類のピザが用意されています。彼らは1種類以上のピザを選ぶことができ、選択したピザは選択肢の変数で示されます。R mlogitはsolve.defaultでエラーをスローします(H、g [!fixed]):システムは計算上単数です:逆数条件番号

pizza cust choice pan thin pineapple veggie sausage romano mozarella oz 
1  1 Cust1  0 1 0   1  0  0  1   0 1 
2  2 Cust1  1 0 1   1  0  0  0   0 0 
3  3 Cust1  0 0 0   1  0  0  0   1 1 
4  4 Cust1  1 0 1   1  0  0  0   0 0 
5  5 Cust1  1 1 0   0  1  0  0   0 1 
6  6 Cust1  0 0 1   0  1  0  1   0 0 
7  7 Cust1  0 0 0   0  1  0  0   0 1 
8  8 Cust1  1 0 1   0  1  0  0   1 0 
9  9 Cust1  0 1 0   0  0  1  0   1 0 
10 10 Cust1  1 0 1   0  0  1  0   0 1 
11 11 Cust1  0 0 0   0  0  1  1   0 0 
12 12 Cust1  0 0 1   0  0  1  0   0 1 
13 13 Cust1  0 1 0   0  0  0  0   0 0 
14 14 Cust1  1 0 1   0  0  0  0   1 1 
15 15 Cust1  0 0 0   0  0  0  0   0 0 
16 16 Cust1  0 0 1   0  0  0  1   0 1 
17  1 Cust10  0 1 0   1  0  0  1   0 1 
18  2 Cust10  0 0 1   1  0  0  0   0 0 
19  3 Cust10  0 0 0   1  0  0  0   1 1 
20  4 Cust10  0 0 1   1  0  0  0   0 0 

私は以下のコマンドを使用してデータを変換します。私はchid.var = "chid"alt.levels=c(1:16)を追加するようなここではほとんど変更を加えませんでした。 alt.levelsalt.varの両方を使用すると、ピザが既に存在し、置き換えられるというエラーが表示されます。しかし、どちらかを使用するとエラーは発生しません。

最後に、mlogitコマンドを使用すると、このエラーが発生します。

mlogit(choice ~ pan + thin + pineapple + veggie + sausage + romano + mozarella + oz, pz) 

Error in solve.default(H, g[!fixed]) : 
    system is computationally singular: reciprocal condition number = 8.23306e-19 

これは私の最初の投稿です。私はこのサイトを頻繁に訪れ、すでに解決策が見つかったので投稿する必要はありませんでした。私はthisのようなほぼすべての同様の投稿を1つしか無駄に使いました。私は離散選択モデリングを初めて学んだので、私がここで根本的な間違いをしているかどうかは分かりません。

また、私は実際にchid.varが何をしているのかは分かりません。

答えて

0

この問題を解決できませんでした。あなたはnnetパッケージから多機能を利用することができます。それは動作するようです。答えを確認しました。問題の選択はあなたが予測したい従属カテゴリ変数である任意の変換

library("nnet") 
pizza_model <- multinom(choice ~ Price + IsThin + IsPan ,data=pizza_all) 
summary(pizza_model) 

ためので必要性を全く示さないよう

データセットは同じまま。 Price、IsThin、およびIsPanは独立変数です。出力は以下の通りです

Call: 
multinom(formula = choice ~ Price + I_cPan + I_cThin, data = pizza_all) 

Coefficients: 
        Values Std. Err. 
(Intercept) 0.007192623 1.3298018 
Price  -0.149665357 0.1464976 
I_cPan  0.098438084 0.3138538 
I_cThin  0.624447867 0.2637110 

Residual Deviance: 553.8519 
AIC: 561.8519 
+0

お手数ですが、お手数ですが下記の通りです。将来のあなたの質問の読者のためにあなたのモデルが 'nnet'で動作するようになっていく方法の詳細を投稿できますか? – DirtStats

関連する問題