第四章 多组的资讯均数比较的SAS编程实现

2022-02-07 02:02:40 来源:
分享:
各种类型资讯均多达较为我们一般选用方差系统性的方法,SAS之前方差系统性的机制极为全盘,能付诸方差系统性机制的流程有anova流程和glm流程。

一、anova流程和glm流程简介

1. anova流程

anova流程依赖于于stat模块,主要的机制就是顺利进行方差系统性。anova流程意在对平衡实验新设计资讯(各分一组考量各素质的所有两一组具有不同的采样量或观察个数)顺利进行系统性,不能用于对非平衡实验新设计资讯的方差系统性。它比glm流程的行驶反应速度要快,尽快的存贮空间也要小一些。

anova流程的一般播放器如下:

PROC ANOVA < 的系统沙罗> ;

CLASS 给定名沙罗 ;

MODEL 理应给定=自给定运算符 < / 的系统沙罗> ;

ABSORB 给定沙罗 ;

BY 给定名 ;

FREQ 给定名 ;

MANOVA < 验的系统>< / 详细的系统> ;

MEANS 自给定运算符 < / 的系统> ;

REPEATED 给定名 素质多达 <(各素质个数)>< / 的系统> ;

TEST < H=给定运算符> E=给定运算符;

proc anova运算符启动anova流程的行驶,年初的的系统沙罗可成份表列计划:

DATA=多达据集名

MANOVA

MULTIPASS

NAMELEN=n

NOPRINT

ORDER=DATA | FORMATTED | FREQ | INTERNAL

OUTSTAT=SAS-data-set

其之前的“order=”的系统原则上anova流程对归纳给定(class运算符原则上的给定)各素质的以此类推方式,可选的四个条目分别问到按照原始多达据之前的以此类推、控制器播放器个数的以此类推、各素质观察个数频多达多少的以此类推、内部个数的以此类推顺利进行以此类推。其余的系统和过去的具体内容不同或较少加到,这里才会行先为前解释。

年初的class运算符、by运算符、freq运算符等和过去的具体内容只不过不同,我们这里着力知晓一下anova流程之前特有的决定性运算符。

model运算符:用来原则上理应给定和自给定,并且通过特定的运算符规定自给定的作用方式。如果从未原则上任何自给定,则假设之前极为少包涵常多达项,此时验的具体内容是理应给定的均多达是否为零。model运算符之前原则上的自给定需要是class运算符之前问到遗憾过的归纳给定,anova流程不允许自给定之前有连续型给定(最大个数给定),而理应给定则需要是最大个数型给定。

自给定运算符最简单来表达三种有所不同的效理应假设,即主效理应假设、交互效理应假设、结构体新设计效理应假设。我们假定理应给定为y,a, b, c为假设之前包涵的三个自给定,各种效理应假设的表达方式如下:

(1)主效理应假设:y=a b c;

(2)交互效理应假设:y=a b c a*b a*c b*c a*b*c;

(3)结构体新设计假设:y=a b c(a b); 其之前c考量为a、b两考量各种两一组下的二级考量。

model运算符末尾的的系统有“intercept”和“nouni”两项,分别原则上SAS顺利进行关于常多达项的假设验和在多给定方差系统性(或单调量度资讯方差系统性)时禁止单给定统计结果的控制器。

absorb运算符:对于极为少发挥主效理应的考量,absorb运算符原则上SAS消除此给定的作用,只对其余给定顺利进行系统性,起到略为减少量度机自然资源和等待时间消耗的作用。理技术的发展此机制时,仍须先为对原则上给定以此类推,且此给定不能再行注意到在class运算符和model运算符之前。

manova运算符:当反理应给定有多个时,此运算符控制anova流程进入多元方差系统性模式,年初的的系统意在原则上多元方差系统性时的各项指标。

means运算符:原则上anova流程量度原则上给定各素质下反理应给定的均个数、均个数,并顺利进行一组间的多重较为。

repeated运算符:如果反理应给定为单调量度多达据,此运算符意在原则上anova流程进入单调量度多达据方差系统性模式。其之前的给定名都是单调量度考量(如量度等待时间等),年初素质多达都是单调量度的次多达,如果必需原则上单调量度各次的具体标签,可在年初按以此类推列出,并用圆括号括上去。

test运算符:意在顺利进行其他类型的f验,这种验有所不同于并不一定方差系统性之前以标准差均均是由小多达的f验,其之前的的系统“h=”意在原则上作为小分子的给定运算符(需要在model运算符之前注意到过),“e=”意在原则上一个作为小多达的给定。

2. glm流程

glm流程也依赖于于stat模块之前,它执行以很小二加法顺利进行假设拟合的机制。意在流程可以付诸的统计学方法有回归系统性、方差系统性、MSE系统性、多元方差系统性以及偏相关系统性。glm流程对多达据的系统性处理均在一般线性假设的框架下顺利进行,反理应给定可以为一个或多个连续型给定,自给定可为连续型也可为离散型。

glm流程的一般播放器如下:

PROC GLM < 的系统沙罗> ;

CLASS 给定名沙罗 ;

MODEL 理应给定=自给定沙罗 < / 的系统沙罗> ;

ABSORB 给定名沙罗 ;

BY 给定名沙罗 ;

FREQ 给定名沙罗 ;

ID 给定名沙罗 ;

WEIGHT 给定名沙罗 ;

CONTRAST '上面' 效理应运算符 常多达乘积 < ... 效理应运算符 常多达乘积> < / 的系统沙罗> ;

ESTIMATE '上面' 效理应运算符 常多达乘积 < ... 效理应运算符 常多达乘积> < / 的系统沙罗> ;

LSMEANS效理应运算符< / 的系统沙罗> ;

MANOVA < 验的系统>< / 详细的系统> ;

MEANS 效理应运算符 < / 的系统沙罗> ;

OUTPUT < OUT=多达据集名> keyword=names < / 的系统> ;

RANDOM效理应运算符 < / 的系统沙罗> ;

REPEATED 考量运算符 < / 的系统沙罗> ;

TEST < H=给定运算符> E=给定运算符< / 的系统沙罗> ;

Proc glm运算符上面glm流程的开始,此句后的的系统可有表列计划。

DATA=多达据集名

ALPHA=p

MANOVA

MULTIPASS

NAMELEN=n

NOPRINT

ORDER=DATA | FORMATTED | FREQ | INTERNAL

OUTSTAT=SAS-data-set

可以见到此的系统沙罗极为少比anova流程的多了一个计划,即“alpha=p”的系统,此的系统用来原则上量度流程之前所选用的显著性素质。其余的系统的词汇和表达方式与anova流程的不同,这里才会行多说了。

大家可以较为一下glm流程和anova流程所限于的运算符,anova流程之前限于的所有运算符都包涵在glm流程所限于的运算符之前,其词汇和机制也都是只不过不同的,这一部分运算符这里就才会行赘述。

Contrast运算符使你可以用自定义的方式顺利进行假设验,它需要注意到在model运算符之后,如果加到manova运算符、repeated运算符、random运算符或test运算符,contrast运算符需要注意到在这些运算符之前。上面用来标签所顺利进行的验,意在标签的文字或符号必需用单引号括上去。效理应运算符意在原则上假设验的考量(两一组),这些考量(两一组)需要是model运算符之前注意到过的。效理应运算符后的常多达乘积意在原则上相理应考量(两一组)各素质的个数,在原则上各素质的情况下顺利进行相关考量的系统性。

Estimate运算符可付诸对线性方程的估计,它也需要注意到在model运算符之后,运用于的规章和contrast运算符只不过不同。其之前的运算符金属元素的表达方式和词汇也与contrast运算符不同。

Lsmeans运算符意在指示SAS对原则上的考量(两一组)量度理应给定的很小二乘均多达并控制器到结果之前。

Output运算符我们在过去的具体内容之前接触过,其机制和词汇和过去的具体内容也只不过不同,各位再行三参考过去的具体内容运用于。

Random运算符意在原则上哪些考量(两一组)是随机给定,即近似于这些考量(两一组),采样是来自于正态基本上的随机采样,这样SAS可对相理应考量(两一组)作随机效理应假设的系统性。

对以上几条运算符的理解只不过来自SAS的帮助文件,由于专业知识最深处不够以及英语素质极小的缘故,以上关于这几种运算符的具体内容亦会依赖于很多的出错与偏差,还再行三各位时也体谅,也期望各位能时也指点。全盘的具体内容我亦会给出及时的预览。

二、有所不同新设计类型资讯的方差系统性

必需顺利进行各种类型均多达较为的资讯新设计类型有只不过随机新设计、随机计量一组新设计、析因新设计、单调量度新设计等,下面我们分别以程序中来科学实验SAS之前这类多达据的系统性方法。

1. 只不过随机新设计资讯的方差系统性(含均多达间的多重较为)

例4-1 为知晓烫伤后有所不同时期切痂对胰脏三甲基腺苷(简写为ATP)含量的严重影响,将30只雄性大鼠随机分3一组, 每一组10只:A一组为烫伤大鼠,B一组为烫伤后24天内(休克期)切痂一组,C一组为烫伤后96天内(非休克期)切痂一组,全部鸟类统一在烫伤后168天内斩首并量度其胰脏的ATP含量,结果见表4.1。再行三验三一组大鼠胰脏的ATP含量所谓不同点。

表4.1 大鼠烫伤后胰脏ATP含量的量度结果(计量:mg)

A一组

B一组

C一组

7.76

11.14

10.85

7.71

11.60

8.58

8.43

11.42

7.19

8.47

13.85

9.36

10.30

13.53

9.59

6.67

14.16

8.81

11.73

6.94

8.22

5.78

13.01

9.95

6.61

14.18

11.26

6.97

17.72

8.68

针对以上资讯,我们运用于anova流程来付诸方差系统性的机制。编制SAS程序在如下:

data temp;

infile 'e:datadata4_1.txt';

do i=1 to 10;

do g=1 to 3;

input x@@;

output;

end;

end;

proc anova;

class g;

model x=g;

Run;

quit;

本段程序在之前加到SAS的尿素运算符,其播放器和时也达量度机程序员理技术的发展软件之前的不同,我近似于它的词汇和表达方式不用再行多罗嗦了。变化多端尿素运算符在SAS程序员流程之前可以节省很多的精力。此处因为多达据位于有所不同的列之前,每行包涵每两一组的一个多达据,用尿素运算符正好可以产生一个都是有所不同分一组的给定,省去了一齐输入带来的麻烦。无必需提醒各位要特别强调的是,尿素运算符的每一个do运算符,需要有一个相理应的end运算符来告一段落它,否则就亦会出错。

此处无必需系统性的给定(理应给定)是x,干预考量(自给定)为分一组给定g,因而model运算符之前效理应运算符为x=g。

呈交上述程序在,结果如下。

The SAS System 10:49 Monday, October 14, 2002 5

The ANOVA Procedure

Class Level Information

Class Levels Values

g 3 1 2 3

Number of observations 30

The SAS System 10:49 Monday, October 14, 2002

分享:
365整形网 整形医院哪家好 五官整容整形 整形医院咨询 整形知识 整形医生 美容整形 整形医院排名 整形医院咨询 整形专业知识 济南整形医院 C店求购 快手买卖平台 抖音号出售 出售公众号网站 买卖小红书网站 出售快手网 出售小红书网站