[程式] SAS 一處理在一區集中不只一重複

看板Statistics作者Tchieh (默契)時間14年前 (2010/01/20 00:00)推噓2(2推 0噓 3→)

留言5則, 3人參與討論串1/1

這是關於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)

推

andrew43

01/20 01:58, , 1^F

01/20 01:58, 1^F

→

andrew43

01/20 01:58, , 2^F

01/20 01:58, 2^F

→

Tchieh

01/20 10:57, , 3^F

01/20 10:57, 3^F