2016-11-17 31 views
0

問題是,我正在比較兩種劑量水平的藥物與安慰劑的有效性。安慰劑是對照組,我有一個低劑量組和高劑量組。我正在嘗試爲每個單獨的組獲得單獨的直方圖和正常的分位數圖。我沒有收到任何錯誤,但我也沒有從SAS本身找到任何有用的信息。有人知道嗎?如何獲得SAS中每個劑量組的單獨圖表?

這裏是我的代碼:

Data Memory; 
input dosegrp $ funct @@; datalines; 
Pla 6 Pla 5 Pla 6 Pla 8 Pla 3 Pla 4 Pla 7 Pla 4 Pla 7 Pla 6 Pla 7 Pla 8 Pla 5 Pla 6 Pla 5 Pla 5 Pla 7 Pla 8 Pla 5 Pla 9 Pla 11 Pla 4 Pla 7 
lo 8 lo 12 lo 7 lo 8 lo 5 lo 6 lo 6 lo 5 lo 3 lo 6 lo 9 lo 6 lo 11 lo 8 lo 6 lo 9 lo 11 lo 5 
hi 11 hi 7 hi 7 hi 11 hi 9 hi 9 hi 12 hi 13 hi 9 hi 13 hi 10 hi 12 hi 9 hi 15 hi 12 hi 14 hi 15 hi 12 
; 
proc sort data = Memory; by dosegrp; 
proc means mean std n data = Memory; 
by dosegrp; 
var funct; 
title1 'one-way ANOVA'; 
title2 'Homework 8 Memory function'; 
RUN; 
title 'Tests for normality'; 
proc univariate normal data = Memory; 
by dosegrp; 
var funct; 
run; 
ods graphics on; 
proc glm data = Memory plots=all; 
class dosegrp; 
model funct=dosegrp; means dosegrp/hovtest welch t dunnett('Pla'); 
contrast 'active vs placebo' dosegrp 0.5 0.5 -1; 
run; 

回答

1

我覺得你的問題是如何獲得直方圖和分位數圖的一個變量單獨組。

如果是這種情況,那麼您確實擁有正確的proc,並且可以使用PROC UNIVARIATEhistogramqqplot語句。

title 'Tests for normality'; 

proc univariate normal data=Memory; 
    by dosegrp; 
    histogram funct; 
    qqplot funct; 
run; 

有例子中的SAS文檔中更多的使用實例(例如見4.28): http://support.sas.com/documentation/cdl/en/procstat/66703/HTML/default/viewer.htm#procstat_univariate_syntax30.htm

+0

Reeza,這似乎已經奏效。只要真的很快,如果我想幹樹葉圖,我需要指定'ods graphics off',對吧? – Jabernet

+1

是的,這些是使用打印機繪圖生成的,類似於點陣打印輸出。 – Reeza

相關問題