Re: [問題] logistic regression model

[west1996 (焦了六年變脆了)]

※ 引述《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
推 wlsherica:descending 的問題??
推 west1996:應該是minitab和SAS fit 的model不一樣，一個是fit
→ west1996:P(Y=0|X)，一個是fit P(Y=1|X)，這樣會剛好差一個負號
→ forever0322:可是SPSS跑出來的也跟MINITAB一樣
→ west1996:因為SPSS和minitab fit一樣的model阿，你外部檔book1.xls
→ west1996:的第一筆資料Y是多少？
→ forever0322:我剛剛改一下將01顛倒在MINITAB跑的確是差正負號
→ forever0322:第一筆是0
推 west1996:把第一筆跟後面換掉讓第一筆是1的話SAS就會跑出SPSS的那
→ west1996:組參數了
→ forever0322:這樣是fit P(Y=1|X)嗎
推 west1996:恩  SAS中的預設值是去fit讀進來的第一個
→ west1996:要改順序的話可以下order=來改，可以參考help
→ forever0322:可是我剛剛把後面一筆Y=1調到第一筆跑出來還是一樣
→ west1996:@@
用回的好了
你用同樣的資料跑兩種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