作者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