求解答，关于category和dummy

blliufan

4193

收藏 2012-05-21

在做回归分析时，关于category variable，是否要转换成dummy variable。
用class和用转换过的dummy variable求出的系数不一样啊。

我该何去何从呢

扫码加我拉你入群

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

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

全部回复

blliufan

2012-5-21 13:25:51

比如marriage是category，single，married，divorce。
方法一
proc genmod data=info;
class marriage;
model y=marriage;
run;
方法二
proc genmod data=info;
model y=married single;
run;
这里married和single是转换过的dummy variable。
两者结果明显不一，这是为何

扫码加我拉你入群

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

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

hcydlee

2012-5-21 14:18:02

两者的reference group应该是不同的吧，你手动产生的dummy variables 和系统自动产生的 reference group应该不一样。

扫码加我拉你入群

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

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

blliufan

2012-5-21 14:34:34

hcydlee 发表于 2012-5-21 14:18
两者的reference group应该是不同的吧，你手动产生的dummy variables 和系统自动产生的 reference group应 ...

那用哪个好呢？还是说其实没有什么大的差别

扫码加我拉你入群

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

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

hcydlee

2012-5-21 14:37:11

比较喜欢用手动的，自己设置的，解释起来也轻松。

扫码加我拉你入群

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

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

blliufan

2012-5-21 15:23:33

hcydlee 发表于 2012-5-21 14:37
比较喜欢用手动的，自己设置的，解释起来也轻松。

谢谢啊。还有问题请教啊，我用sample1得到自变量的parameter，然后要预测sample2，这个该怎么做。

扫码加我拉你入群

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

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

点击查看更多内容…

hcydlee

2012-5-21 16:36:38

blliufan 发表于 2012-5-21 15:23
谢谢啊。还有问题请教啊，我用sample1得到自变量的parameter，然后要预测sample2，这个该怎么做。

http://support.sas.com/kb/33/307.html 这篇文章的Example 1和你的几乎一样。你可以参考一下，基本方法就是用training set 拟合好模型，然后把validation set 的因变量变为missing，和原training set 合并在一起,用合并后的set再拟合一次（validation set的因变量是missing，所以不会用到，模型结果和traning set 的一样），通过设置把来自于validation set的observation做个预测。具体内容你可以仔细的看下这篇文章。

记得帮我评分哦

扫码加我拉你入群

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

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

hcydlee

2012-5-21 17:50:48

hcydlee 发表于 2012-5-21 14:18
两者的reference group应该是不同的吧，你手动产生的dummy variables 和系统自动产生的 reference group应 ...

看了下user guide.是不一样的。GENMOD, GLMSELECT, LOGISTIC, and PHREG.产生哑变量的方法是由class 的param= option控制的，默认是PARAM=GLM，我们一般手动产生哑变量的方法对应 param=REF.(默认REF=last)

扫码加我拉你入群

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

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

pl172231

2012-5-21 19:30:36

hcydlee 发表于 2012-5-21 17:50
看了下user guide.是不一样的。GENMOD, GLMSELECT, LOGISTIC, and PHREG.产生哑变量的方法是由class 的 ...

厉害

扫码加我拉你入群

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

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

blliufan

2012-5-21 19:54:46

hcydlee 发表于 2012-5-21 16:36
http://support.sas.com/kb/33/307.html 这篇文章的Example 1和你的几乎一样。你可以参考一下，基本方法 ...

如果在的话，请教下，做genmod会在result窗口输出结果，但我要做cross validation，不想每一次都要在results窗口处结果，只要能输出data set就行。该怎么办？

扫码加我拉你入群

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

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

hcydlee

2012-5-21 21:05:15

blliufan 发表于 2012-5-21 19:54
如果在的话，请教下，做genmod会在result窗口输出结果，但我要做cross validation，不想每一次都要在resu ...

一般来说你可以在程序上加
ods trace on；

然后程序结束加上
ods trace off；

然后在log里面看到
Output Added:
-------------
Name: *****
*****
你只要再跑一次程序，程序结尾写上：
ods output Name后面的名字=你想存储的data set名字；
就可以把结果输出到data里面了。

建议你看下the little sas book.

扫码加我拉你入群

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

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

blliufan

2012-5-21 21:23:01

hcydlee 发表于 2012-5-21 21:05
一般来说你可以在程序上加
ods trace on；

好像不行啊，我做手动的backwards selection。每次结果就输出在results窗口了，我这要做1000validation，会不会出问题啊

扫码加我拉你入群

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

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

hcydlee

2012-5-21 21:38:36

blliufan 发表于 2012-5-21 21:23
好像不行啊，我做手动的backwards selection。每次结果就输出在results窗口了，我这要做1000validation， ...

我这个方法是在你不熟悉sas输出结果的时候，用trace追踪输出的目标在哪里，然后把这个目标输出到数据集里面去。第一步是了解SAS把结果输出到哪里去了，第二部是让SAS把你要的数据输出到你指定的地方去。
你的问题的话：
ods output ParameterEstimates=dataset（你想要输出的数据集的名字）; 这个是参数估计结果。还有其他很多结果你可以看看help documents.或者看看http://support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#statug_genmod_sect058.htm这个表里面是所有genmod可以产生的结果的汇总。

扫码加我拉你入群

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

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

blliufan

2012-5-22 17:48:41

hcydlee 发表于 2012-5-21 21:38
我这个方法是在你不熟悉sas输出结果的时候，用trace追踪输出的目标在哪里，然后把这个目标输出到数据集里 ...

你好啊，谢谢你昨天的帮助。我自己没能描述好问题。其实，我的疑问是这样的。我手动做backwards selection。每次去掉一个variable，然后跑回归，都会在results viewer中显示这次回归的结果，包括p-value的表啊什么的。
看图。

能不让SAS输出么

附件列表

QQ截图20120522174730.jpg

原图尺寸 41.68 KB

results

扫码加我拉你入群

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

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

hcydlee

2012-5-22 18:37:47

blliufan 发表于 2012-5-22 17:48
你好啊，谢谢你昨天的帮助。我自己没能描述好问题。其实，我的疑问是这样的。我手动做backwards selectio ...

这个还是会输出的，实际上你这种情况最好还是写个宏，然后批量把结果输出到数据集里面去。然后你直接看结果数据集，而不是去系统的结果。

扫码加我拉你入群

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

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

blliufan

2012-5-22 19:25:32

hcydlee 发表于 2012-5-22 18:37
这个还是会输出的，实际上你这种情况最好还是写个宏，然后批量把结果输出到数据集里面去。然后你直接看结 ...

是的，我是不看那里的结果，我担忧的是大量的results会不会造成死机

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群