2013-04-12 50 views
0

我目前正在学习多项式对数估计量,我想用mlogit包来估计R中的一个模型。在广泛阅读了该主题后,很明显,该过程中的一个重要步骤是使用mlogit.data()函数重新整理数据。我的数据帧包含以下信息:整形多项式回归的数据

  1. satisfaction:这是从1(非常不满意)分类变量到5(非常满意
  2. education:这是多年的个人
  3. 教育的数量
  4. country:这是个别
  5. average_income原产国:这是该国的平均收入

这里是数据帧的视觉表示形式:

> dat 
    country satisfaction education income 
1  1   3  12 750 
2  1   5  13 750 
3  1   2  10 750 
4  3   4  13 675 
5  3   5  14 675 
6  3   4  11 675 
7  3   1  14 675 
8  2   1  11 820 
9  2   5  14 820 
10  2   3  12 820 

似乎该mlogit()函数不喜欢它以这种形式。我试图用satisfaction作为因变量并且education作为自变量来估计模型。

我该如何重塑它才能使它工作?

+0

请发布您的代码和错误消息。 –

+1

我不清楚你的回应变量是什么。它是“满意”吗?如果是这样,似乎[*序数*逻辑回归](http://en.wikipedia.org/wiki/Ordered_logit)比多项式更合适。加州大学洛杉矶分校统计网站有一个[有序逻辑回归在R](http://www.ats.ucla.edu/stat/r/dae/ologit.htm),也可能对你有所帮助。 – gung

+1

请考虑在“nnet”包中使用“multinom”函数。 'multinom(satisfaction〜...,data = dat)'会自动将'满意度'重构为一个多项式数组,其中* minimum *被视为参考级别(与SAS使用数据中的第一个值相反)。 – AdamO

回答

0

我将包括ALT变量将表明替代1至5因变量是“模式”变量应该是TRUEFALSE你重塑你的数据到mlogit格式后。