這是關於SAS的指令撰寫問題 T^T 我快瘋掉了..... ------------------------------------------------------------------------ [軟體程式類別]: SAS [程式問題]: 區集設計分析 指令的撰寫 [軟體熟悉度]: 新手（不到一個月） [問題敘述]: 不好意思，我在處理作業時遇到了重大打擊 orz 不知道怎樣輸入指令才能跑出有意義的東西 我的資料如下，假設是關於四種飼料養雞的試驗 由於公雞跟母雞可視為異質材料，因此試驗材料主要可分為兩大區集（公、母雞） 而每種飼料在單一區集中有三個重複 A B C D A B C D 雌 Ⅰ 5 2 3 6 雄 Ⅰ 6 4 4 8 雞 Ⅱ 2 4 7 8 雞 Ⅱ 6 8 9 10 Ⅲ 6 3 4 5 Ⅲ 7 2 5 6 可是我們上課只有教RCBD(隨機完全區集設計)的指令 我不會寫區集內重複數不為一的指令阿阿 orz 因為重複數不為一好像就不是RCBD了.... 我到底該怎樣寫>"< 要把它視為什麼資料阿.... 若把資料稍微整合一下，大概是這樣 區　集 雌雞 雄雞 小樣本 Ⅰ Ⅱ Ⅲ Ⅰ Ⅱ Ⅲ A 5 2 6 6 6 7 品 B 2 4 3 4 8 2 種 C 3 7 4 4 9 5 D 6 8 5 8 10 6 [程式範例]: 我嘗試的解決方法有二 一是直接用課本教的RCBD指令改寫 寫出了 title'RCBD with subsampling final'; data subsamp; input block $ var $ subsamp $ obs @@; cards; female A 1 5 female A 2 2 female A 3 6 male A 1 6 male A 2 6 male A 3 7 female B 1 2 female B 2 4 female B 3 3 male B 1 4 male B 2 8 male B 3 2 female C 1 3 female C 2 7 female C 3 4 male C 1 4 male C 2 9 male C 3 5 female D 1 6 female D 2 8 female D 3 5 male D 1 8 male D 2 10 male D 3 6 proc glm; class block var; model obs=block var ; means var/lsd; means var/duncan ; means var/duncan alpha=0.01; run; 可以跑，而且顯示區集跟處理都有5%顯著差異 但是我不知道這樣可不可以orz 或我該怎樣解釋 有區集內重複數不為一這種試驗嗎？ （我想應該有，但是我不知道他應該稱為什麼>"<） 再來 我嘗試將資料視為二次分段抽樣的RCBD 寫出了 title'RCBD with subsampling final'; data subsamp; input block $ var $ subsamp $ obs @@; cards; female A 1 5 female A 2 2 female A 3 6 male A 1 6 male A 2 6 male A 3 7 female B 1 2 female B 2 4 female B 3 3 male B 1 4 male B 2 8 male B 3 2 female C 1 3 female C 2 7 female C 3 4 male C 1 4 male C 2 9 male C 3 5 female D 1 6 female D 2 8 female D 3 5 male D 1 8 male D 2 10 male D 3 6 proc glm; class block var; model obs=block var block*var; test h=block e=block*var; test h=var e=block*var; means var/duncan e=block*var; means var/duncan e=block*var alpha=0.01; run; 也可以跑 但是結果跟方法一差不少阿阿阿orz 因為就分析來說，這個指令好像讓試驗莫名損失了自由度 也就是他考慮了二次抽樣誤差，但實際上並沒有這個問題 ----------------------------------------------------------------------------- 我真的快瘋掉了 可不可以幫幫我 orz 感激不盡>"< --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 114.45.182.240 ※ 編輯: Tchieh 來自: 114.45.182.240 (01/20 00:04)