【文章內(nèi)容簡介】 3 85 163 4 95 258 二維頻數(shù)表 A=4 Frequency| Percent | C Row Pct | Col Pct | 1 | 2 | 3 | 4 | Total +++++ 1 | 0 | 6 | 37 | 26 | 69 B | | | | | | | | | | | . | | | | +++++ 2 | 0 | 2 | 18 | 6 | 26 | | | | | | | | | | | . | | | | +++++ Total 0 8 55 32 95 頻數(shù) 構(gòu)成比 行構(gòu)成比 列構(gòu)成比 定性資料統(tǒng)計描述 —— 非編程方式（分析員模塊） 菜單欄 Globals Analyze Statistical Analysis 菜單欄 File Open by SAS name ?作一維頻數(shù)表 菜單欄 Statistics Descriptive Freguency counts 把變量選入 Frequencies 窗口， OK ?作二維頻數(shù)表 菜單欄 Statistics Table Analysis 分別把變量選入 Row、 Column、 Strata窗口， OK 用 Tables可以選每個格中出現(xiàn)更多的結(jié)果 三、定量資料頻數(shù)表 ? 先用 MEANS過程找出最大值和最小值 MIN= MAX= RANGE= ? 確定組距 (I=)和第一組段的下限 (L=) ? 然后用數(shù)據(jù)步分組 賦值新的變量： Y=INT((XL)/I)*I+L ? 用 FREQ過程作頻數(shù)表： DATA AA。 SET AA5。 Y=INT(()/)*+。 PROC FREQ。 TABLES Y。 RUN。 定量資料作直方圖 PROC GCHART。 VBAR X/LEVELS=11 SPACE=0 。 RUN。 GCHART過程作圖 VBAR 作縱向條圖 ； HBAR作橫向條圖 常用選項有： LEVELS 分組的數(shù)目 WIDTH 指定豎條的寬度，以字符寬度為單位 SPACE 指定豎條之間的距離，直方圖 =0 MIDPOINTS= TO BY 指定組中值 直方圖 —— 非編程方式（分析員模塊） 菜單欄 Globals Analyze Statistical Analysis 菜單欄 File Open by SAS name 菜單欄 Graphs Histogram 把變量選入 Analysis窗口， OK 用 Display 可以增加更多的選項 第四節(jié) 定量資料假設(shè)檢驗(yàn) 一、樣本均數(shù)與總體均數(shù)比較 以往新生兒出生體重平均為 3000g，某醫(yī)院對產(chǎn)婦開展孕期保健指導(dǎo)后，測得 15名新生兒出生體重如下，問開展孕期保健指導(dǎo)后對新生兒體重是否有影響？ 2800 3400 2900 2750 3000 3100 4600 3900 2350 3500 2600 3700 3800 2950 3200 樣本均數(shù)與總體均數(shù)比較用 MEANS過程 DATA AA6。 INPUT X。 Y=X3000。 CARDS。 2800 3400 2900 2750 3000 3100 4600 3900 2350 3500 2600 3700 3800 2950 3200 。 PROC MEANS MEAN STD STDERR T PRT。 VAR Y。 RUN。 結(jié)果 Mean Std Dev Std Error T Prob|T| P，不拒絕 H0，尚不能認(rèn)為開展孕期保健指導(dǎo)對新生兒體重有影響 。 二、配對設(shè)計 t 檢驗(yàn) 為了比較新舊兩個測聲計的測定效果，隨機(jī)選定 10個場地，在每地同時用新舊兩個測聲計測量數(shù)據(jù)如下，問新舊兩個測聲計的測量結(jié)果是否不同？ 場地 1 2 3 4 5 6 7 8 9 10 舊： 87 65 74 95 65 55 63 88 61 54 新： 86 66 77 95 60 53 62 85 59 55 配對設(shè)計 t 檢驗(yàn)用 MEANS過程 DATA AA7。 INPUT X1 X2 。 Y=X1X2。CARDS。 87 86 65 66 74 77 95 95 65 60 55 53 63 62 88 85 61 59 54 55 。 PROC MEANS MEAN STD STDERR T PRT。 VAR Y。 RUN。 結(jié)果 Mean Std Dev Std Error T Prob|T| P，不拒絕 H0，尚不能認(rèn)為新舊兩個測聲計的測定結(jié)果有差別。 三、兩樣本均數(shù)比較 試比較克山病人與健康人的血磷值 (mmol/L)是否不同？ 克山病人 健康人 兩樣本均數(shù)比較用 TTEST過程 DATA AA8。 INPUT X。 IF _N_11 THEN G=1。 ELSE G=0。 CARDS。 。 PROC TTEST。 CLASS G。 VAR X。 RUN。 在 TTEST過程中，分類變量 G只能是兩分類。 條件語句 IF – THEN / ELSE IF 后面是條件的表達(dá)式，如果滿足條件就執(zhí)行 THEN 后面操作，否則 (不滿足條件 )，就執(zhí)行 ELSE 后面的操作 (ELSE可選 )。 IF X4=0 THEN X4=. 。 IF CJ60 THEN S+1。 IF X42 AND X5=1 THEN Y=1。 ELSE Y=2。 _N_ 是 SAS數(shù)據(jù)步中自定義的變量，它的值表示觀測的序號。 IF _N_11 THEN G=1。 ELSE G=0。 也可用 G=(_N_11)。賦值分組變量 G：前 11個觀測 G=0，超過 11時 G=1。 t 檢驗(yàn)結(jié)果 (1) G N Mean Std Dev Std Error 1 11 2 13 (2) Variances方差 T DF Prob|T| Unequal 不齊 Equal 齊 (3) For H0: Variances are equal, 方差齊性檢驗(yàn) F39。 = DF = (12,10) ProbF39。 = 結(jié)論： P，拒絕 H0，接受 H1，可以認(rèn)為兩組人群血磷值不同，克山病人血磷值較高。 兩大樣本均數(shù)比較時，如只有均數(shù)標(biāo)準(zhǔn)差，而無原始數(shù)據(jù)時，無法直接調(diào)用 SAS軟件中的某個過程進(jìn)行分析，只能在數(shù)據(jù)步中按照計算公式編寫程序完成。 調(diào)查健康人群紅細(xì)胞數(shù) (1012/L)資料如下： 男性 360人，均數(shù) ，標(biāo)準(zhǔn)差 ； 女性 255人，均數(shù) ，標(biāo)準(zhǔn)差 。 試比較男女紅細(xì)胞數(shù)有無差別？ DATA BB。 INPUT N1 X1 S1 N2 X2 S2。 SE=SQRT(S1*S1/N1+S2*S2/N2)。 U=ABS(X1X2)/SE。 P1=1PROBNORM(U)。 P=2*P1。 U=ROUND(U,)。 P=ROUND(P,)。 IF P THEN P=。 FILE PRINT。 PUT 2 10 U= 25 P= 。 CARDS。 360 255 。 RUN。 定量資料假設(shè)檢驗(yàn) —— 非編程方式（分析員模塊） 菜單欄 Globals Analyze Statistical Analysis 菜單欄 File Open by SAS name (1)樣本均數(shù)與總體均數(shù)比較 菜單欄 Statistics Hypothesis Tasts OneSample ttest for Means 往 Variable窗口選入變量， Mean 窗口輸入 μ,OK ?μ和 σ己知時可選 OneSample ztest for Means Mean 窗口輸入 μ, Std Dev 窗口輸入 σ。 (2)配對 t檢驗(yàn) 菜單欄 Statistics Hypothesis Tasts TwoSample Paired ttest for Means 分別選入第一組和第二組變量 ， OK (3)兩樣本均數(shù)比較 t檢驗(yàn) 菜單欄 Statistics Hypothesis Tasts TwoSample ttest for Means 往 Dependent窗口選入要分析的變量 往 Group窗口選入分類變量， OK 四、單因素方差分析 用于完全隨機(jī)設(shè)計（成組設(shè)計）資料 例：比較三組人群載脂蛋白測定值有無差別？ 糖尿病 IGT異常 健康人 單因素方差分析用 ANOVA過程 DATA AA9。 INPUT X。 IF _N_12 THEN G=1。 ELSE IF _N_20 THEN G=3。 ELSE G=2。 CARDS。 …… ； PROC MEANS。 CLASS G。 VAR X。 PROC ANOVA。 CLASS G。 MODEL X=G。 MEANS G/SNK 。 RUN。 在 SNK后面加 HOVTEST可以作方差齊性檢驗(yàn) Analysis of Variance Procedure Dependent Variable: X Source DF Sum of Squares Mean Square F Value Pr F 變異來源 自由度