實(shí)習(xí)二定量資料的統(tǒng)計(jì)描述1

定量資料的統(tǒng)計(jì)描述

SAS中的univariate過程、means過程、summary過程tabulate過程可用于定量資料的統(tǒng)計(jì)描述。

相同點(diǎn)

不同點(diǎn)一、UNIVARIATE過程基本格式

procunivariate選項(xiàng)；

by變量名稱（分組變量）；

class變量名稱（分組變量）；

freq變量名稱（數(shù)值變量，用以表示相應(yīng)記錄出現(xiàn)的頻數(shù)）

weight變量名稱（數(shù)值變量，用以表示相應(yīng)記錄的權(quán)重系數(shù)）

histogram變量名稱/繪圖選項(xiàng)

output<out=數(shù)據(jù)集名><統(tǒng)計(jì)量關(guān)鍵字=自定義變量名>

var變量名稱（待分析的數(shù)值變量）；run；選擇項(xiàng)DATA=

指定要分析的數(shù)據(jù)集名NOPRINT

禁止統(tǒng)計(jì)結(jié)果在OUTPUT視窗中輸出NORMAL

對變量進(jìn)行正態(tài)性檢驗(yàn)PLOT

生成統(tǒng)計(jì)圖

注：在proc語句中的選項(xiàng)直接跟在proc名后，以空格分隔，中間不加“/”，其他命令語句的選項(xiàng)需以“/”分隔。二、MEANS過程的基本格式

procmeans選項(xiàng)列表；

by變量名稱（分組變量）；

class變量名稱（分組變量）；

freq變量名稱（數(shù)值變量）

weight變量名稱（數(shù)值變量）

var變量名稱（待分析的數(shù)值變量）；run；選擇項(xiàng)data=

指定要分析的數(shù)據(jù)集名maxdec=

指定輸出結(jié)果的小數(shù)位數(shù),默認(rèn)為7位noprint

禁止統(tǒng)計(jì)結(jié)果在output窗輸出alpha=

設(shè)定計(jì)算可信區(qū)間的α水平，默認(rèn)為0.05當(dāng)無選擇項(xiàng)時(shí)，MEANS過程默認(rèn)為每個(gè)變量輸出五個(gè)統(tǒng)計(jì)值：N、MIN、MAX、MEAN、STD。OUTPUT<OUT=數(shù)據(jù)集名>[關(guān)鍵字=<新變量名列表>]

/*將分析結(jié)果輸出到SAS數(shù)據(jù)集*/OUT=數(shù)據(jù)集名

/*指定SAS數(shù)據(jù)集名稱*/關(guān)鍵字=<新變量名列表>

/*指定要輸出的統(tǒng)計(jì)量*/

可輸出的有N,MEAN,STD,MIN，MAX，RANGE,SUM,VAR,USS(平方和),CSS(離均差平方和)，CV,STDERR,T,PRT,SKEWNESS,KURTOSIS,LCLM,UCLM（同時(shí)指定輸出為雙側(cè)；如只選一個(gè)，輸出單側(cè)可信區(qū)間）。

可以重新對統(tǒng)計(jì)量命名，也可用原名：格式n=nmean=meanprt=p三、SUMMARY過程的基本格式

procsummary選項(xiàng)列表；

by變量名稱（分組變量）；

class變量名稱（分組變量）；

freq變量名稱（數(shù)值變量）

weight變量名稱（數(shù)值變量）

output<out=數(shù)據(jù)集名><統(tǒng)計(jì)量關(guān)鍵字=自定義變量名>

var變量名稱（待分析的數(shù)值變量）；run；四、TABULATE過程的基本格式

proctabulate選項(xiàng)列表；

by變量名稱（分組變量）；

class變量名稱（分組變量）；

freq變量名稱（數(shù)值變量）

weight變量名稱（數(shù)值變量）

table<<頁變量表達(dá)式>，<行變量表達(dá)式>，<列變量表達(dá)式>></表格選項(xiàng)>

var變量名稱（待分析的數(shù)值變量，統(tǒng)計(jì)量列入相應(yīng)的表單元格）；run；VAR語句如果只想分析某一個(gè)或幾個(gè)特定的變量，則可用VAR語句指定它們。例如只想顯示變量x的列表，則PRINT過程如下：

procprint;

varx;run;如果想顯示x和y兩個(gè)變量，則將VAR語句改為如下形式即可：

varxy;

Where語句如果你想處理的不是整個(gè)數(shù)據(jù)集而只是其中符合某種條件的子集，可以應(yīng)用WHERE語句，如上例中我們只想顯示大于50的x變量的值，則在PRINT過程中加入where語句如下：

wherex>50;如果條件變?yōu)閤、y兩個(gè)變量的值都要大于50，則where語句改為：wherex>50andy>50;

By語句如果你需要分組處理數(shù)據(jù)，例如要按性別分組輸出統(tǒng)計(jì)結(jié)果，可以采用by語句。如上例我們想按不同的y值輸出x值，則加入by語句如下：

byy;使用BY語句要求數(shù)據(jù)集已經(jīng)按BY語句中指定的變量排序。如果沒有排序，則程序無法正確運(yùn)行?？梢杂肧ORT過程來排序，語法結(jié)構(gòu)如下：

PROCSORTDATA=數(shù)據(jù)集名；

BY變量名列；RUN；

相同點(diǎn)：他們均可計(jì)算出均數(shù)、標(biāo)準(zhǔn)差、方差、標(biāo)準(zhǔn)誤、總和、加權(quán)值的總和、最大值、最小值、全距、校正的和未校正的離差平方和、變異系數(shù)、樣本分布位置的t檢驗(yàn)統(tǒng)計(jì)量、遺漏數(shù)據(jù)和有效數(shù)據(jù)個(gè)數(shù)等，均可應(yīng)用by語句將樣本分割為若干個(gè)更小的樣本，以便分別進(jìn)行分析。不同點(diǎn)：①means、summary、univariate過程可以計(jì)算樣本的偏度（skewness）和峰度（kurtosis），而tabulate過程不計(jì)算；②univariate過程可以計(jì)算出樣本的眾數(shù)（mode）；③summary過程執(zhí)行后不會(huì)自動(dòng)給出分析的結(jié)果，須引用output語句和print過程來顯示分析結(jié)果；④univariate過程具有統(tǒng)計(jì)制圖的功能。dataex2;inputxf@@;y=log(x);cards;1122.51410184022160126406;proc

means;vary;freqf;outputout=ex2gmean=;dataex2a;setex2g;g=exp(y);proc

printdata=ex2a;varg;run;P20例2-3datahu1;inputx@@;cards;119.3121.2116.6126.4120.0115.5119.9116.5121.6116.9122.1123.4119.2117.3114.4119.9118.7123.2119.8122.3122.6127.6110.4115.9129.4115.2115.9113.0126.3125.3114.6122.4125.7124.0130.4118.9122.1118.6117.2116.7121.3115.6120.5125.2122.2112.8117.3125.3116.7132.2119.3116.4118.2121.7118.6115.2118.0109.7120.5110.3117.2113.7120.0117.0121.3112.9110.2116.1118.5120.0120.7121.4119.6114.2121.7123.6120.6112.7119.8123.4122.3113.2120.1124.3120.5127.8123.7125.9108.9119.7114.3127.8120.3114.5130.6120.8120.5114.6119.1114.5121.9;proc

meansmaxdec=2meanstdminmax;varx;run;P20例2-3datahu1;inputx@@;cards;119.3121.2116.6126.4120.0115.5119.9116.5121.6116.9122.1123.4;proc

univariate;outputout=abn=nmean=xbarmedian=m;varx;histogramx;proc

printdata=ab;run;P362