「統計学関連なんでもあり」の過去ログ---028

★ influential dataの検索 ★

2448.　influential dataの検索　鈴木　2004/02/23 (月) 20:29
└2450.　Re: influential dataの検索　青木繁伸　2004/02/23 (月) 22:33
　└2473.　Re^2: influential dataの検索　鈴木　2004/02/24 (火) 20:43

2448.　influential dataの検索　鈴木　　2004/02/23 (月) 20:29

いつもお世話になります。

influential dataの検索というのを丹後俊郎先生他の本からみつけSASでやってみたのですが・・SASは無反応でした。
よくわからないのですが，プログラムが間違っているとSASが教えてくれるのですが・・・
反応無しです。何にもしてくれない！！

以下のプログラムです。
離れ値やおかしなデータを簡単に見つける方法を教えていただけませんか？
%macro inf(data,vard,vlist,ilist); proc catmod data&data; model &vard=&vlist &ilit/ nodesign nogls noiter noprofile pred=freq; response/ out=sasout; keep _sample_ _type_ _obs_ _pred_ &vard &vlist; if &vard=. then delete; if &vard=2 then _pred_=.; proc sort; by _sample_ &vlist; data Tl1; set Tl1; if _pred_=. then delete; M=col1+col2; P=_pred_/M; Y=Log(P/(1.-p)) + (col1-M*p)/(M*P*(1.-p)) ; W=M*P*(1.-p); Drop _sample_ col2 _name_ _pred_; proc glm; class &vlist &ilist ? solition ; output out=gl student=ch cook=cook H=hat; weight w; data gl ; set gl; ch2=ch**2; Proc print; proc plot; plot (ch2=ch**2; run; %mend;

　　　　　[このページのトップへ]

2450.　Re: influential dataの検索　青木繁伸　　2004/02/23 (月) 22:33

私は SAS は知らないのですが，プログラムを見たところ？？のところを

> model &vard=&vlist &ilit/ nodesign nogls noiter noprofile pred=freq;

&ilist では？

> proc glm; class &vlist &ilist ? solition ;

solution では？

> proc plot; plot (ch2=ch**2;

2 と ; の間に ) はなくていいの？

昔からそうですが，ユーザが打ち込んだプログラムが動かないというのは，そのほとんどが，入力ミスです。一文字でも違っていたらアウトです。

SAS がエラーメッセージを出してくれるのは，文法的におかしいときだけです。
文法的に正しくても，論理的におかしければ，期待される結果は得られません。
そのプログラムの論理的正しさから言えば，なんにも出力しないと言うのがそのプログラム正しい動き方なのでしょう。

　　　　　[このページのトップへ]

2473.　Re^2: influential dataの検索　鈴木　　2004/02/24 (火) 20:43

青木先生
お返事ありがとうございました。

ご助言どおりSAS再チャレンジしてみます。

　　　　　[このページのトップへ]