SAS 6.12教程

b · 发表于 2004-5-30 05:32:17

<

align=center>§4.1 GCHART过程<

align=justify>Gchart过程可以绘制水平/竖直条图（直方图）、立体直方图、饼图和星状图。我们可以用这些图来了解单个变量的分布或者多个变量之间的关系。<

>4.1.1 语法格式<DIV align=left><TABLE 3px outset; BORDER-TOP: 3px outset; BORDER-LEFT: 3px outset; BORDER-BOTTOM: 3px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top width=351>PROC GCHART [ DATA=<数据集名> [选项] ] ;</TD><TD vAlign=top width=239>指定要分析的数据集名及一些选项</TD></TR><TR><TD vAlign=top width=351>　HBAR <变量名列> / [选项] ;</TD><TD vAlign=top width=239>绘出条形图</TD></TR><TR><TD vAlign=top width=351>　VBAR <变量名列> / [选项];</TD><TD vAlign=top width=239>绘出水平条形图</TD></TR><TR><TD vAlign=top width=351>　BLOCK <变量名列> / [选项];</TD><TD vAlign=top width=239>绘出三维直方图</TD></TR><TR><TD vAlign=top width=351>　PIE <变量名列> / [选项];</TD><TD vAlign=top width=239>绘出饼图</TD></TR><TR><TD vAlign=top width=351>　STAR <变量名列> / [选项];</TD><TD vAlign=top width=239>绘出星状图</TD></TR><TR><TD vAlign=top width=351>　AXISn [选项];</TD><TD vAlign=top width=239>控制坐标轴的形状和颜色</TD></TR><TR><TD vAlign=top width=351>　BY <变量名列>;</TD><TD vAlign=top width=239>按该变量取值分层绘制，要求数据集已按该变量排序</TD></TR></TABLE></DIV><img src="http://medstatstar.myetang.com/sas/coach6/image11.gif">方括号中的为可选项，如果不写，系统会以默认值代替或者忽略该功能；而尖括号中的为必选项，表示如果写了该语句，则必须提供尖括号中的信息，否则可能出错。如[ DATA= <数据集名> ]，表示你可以不写这一项，此时系统默认处理的是最新建立的数据集；而如果写了DATA=这一项，你就必须指定一个数据集名，不然SAS跟你没完。4.1.2 语法说明事实上，绘图程序步中的选项极为繁杂，如果合理地使用这些选项可以做出非常漂亮的统计图来。但这些选项也给大家的学习带来很大的不便，因此在这里我们只简要介绍几个常用的绘图语句中的选项。【绘图语句选项】<UL><LI>MISSING 指定绘图时要将变量的缺失值也包括在内。 <LI>TYPE=做图类型关键字 指定要做图的类型，即图中条块代表的含义：缺省值是频数（FREQ）；如果指定了选择项SUMVAR，则缺省值为总和（SUM）。可选的关键字有： <LI>FREQ 要求按指定变量的频数做图。 <LI>PERCENT 要求按在横轴刻度表示范围内出现的频数占总数的百分比做图。 <LI>CFREQ 按累计频数做图。 <LI>CPERCENT 按累计百分比做图。 <LI>SUM 只能与SUMVAR选项同时使用，要求图中的每一条代表：变量在横轴表示的取值范围内时，SUMVAR指定变量的总和。 <LI>MEAN 只能与SUMVAR选项同时使用，要求图中的每一条代表：变量在横轴表示的取值范围内时，SUMVAR指定变量的均数。 <LI>SUMVAR=求和变量 指定使用TYPE=SUM或MEAN时，用于求总和、均值的变量。 <LI>LEVAL=n 如果绘图变量是连续变量，用该选项产生有N个组段的图形。 <LI>GROUP=分组变量 要求产生以分组变量的值分组的并排图。 <LI>SUBGROUP=亚组变量 要求每个图形内部再按亚组变量的值分块。 <LI>CAXIS=颜色 指定坐标轴的颜色。 <LI>CTEXT=颜色 指定坐标轴文本的颜色。
</LI></UL>4.1.3 应用实例例4.1 绘制《卫生统计学》第三版习题1.1的直方图，数据存在c:\user\wt1_1.dat中。解：程序如下：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD width=239>libname a 'c:\user'; </TD><TD width=343>指定c:\user文件夹为数据库a</TD></TR><TR><TD width=239>data a.wt1_1;</TD><TD width=343>数据步开始，指定要建立的数据集为a库的wt1_1</TD></TR><TR><TD width=239>　infile 'c:\user\wt1_1.dat';</TD><TD width=343>采用外部文件读入方式，文件名为c:\user\WT1_1.dat</TD></TR><TR><TD width=239>　input x @@; </TD><TD width=343>输入的变量为x，采用连续输入的格式</TD></TR><TR><TD width=239>proc gchart data=a.wt1_1;</TD><TD width=343>调用绘图程序步gchart</TD></TR><TR><TD width=239>　vbar x / levels=10;</TD><TD width=343>绘出直方图，用于绘图的变量为x，分为10组</TD></TR><TR><TD width=239>run; </TD><TD width=343>开始运行以上程序</TD></TR></TABLE><img src="http://medstatstar.myetang.com/sas/coach6/image10.gif">事实上ASSIST视窗的输出图形是彩色的，和我们用程序做出的有些区别，这是因为ASSIST视窗自动提交了一些绘图环境语句，对输出图形作了修饰的缘故。<img src="http://medstatstar.myetang.com/sas/coach6/image11.gif">实际上libname语句只需要运行一次就可以了，但在开始的几道例题里我们都写上了这一句，让大家逐步养成启动SAS后就运行这一句的习惯。请将C：\USER文件夹指定为A库，以后我们实习的数据集均放在该库中。

b · 发表于 2004-5-30 05:32:33

<

align=center>§4.2 GPLOT过程<

align=justify>GPLOT过程用于绘出散点图。图中的横、纵坐标分别代表两个变量。<

>4.2.1 语法格式<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top width=352>PROC GPLOT [ DATA=<数据集名> [选项] ] ;</TD><TD vAlign=top width=240>　</TD></TR><TR><TD vAlign=top width=352>　PLOT <纵坐标变量*横坐标变量[=分层变量名]...> / [选项];</TD><TD vAlign=top width=240>指定绘图变量和选项</TD></TR><TR><TD vAlign=top width=352>　PLOT2 <纵坐标变量*横坐标变量[=分层变量名]...> / [选项];</TD><TD vAlign=top width=240>在原图基础上重叠绘制第二幅散点图</TD></TR><TR><TD vAlign=top width=352>　SYMBOLn [选项]</TD><TD vAlign=top width=240>定义符号、添加趋势线、定义点和线的颜色</TD></TR><TR><TD vAlign=top width=352>　BY <变量名列>;</TD><TD vAlign=top width=240>按该变量取值分层绘制，要求数据集已按该变量排序</TD></TR></TABLE><img src="http://medstatstar.myetang.com/sas/coach6/image10.gif">实际上，SYMBOL语句中定义的选项为系统环境控制选项，这意味着一次定义，终生使用。除非重新定义，否则以后的输出图形都将会是第一个图形的“孪生兄弟”。<img src="http://medstatstar.myetang.com/sas/coach6/image11.gif">在这里SYMBOL后面紧跟了一个n，表示任意自然数。因为PLOT语句可以在同一坐标系内重叠绘制许多层图，而这个数字就表示SYMBOL语句是控制的哪一层图。4.2.2 语法说明【GPLOT过程的选项】<UL><LI>UNIFORM 要求用BY语句分组打印的散点图的坐标刻度相同，便于比较。 <LI>VTOH=数值 指定纵横坐标的比例。</LI></UL>【PLOT语句说明】<UL><LI>分层变量</LI></UL>表示所做的散点图按指定变量的取值分层，默认按不同颜色来区分。<UL><LI>语句选项 <LI>OVERLAY 同一语句做的图重叠在同一个坐标系中显示。 <LI>HAXIS=数值 定义横坐标的刻度。 <LI>VAXIS=数值 定义纵坐标的刻度。 <LI>CAXIS=颜色 定义坐标轴的颜色。 <LI>CTEXT=颜色 定义坐标轴文本的颜色。</LI></UL>【SYMBOL语句选项】<UL><LI>VALUE=符号 可用的符号及相应名称有： ＋ PLUS ★ STAR ■ SQUARE ◆ DIAMOND ▲ TRIANGLE <LI>I=连线方式 　① JOIN 用直线连接。</LI></UL>　　　　　　　　② SPLINE 用光滑的曲线连接。<OL><OL><OL><LI>NEEDLE 向横坐标画垂线。 <LI>RL 添加回归直线。</LI></OL></OL></OL><UL><LI>WIDTH=宽度 定义数据点和连线的宽度。 <LI>COLOR=颜色 定义数据点和连线的颜色。</LI></UL>4.2.3 应用实例例4.2 10名20岁男青年身高(cm)与前臂长(cm)如下，请绘出散点图（卫统p236 5.2题）。<DIV align=center><CENTER><TABLE cellSpacing=0 cellPadding=0 border=0><TR><TD vAlign=top>身 高</TD><TD vAlign=top>170</TD><TD vAlign=top>173</TD><TD vAlign=top>160</TD><TD vAlign=top>155</TD><TD vAlign=top>173</TD><TD vAlign=top>168</TD><TD vAlign=top>178</TD><TD vAlign=top>183</TD><TD vAlign=top>180</TD><TD vAlign=top>165</TD></TR><TR><TD vAlign=top>前臂长</TD><TD vAlign=top>45</TD><TD vAlign=top>42</TD><TD vAlign=top>44</TD><TD vAlign=top>41</TD><TD vAlign=top>47</TD><TD vAlign=top>50</TD><TD vAlign=top>47</TD><TD vAlign=top>46</TD><TD vAlign=top>49</TD><TD vAlign=top>43</TD></TR></TABLE></CENTER></DIV>解：程序如下（注意程序中对输出图形作了些修饰）。<TABLE 3px outset; BORDER-TOP: 3px outset; BORDER-LEFT: 3px outset; BORDER-BOTTOM: 3px outset" cellSpacing=0 cellPadding=0 width=603 border=0><TR><TD width=298>Libname a ’c:\user’; </TD><TD width=295>指定c:\user文件夹为数据库a</TD></TR><TR><TD width=298>data a.wt5_2;</TD><TD width=295>数据步开始，指定要建立的数据集为a库的wt5_2</TD></TR><TR><TD width=298>　input x y @@;</TD><TD width=295>输入的变量为x和y，采用连续输入的格式</TD></TR><TR><TD width=298>　cards; </TD><TD width=295>数据块开始</TD></TR><TR><TD width=298>　170 45 173 42 160 44 155 41 173 41 </TD><TD width=295>数据块</TD></TR><TR><TD width=298>　188 50 178 47 183 46 180 49 165 43 </TD><TD width=295>数据块</TD></TR><TR><TD width=298>; </TD><TD width=295>数据块结束</TD></TR><TR><TD width=298>proc gplot data=a.wt5_2; </TD><TD width=295>调用绘图程序步gplot</TD></TR><TR><TD width=298>　plot y*x/ ctext=blue;</TD><TD width=295>绘出散点图，用于绘图的变量为x，y,坐标轴文本颜色为蓝色</TD></TR><TR><TD width=298>　symbol value=star color=pink width=2; </TD><TD width=295>定义散点标记为粉红色星号，大小为2</TD></TR><TR><TD width=298>run; </TD><TD width=295>开始运行以上程序</TD></TR></TABLE><img src="http://medstatstar.myetang.com/sas/coach6/image26.gif"> 绘制高分辨率统计图非常耗费系统资源，因此SAS在执行作图语句时实际上是将这些语句存入程序缓冲区，当用户打开GRAPH视窗看结果时才正式执行，从而绘出所需的统计图。如果用户没有看所绘制的图形，则绘图语句就一直留在缓冲区内，此时用户再执行任何程序，程序语句就全部在缓冲区内排队，等待绘图语句执行完毕，因此就无法得到相应的输出结果，这是同学们在上机时常犯的一个错误。适当地使用统计图可以帮助我们更加深入的了解数据的分布规律。这里的例子只是绘图语句的简单用法，而绘图语句的复杂用法将结合例题在以后的各章中逐步介绍。

b · 发表于 2004-5-30 05:33:23

<

align=center> <

align=center>第五章 定量资料的统计描述和t、u检验 <

align=center> 从本章开始，我们将正式开始使用SAS解决我们的统计问题。从前面的几章可知，SAS的主要功能是由不同的程序步来体现的。因此在以后的各章中，我们将对每种问题重点介绍一些常用的程序步，以及它们的输出结果的解释。对于定量资料的统计描述和简单推断，SAS提供了三个强有力的程序步，它们是：<UL><LI>UNIVARIATE过程 提供单个变量的详细描述和对其分布类型的检验。 <LI>MEANS过程 提供单个或多个变量的简单描述，对于多个变量，它的输出格式紧凑，便于阅读。 <LI>TTEST过程对变量进行t/u检验。</LI></UL>

b · 发表于 2004-5-30 05:33:43

<

align=center>§5.1 引 例<

align=justify>例5.1 文本文件“C:\USER\WTLI1_1.DAT”中已存入某市110名7岁男子童的身高资料(cm) ，请计算均数<img src="http://medstatstar.myetang.com/sas/coach6/image29.gif">、标准差s、变异系数CV（卫统第三版p6例2.1）。<

align=justify>解：该题应首先用数据步建立一个数据集，然后调用UNIVARIATE过程或MEANS过程来求出所需要的统计量。具体的程序如下：① 设定数据库环境：LIBNAME A 'C:\USER';② 数据步，建立数据集：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 border=0><TR><TD vAlign=top>DATA A.WTLI1_1;</TD></TR><TR><TD vAlign=top>　INFILE 'C:\USER\WTLI1_1.DAT';</TD></TR><TR><TD vAlign=top>　INPUT X @@;</TD></TR><TR><TD vAlign=top>RUN;</TD></TR></TABLE>③ UNIVARIATE或MEANS过程，求出所需要的统计量：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top>PROC UNIVARIATE DATA=A.WTLI1_1;</TD><TD vAlign=top>PROC MEANS DATA=A.WTLI1_1N MEAN STD CV ;</TD></TR><TR><TD vAlign=top width="52%">　VAR X;</TD><TD vAlign=top width="48%">　VAR X;</TD></TR><TR><TD vAlign=top width="52%">RUN;</TD><TD vAlign=top width="48%">RUN;</TD></TR></TABLE>例5.2 某医生测得18例慢支炎患者与16例健康人的尿17酮类固醇排出量 (mg/dl)分别为X1和X2，问两者均数是否不同（医统第二版P19例2.17）？解：这是成组设计的两样本均数比较的t检验，程序应首先建立数据集，然后调用TTEST过程进行检验，在检验的同时也可以得到两个样本的简单描述。① 设定数据库环境：LIBNAME A 'C:\USER';② 数据步，建立数据集,这里采用直接输入数据的方法：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 border=0><TR><TD vAlign=top>DATA A.YTLI2_17;</TD></TR><TR><TD vAlign=top>　INPUT GROUP VALUE @@;</TD></TR><TR><TD vAlign=top>　CARDS;</TD></TR><TR><TD vAlign=top>　1 3.14 1 5.83 1 7.35 1 4.62 1 4.05 1 5.08 1 4.98 1 4.22 1 4.35 1 2.35</TD></TR><TR><TD vAlign=top>　1 2.89 1 2.16 1 5.55 1 5.94 1 4.4 1 5.35 1 3.8 1 4.12</TD></TR><TR><TD vAlign=top>　2 4.12 2 7.89 2 3.24 2 6.36 2 3.48 2 6.74 2 4.67 2 7.38 2 4.95 2 4.08</TD></TR><TR><TD vAlign=top>　2 5.34 2 4.27 2 6.54 2 4.62 2 5.92 2 5.18</TD></TR><TR><TD vAlign=top>　;</TD></TR><TR><TD vAlign=top>RUN;</TD></TR></TABLE>③ TTEST过程，进行两样本的t检验。<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 border=0><TR><TD vAlign=top>PROC TTEST DATA=A.YTLI2_17;</TD></TR><TR><TD vAlign=top>　VAR VALUE;</TD></TR><TR><TD vAlign=top>　CLASS GROUP;</TD></TR><TR><TD vAlign=top>RUN;</TD></TR></TABLE>

b · 发表于 2004-5-30 05:33:58

<

align=center>§5.2 UNIVARIATE过程<

align=justify>Univariate过程对数值变量给出比较详细的变量分布的描述，其中包括：<UL><LI>变量的极端值。 <LI>常用的百分位数，包括四分位数和中位数。 <LI>用几个散点图描绘变量的分布。 <LI>频数表。 <LI>确定数据为正态分布的检验。</LI></UL><

>5.2.1 语法格式Univariate过程的语法格式如下：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top width="50%">PROC UNIVARIATE [ DATA= <数据集名> [选项] ];</TD><TD vAlign=top width="50%">指定要分析的数据集名及选项</TD></TR><TR><TD vAlign=top width="50%">　[ VAR <变量名列> ;</TD><TD vAlign=top width="50%">指定要分析的变量名列</TD></TR><TR><TD vAlign=top width="50%">　BY <变量名列> ;</TD><TD vAlign=top width="50%">按变量名列分组统计，要求数据集已按该变量名列排序</TD></TR><TR><TD vAlign=top width="50%">　FREQ <变量名> ;</TD><TD vAlign=top width="50%">表明该变量为分析变量的频数</TD></TR><TR><TD vAlign=top width="50%">　WEIGHT <变量名> ;</TD><TD vAlign=top width="50%">表明分析变量在统计时要按该变量权重</TD></TR><TR><TD vAlign=top width="50%">　ID <变量名> ;</TD><TD vAlign=top width="50%">输出时加上该变量作为索引</TD></TR><TR><TD vAlign=top width="50%">　OUTPUT OUT= <数据集名></TD><TD vAlign=top width="50%">指定统计量的输出数据集名</TD></TR><TR><TD vAlign=top width="50%">　关键字= <新变量名列>...</TD><TD vAlign=top width="50%">指定统计量对应的新变量名</TD></TR><TR><TD vAlign=top width="50%">　pctlpts=<百分位数, ...></TD><TD vAlign=top width="50%">指定需要的百分位数</TD></TR><TR><TD vAlign=top width="50%">　pctlpre=<新变量名列>] ;</TD><TD vAlign=top width="50%">指定所需百分位数对应的输出变量名</TD></TR></TABLE>如果省略所有非必需的语句和选项，则UNIVARIATE过程按默认情况输出全部变量的全部常用统计量。

b · 发表于 2004-5-30 05:34:11

<

>5.2.2 语法说明<

align=justify>【选项】<

align=justify>Univariate过程常用的选项如下：<UL><LI>NOPRINT 禁止统计报告在OUTPUT视窗中输出 <LI>PLOT 绘出茎叶图、箱式图和正态概率图 <LI>FREQ 给出频数表 <LI>NORMAL 对变量进行正态性检验</LI></UL>【关键字】SAS中用关键字来指定所需要的统计量，事实上结果输出中用的就是各种关键字，常用的关键字有：<UL><LI>基本统计量 N MEAN STD(标准误) CV SUM VAR(方差) RANG <LI>百分位数描述 MIN P1 P5 P10 Q1 MEDIAN Q3 P90 P95 P99 MAX <LI>与假设检验有关的统计量 STDMEAN(标准误) T</LI></UL>

b · 发表于 2004-5-30 05:34:22

<

>5.2.3 结果解释<

align=justify>在默认的情况时，Univariate过程会输出绝大部分统计量，此时的输出结果如下：<

RE>Variable=变量名变量标签</PRE><PRE> Moments 和矩有关的统计量 Quantiles(Def=5) 分位间距统计量</PRE><PRE>N 样本量 Sum Wgts 权重总和 100% Max 最大值 99% 99%百分位数Mean 均数 Sum 总和 75% Q3 75%百分位数 95% 95%百分位数Std Dev 标准差 Variance 方差 50% Med 50%百分位数 90% 90%百分位数Skewness 偏度系数 Kurtosis 峰度系数 25% Q1 25%百分位数 10% 10%百分位数USS 未校正平方和 CSS 校正平方和 0% Min 最小值 5% 5%百分位数CV 变异系数 Std Mean 标准化均数 1% 1%百分位数 T:Mean=0 变量总体均数为0的t检验 Pr>|T| t检验的p值 Range 全距Num ^= 0 变量值非0的例数 Num > 0 变量值大于0的例数 Q3-Q1 四分位间距M(Sign) 变量总体均数为0的符号检验 Pr>=|M| 符号检验的p值 Mode 众数Sgn Rank 变量总体均数为0的秩和检验 Pr>=|S| 秩和检验的p值</PRE><PRE> Extremes 极端值统计</PRE><PRE> Lowest Obs Highest Obs 老幺 (观察值序号) 五大值(观察值序号) 次小值 (观察值序号) 四大值(观察值序号) 三小值 (观察值序号) 三大值(观察值序号) 四小值 (观察值序号) 次大值(观察值序号) 五小值 (观察值序号) 大哥大(观察值序号)</PRE>

b · 发表于 2004-5-30 05:35:13

<

>5.2.4 应用实例<

align=justify>例5.3 某地101例健康男子血清总胆固醇值测定结果已存入文本文件“c:\user\WT1_1.dat”中，请绘制直方图，计算均数<img src="http://medstatstar.myetang.com/sas/coach6/image29.gif">、标准差s、变异系数CV、中位数M、p2.5和p97.5（卫统p233 1.1题）。<

align=justify>解：UNIVARIATE过程的默认输出中并不给出p2.5和p97.5，因此程序中要加以相应修改，最后在OUTPUT视窗中只会输出所需的几个统计量，具体程序如下：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top>libname a 'c:\user';</TD><TD vAlign=top>指定c:\user文件夹为数据库a</TD></TR><TR><TD vAlign=top>data a.wt1_1;</TD><TD vAlign=top>数据步开始，指定要建立的数据集为a库的wt1_1</TD></TR><TR><TD vAlign=top>　infile 'c:\user\wt1_1.dat';</TD><TD vAlign=top>采用外部文件读入方式，文件名为c:\user\WT1_1.dat</TD></TR><TR><TD vAlign=top>　input x @@;</TD><TD vAlign=top>输入的变量为x，采用连续输入的格式</TD></TR><TR><TD vAlign=top>proc gchart data=a.wt1_1;</TD><TD vAlign=top>调用绘图程序步gchart，所用数据集为a.wt1_1</TD></TR><TR><TD vAlign=top>　vbar x ;</TD><TD vAlign=top>绘出竖直条图，用于绘图的变量为x</TD></TR><TR><TD vAlign=top>proc univariate data=a.wt1_1 noprint;</TD><TD vAlign=top>调用程序步univariate，并且禁止在OUTPUT视窗中输出</TD></TR><TR><TD vAlign=top>　var x;</TD><TD vAlign=top>要分析的变量为x</TD></TR><TR><TD vAlign=top>　output out=temp</TD><TD vAlign=top>指定输出数据集为work.temp,</TD></TR><TR><TD vAlign=top>　n=n mean=xbar std=s cv=cv median=m</TD><TD vAlign=top>将n、mean、std、cv、median按指定变量名存入</TD></TR><TR><TD vAlign=top>　pctlpts=2.5,97.5 pctlpre=per;</TD><TD vAlign=top>指定输出p2.5和p97.5，其输出变量名以per开头。</TD></TR><TR><TD vAlign=top>proc print data=temp;</TD><TD vAlign=top>将数据集work.temp的内容打印输出</TD></TR><TR><TD vAlign=top>run;</TD><TD vAlign=top>开始运行以上程序</TD></TR></TABLE>

b · 发表于 2004-5-30 05:35:33

<

align=justify>例5.4 50例链球菌咽峡炎患者的潜伏期如下，计算其均数、中位数和几何均数（卫统p233 1.3题）。<DIV align=center><CENTER><TABLE cellSpacing=0 border=0><TR><TD vAlign=top align=middle width=50><

align=center>12~</TD><TD vAlign=top align=middle width=50><

align=center>24~</TD><TD vAlign=top align=middle width=50>36~</TD><TD vAlign=top align=middle width=50>48~</TD><TD vAlign=top align=middle width=50>60~</TD><TD vAlign=top align=middle width=50>72~</TD><TD vAlign=top align=middle width=50>84~</TD><TD vAlign=top align=middle width=50>96~</TD><TD vAlign=top align=middle width=50>108~120</TD></TR><TR><TD vAlign=top align=middle width=50>1</TD><TD vAlign=top align=middle width=50>7</TD><TD vAlign=top align=middle width=50>11</TD><TD vAlign=top align=middle width=50>11</TD><TD vAlign=top align=middle width=50>7</TD><TD vAlign=top align=middle width=50>5</TD><TD vAlign=top align=middle width=50>4</TD><TD vAlign=top align=middle width=50>2</TD><TD vAlign=top align=middle width=50>2</TD></TR></TABLE></CENTER></DIV>解：由于几何均数无法直接得到，因此将数据集加以对数变换，求出均数后再行反对数变换得到几何均数，程序如下：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD 0px" vAlign=top>libname a 'c:\user';</TD><TD 0px" vAlign=top>指定c:\user文件夹为数据库a</TD></TR><TR><TD 0px" vAlign=top>data a.wt1_3;</TD><TD 0px" vAlign=top>数据步开始，指定要建立的数据集为a库的wt1_3</TD></TR><TR><TD 0px" vAlign=top>　input x f @@;</TD><TD 0px" vAlign=top>输入的变量为x和f，采用连续输入的格式</TD></TR><TR><TD 0px" vAlign=top>　x=x+6;</TD><TD 0px" vAlign=top>将变量x的值更正到每个组段的组中值处</TD></TR><TR><TD 0px" vAlign=top>　logx=log(x);</TD><TD 0px" vAlign=top>定义新变量logx为变量x的自然对数，用于算出几何均数</TD></TR><TR><TD 0px" vAlign=top>　cards;</TD><TD 0px" vAlign=top>数据块开始</TD></TR><TR><TD 0px" vAlign=top>　12 1 24 7 36 11 48 11 60 7 72 5 84 4 96 2 106 2</TD><TD 0px" vAlign=top>数据块</TD></TR><TR><TD 0px" vAlign=top>;</TD><TD 0px" vAlign=top>数据块结束</TD></TR><TR><TD 0px" vAlign=top>proc print;</TD><TD 0px" vAlign=top>将数据集a.wt1_3的内容打印输出</TD></TR><TR><TD 0px" vAlign=top>proc univariate data=a.wt1_3 noprint;</TD><TD 0px" vAlign=top>调用程序步univariate，并且禁止在OUTPUT视窗中输出</TD></TR><TR><TD 0px" vAlign=top>　var x logx;</TD><TD 0px" vAlign=top>要分析的变量为x和logx</TD></TR><TR><TD 0px" vAlign=top>　freq f;</TD><TD 0px" vAlign=top>指定变量f代表分析变量x的频数</TD></TR><TR><TD 0px" vAlign=top>　output out=temp n=n mean=xbar logxmean median=m ;</TD><TD 0px" vAlign=top>输出数据集和统计量的定义</TD></TR><TR><TD 0px" vAlign=top>data temp2;</TD><TD 0px" vAlign=top>数据步开始，指定要建立的数据集为work.temp2</TD></TR><TR><TD 0px" vAlign=top>　set temp;</TD><TD 0px" vAlign=top>让work.temp2继承work.tmep的全部数据</TD></TR><TR><TD 0px" vAlign=top>　g=exp(logxmean);</TD><TD 0px" vAlign=top>产生新变量g，它等于elogxmean</TD></TR><TR><TD 0px" vAlign=top>　drop logxmean;</TD><TD 0px" vAlign=top>在work.temp2中删除临时变量logxmean</TD></TR><TR><TD 0px" vAlign=top>proc print data=temp2;</TD><TD 0px" vAlign=top>输出数据集work.temp2中的数据</TD></TR><TR><TD 0px" vAlign=top>run;</TD><TD 0px" vAlign=top>开始运行以上程序</TD></TR></TABLE><img src="http://medstatstar.myetang.com/sas/coach6/image11.gif">请注意，这里UNIVARIATE过程有两个分析变量X和LOGX，因此在OUTPUT语句中MEAN=后跟了两个变量名，它们分别存储两个变量的均数，而其余的关键字后只有一个变量名，则它们只存储分析变量序列的第一个变量X的统计结果。

b · 发表于 2004-5-30 05:35:47

<

align=center>§5.3 MEANS过程<

align=justify>Means过程提供单个或多个变量的简单描述。和Univariate过程相比，它更倾向于描述已经明确样本所在总体符合正态分布的变量，因此它不提供百分位数，但可以提供95%可信区间。同时在多个变量输出时，它的输出格式紧凑，便于阅读。<

>5.3.1 语法格式<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 width=600 border=0><TR><TD vAlign=top width="50%">PROC MEANS [ DATA= <数据集名> [选项]</TD><TD vAlign=top width="50%">指定要分析的数据集名及一些选项</TD></TR><TR><TD vAlign=top width="50%">　[统计量关键字列表] ] ;</TD><TD vAlign=top width="50%">列出需要的统计量</TD></TR><TR><TD vAlign=top width="50%">　[VAR <变量名列>;</TD><TD vAlign=top width="50%">要分析的变量名列</TD></TR><TR><TD vAlign=top width="50%">　BY <变量名列>;</TD><TD vAlign=top width="50%">按变量名列分组统计，要求数据集已按变量名列排序</TD></TR><TR><TD vAlign=top width="50%">　CLASS <变量名列>;</TD><TD vAlign=top width="50%">按变量名列分组统计，不要求数据集排序</TD></TR><TR><TD vAlign=top width="50%">　FREQ <变量名>;</TD><TD vAlign=top width="50%">表明该变量为分析变量的频数</TD></TR><TR><TD vAlign=top width="50%">　WEIGHT <变量名>;</TD><TD vAlign=top width="50%">表明分析变量在统计时要按该变量权重</TD></TR><TR><TD vAlign=top width="50%">　ID <变量名列>;</TD><TD vAlign=top width="50%">输出时加上该变量作为索引</TD></TR><TR><TD vAlign=top width="50%">　OUTPUT <OUT= 数据集名></TD><TD vAlign=top width="50%">指定统计量的输出数据集名</TD></TR><TR><TD vAlign=top width="50%">　关键字= <新变量名列>... ] ;</TD><TD vAlign=top width="50%">指定统计量对应的新变量名</TD></TR></TABLE>5.3.2 语法说明【选项】Means过程常用的选项如下：<UL><LI>NOPRINT 禁止统计报告在OUTPUT视窗中输出 <LI>MAXDEC=n 给出列表输出的最大小数位数，缺省值为2</LI></UL>【统计量关键字】MEANS过程中常用的统计量关键字有：<UL><LI>基本统计量 N MEAN STD CV SUM VAR RANG MIN MAX <LI>与假设检验有关的统计量 STDERR(标准误） T PRT(与t对应的p值) LCLM(可信区间下限) UCLM(可信区间上限)</LI></UL><img src="http://medstatstar.myetang.com/sas/coach6/image11.gif">注意Means过程中标准误的关键字是STDERR，而Univariate过程中为STDMEAN，另外LCLM和UCLM这两个关键字也是Univariate过程所没有的。5.3.3 结果解释和Univariate过程不同，MEANS过程在默认情况下只输出样本量、均数、标准差、最小值和最大值，如例5.1的数据，如果MEANS过程不加任何选项，则输出如下：<PRE> Analysis Variable : X 分析变量名为X</PRE><PRE> N Mean Std Dev Minimum Maximum --------------------------------------------------------------- 110 119.7272727 4.7413254 108.2000000 132.5000000 ---------------------------------------------------------------</PRE>可见Means过程的输出结构极为紧凑。5.3.4 应用实例例5.5 给出例5.1的均数<img src="http://medstatstar.myetang.com/sas/coach6/image29.gif">、标准差s、变异系数CV和95%可信区间，并检验其总体均数是否为0（卫统p233 1.1题）。解：如果数据集a.wt1_1已经建立，则程序如下：<TABLE 2px outset; BORDER-TOP: 2px outset; BORDER-LEFT: 2px outset; BORDER-BOTTOM: 2px outset" cellSpacing=0 cellPadding=0 border=0><TR><TD vAlign=top>proc means data=a.wt1_1 n mean std cv lclm uclm t prt ;</TD></TR><TR><TD vAlign=top>　var x;</TD></TR><TR><TD vAlign=top>run;</TD></TR></TABLE>

		自动登录	找回密码
密码			注-册-帐-号