STATA mixlogit collinear 问题

2014-2-23 16:45:52

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-24 03:52:08

if you do

_rmcoll color size firmness sweetness flavor shelf cost

does it show anything?

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-24 11:18:21

jjjj6666 发表于 2014-2-24 03:52
if you do

_rmcoll color size firmness sweetness flavor shelf cost

Hi jjjj6666, 执行了你建议的口令后，返回的是“note: cost omitted because of collinearity”。可是我不想omit cost啊，如果不去掉任何变量，有没有解决的办法？谢谢！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-24 15:01:23

继续等待高手的帮助！我用mixlogit命令分析了其他几个类似的数据，解释变量中cost是两个数值，其他是0和1虚拟变量，解释变量之间的相关系数也类似，STATA都运行了，没有出现任何报错。先谢了！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-25 01:11:16

if force collinearity term, is the model unidentifiable? At least, if force it, the model estimating may not converge.

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

2014-2-25 01:54:37

jjjj6666 发表于 2014-2-25 01:11
if force collinearity term, is the model unidentifiable? At least, if force it, the model estimatin ...

Thanks for the reply! I want to have a try "force the collinearity term", but how do I do it in STATA? Could you let me know the code to fulfill it? Thank you so much!

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-25 02:04:15

the ado file doesn't allow it. You need to modify the ado file (I think stata ml optimization can have the option to allow it), which may take some works. But if looking the likelihood function, it seems the model itself has problem under this case, is this correct?

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-25 03:48:56

jjjj6666 发表于 2014-2-25 02:04
the ado file doesn't allow it. You need to modify the ado file (I think stata ml optimization can ha ...

Sorry but I do not see problems of the model. Could you give me a little more detail? Thanks!

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-25 03:49:32

真是抱歉，我上传的数据中没有count这个变量，请重新下载数据。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-25 14:49:24

你好，我也是用到mixed logit 模型，我想请教一下，用stata可以实现mixlogit，吧，这个里面group()这项是什么意思，我用我自己的数据运行，出现错误，没有结果。请你一定帮帮忙啊

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-28 08:13:54

xiadan1025 发表于 2014-2-25 14:49
你好，我也是用到mixed logit 模型，我想请教一下，用stata可以实现mixlogit，吧，这个里面group()这项是什 ...

我也是第一次用mixed logit，第一次用stata，所以我也不是很明白，现在用的程序是别人编好的。我刚刚查了一下stata对group的解释：“ group(varname) is required and specifies a numeric identifier variable for the choice occasions”. 我的理解是group代表每个人所面临的选择题，所有人一起排序，我们感兴趣的是每个人选的是每个选择题里哪个选项和为什么，才用mixed logit吧，这个group就是选择题的编号。例如我的问题中每个题有两个选项，所以我先生成一个变量egen count=fill (1 1 2 2)，这样count这个变量（赋值是11223344...）就是所有人一起排的选择题的编号（每个题有两个选项，所以每个题号重复两次），那么mixlogit 口令中就是 group（count）。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-2-28 08:23:01

jjjj6666 发表于 2014-2-25 09:07
what I think is if you add all the terms, the model estimating will not converge based on how the ML ...

谢谢你一直的回复！You have partially solved my problem! But since no one else except you came to help, I will give you the credit

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-3-11 10:03:01

ivylee_777 发表于 2014-2-28 08:13
我也是第一次用mixed logit，第一次用stata，所以我也不是很明白，现在用的程序是别人编好的。我刚刚查了 ...

太谢谢你的回答了，我的论文是三分类的，太难懂了，
谢谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-3-12 10:35:31

ivylee_777 发表于 2014-2-28 08:13
我也是第一次用mixed logit，第一次用stata，所以我也不是很明白，现在用的程序是别人编好的。我刚刚查了 ...

还有，你研究的问题是2分类的吗？能不能把你的程序借我参考一下，我知道这个请求有点过，但是实在太着急了，请亲一定要帮帮忙啊，感激感激

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-3-14 03:09:02

xiadan1025 发表于 2014-3-12 10:35
还有，你研究的问题是2分类的吗？能不能把你的程序借我参考一下，我知道这个请求有点过，但是实在太着 ...

我的数据处理是用SAS做的，code大致是这样的：
data YOUR DATA TO BE USED IN STATA; /*输入你自己的数据文件名字*/
set YOUR ORIGINAL DATA; /*你原始数据文件名字*/
RETAIN ID 0;
ID+1;

decision = (s1fp=1); /*S1FP 是你的第一个选择题，1 是第一个选项，下面的2是第二个选项*/
choice1=1;choice2=0;choice3=0;option=1;Qnum=1;
Appearance=0; /*这些是选择1时，虚拟变量的值，这里的APPEARANCE等是我的自变量（都是虚拟变量，你要换成你自己的*/
Size=1;
Firmness=1;
Sweetness=0;
Flavor=1;
Cost=2.19;
output;
decision = (s1fp=2);
choice1=0;choice2=1;choice3=0;option=2;Qnum=1;
Appearance=1;
Size=0;
Firmness=0;
Sweetness=1;
Flavor=0;
Cost=1.69;
output;

decision = (s2fp=1); /*S2FP 是你的第一个选择题，1 是第一个选项，下面的2是第二个选项*/
choice1=1;choice2=0;choice3=0;option=1;Qnum=2;
Appearance=0;
Size=1;
Firmness=0;
Sweetness=1;
Flavor=1;
Cost=2.19;
output;
decision = (s2fp=2);
choice1=0;choice2=1;choice3=0;option=2;Qnum=2;
Appearance=1;
Size=0;
Firmness=1;
Sweetness=0;
Flavor=0;
Cost=1.69;
output;

. /* 有几个选择题，上面的code就重复几遍*/
.
.
run;

然后就是在STATA里：
insheet using "路径\YOUR DATA TO BE USED IN STATA.csv", clear

** generate count
egen count=fill (1 1 2 2)

** mixed logit
mixlogit decision, group(count) id(id) rand(appearance size firmness sweetness flavor shelflife cost)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-3-14 09:18:41

ivylee_777 发表于 2014-3-14 03:09
我的数据处理是用SAS做的，code大致是这样的：
data YOUR DATA TO BE USED IN STATA; /*输入你自己的数据 ...

谢谢你，我把你的数据下载下来，看过了
在我理解的一个人在两个选择里面选择一个，那他其他变量的值应该不变，只有decision s是变得
我知道我理解的你是数据应该是下面
ID decision choice1 choice2 choice3 option Qnum Color Size
1 0 1 0 0 1 6 1 0
1 1 1 0 0 1 6 1 0
我实在不懂了，如果一百个人，做选择，每个选择两个选项，那应该数据是有200行的，那数据除了desion这项，那么其他的变量一定是要重复的呀，求解答，太困扰我了。感激不尽，感激不尽~~~

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2014-3-14 16:59:47

ivylee_777 发表于 2014-3-14 03:09
我的数据处理是用SAS做的，code大致是这样的：
data YOUR DATA TO BE USED IN STATA; /*输入你自己的数据 ...

还想问你数据里面的ID是怎么赋值的难道是配对吗？求解答，
请一定帮帮忙
very very thanks

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝