作者west1996 (焦了六年变脆了)
看板Statistics
标题Re: [问题] logistic regression model
时间Thu May 28 22:22:05 2009
※ 引述《forever0322 (水蓝色)》之铭言：
想请问SAS logistic regression model
跑出来的beta0和beaa1值和我用minitab和SPSS差正负号
检查蛮久也是找不出来
PROC IMPORT OUT= book1 /*this is the name of the data set*/
DATAFILE= "I:\123\book1.xls" /*this is where you saved your excel data file*/
DBMS=EXCEL REPLACE;
SHEET="sheet1"; /*sheet1 is the excel files sheet name*/
GETNAMES=yes; /*yes means the first row of the data is the variables names*/
RUN;
proc LOGISTIC ;
MODEL Y = X ;
run;
这是SAS的程式码
我这样KEY有错嘛?
麻烦是否可以给我指正
拜托大家了
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.135.139.143
1F：推 wlsherica:descending 的问题?? 05/28 21:47
2F：推 west1996:应该是minitab和SAS fit 的model不一样，一个是fit 05/28 21:48
3F：→ west1996:P(Y=0|X)，一个是fit P(Y=1|X)，这样会刚好差一个负号 05/28 21:48
4F：→ forever0322:可是SPSS跑出来的也跟MINITAB一样 05/28 21:52
5F：→ west1996:因为SPSS和minitab fit一样的model阿，你外部档book1.xls 05/28 21:55
6F：→ west1996:的第一笔资料Y是多少？ 05/28 21:56
7F：→ forever0322:我刚刚改一下将01颠倒在MINITAB跑的确是差正负号 05/28 21:56
8F：→ forever0322:第一笔是0 05/28 21:57
9F：推 west1996:把第一笔跟後面换掉让第一笔是1的话SAS就会跑出SPSS的那 05/28 21:59
10F：→ west1996:组参数了 05/28 21:59
11F：→ forever0322:这样是fit P(Y=1|X)吗 05/28 22:06
12F：推 west1996:恩 SAS中的预设值是去fit读进来的第一个 05/28 22:07
13F：→ west1996:要改顺序的话可以下order=来改，可以参考help 05/28 22:08
14F：→ forever0322:可是我刚刚把後面一笔Y=1调到第一笔跑出来还是一样 05/28 22:10
15F：→ west1996:@@ 05/28 22:16
用回的好了
你用同样的资料跑两种code
1. proc logistic;
order=internal;
model Y=X;
run;
2. proc logistic;
order=internal desc;
model Y=X;
run;
这两个应该会跑出刚好相反的beta值
第一个是fit P(Y=0|X)
第二个是fit P(Y=1|X)
看看那个是你要的吧
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.230.72.191