[程式] SAS的nested-factorial experiments

看板Statistics作者 (廣廣)時間13年前 (2012/05/25 23:13), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/1
[軟體程式類別]: SAS [程式問題]: 資料處理 [軟體熟悉度]: 低(1~3個月) [問題敘述]: 在model中,若想要單獨看其中一個factor裡面的level間關係可以用snk 但在nested model裡面,若有三個factor:A.B.C 且B在A下 model大致會寫成y=mu+Ai+Bj(i)+Ck+ACik+BCjk(i)+em(ijk) 現在我遇到的問題是 在這model下 可以對A做snk 但好像不能對B做snk?(不是我想要的) 我在書中還看到一個proc nested; 指令 有辦法解決這個問題嘛? 下面附上原本的程式碼(我是用proc glm;) [程式範例]: data okok; input type$ num$ stock$ y; datalines; H 2 1 3279 H 2 1 3262 H 2 1 3246 H 2 2 3294 H 2 2 2974 H 2 2 3157 H 3 1 3527 H 3 1 3136 H 3 1 3253 H 3 2 3440 H 3 2 3356 H 3 2 3240 P 5 1 1904 P 5 1 2166 P 5 1 2058 P 5 2 2188 P 5 2 2105 P 5 2 2379 P 8 1 2464 P 8 1 2595 P 8 1 2303 P 8 2 2429 P 8 2 2410 P 8 2 2685 ; run; proc glm; class type num stock; model y=type | stock | num(type)/p; test H=type E=num(type); test H=stock type*stock E=stock*num(type); means type/snk; means num(type)/snk; run; (這裡的num是在type下 anova表示num和type有顯著差異 所以我想在這個model下用snk去 探討num以及type在level間的關係,但是num的snk跑出來不是我想要的) ----------------------------------------------------------------------------- 先謝謝各位大大了!! -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.242.15.85
文章代碼(AID): #1Flw6San (Statistics)