作者sinasass (廣廣)
看板Statistics
標題[程式] SAS的nested-factorial experiments
時間Fri May 25 23:13:29 2012
[軟體程式類別]:
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