多元統(tǒng)計(jì)實(shí)驗(yàn)三方差分析

1、桂林電子科技大學(xué) 數(shù)學(xué)與計(jì)算科學(xué)學(xué)院實(shí)驗(yàn)報(bào)告院系數(shù)學(xué)與計(jì)算科學(xué)學(xué)院學(xué)號(hào)姓名成績(jī)課程名稱(chēng)應(yīng)用多元統(tǒng)計(jì)分析實(shí)驗(yàn)實(shí)驗(yàn)項(xiàng)目名 稱(chēng)實(shí)驗(yàn)三 方差分析一 ，實(shí)驗(yàn)?zāi)康挠肧AS軟件進(jìn)行單個(gè)總體均值的推斷二，實(shí)驗(yàn)原理1means過(guò)程 功能：對(duì)數(shù)值型變量給出簡(jiǎn)單描述性統(tǒng)計(jì)值。 (16種統(tǒng)計(jì)量，P51,默認(rèn)的只有5種) 也可以檢驗(yàn)正態(tài)總體均值是否為零。 格式： proc means probt; /* probt 用于檢驗(yàn)正態(tài)總體均值是否為零*/ var x ; run; 2. Univariate過(guò)程 功能：UNIVARIATE過(guò)程除了可以完成與MEANS過(guò)程相同的基本統(tǒng)計(jì)量外，還可以計(jì)算變量的極端值、分位數(shù)，生成

2、頻率表，并支持對(duì)數(shù)據(jù)進(jìn)行正態(tài)性檢驗(yàn)。格式： proc univariate normal ; /* normal 用于檢驗(yàn)總體是否服從正態(tài)*/ var x ; run;3. IML 過(guò)程由于在DATA數(shù)據(jù)步的過(guò)程中， SAS系統(tǒng)并沒(méi)有提供向量或矩陣的輸入，要想建立一個(gè)矩陣的數(shù)據(jù)，就要用到SAS/IML模塊，它有著一套很強(qiáng)大、靈活的運(yùn)算語(yǔ)言。在IML的環(huán)境下不但可以進(jìn)行矩陣的各種運(yùn)算，而且還可以與DATA步一起使用實(shí)現(xiàn)各種復(fù)雜的計(jì)算和操作功能三，實(shí)驗(yàn)內(nèi)容中世紀(jì)英國(guó)的建筑風(fēng)格，早期是羅馬式，后來(lái)有哥特式。羅馬式教堂平均長(zhǎng)度為145.29m，中殿的平均高度為22.69m.n=16個(gè)哥特式教堂的樣本數(shù)

3、據(jù)如下： 序號(hào)長(zhǎng)度/m中殿高度/m序號(hào)長(zhǎng)度/m中殿高度/m1158.1930.499124.0521.95268.5822.861089.9226.82391.4415.851183.2116.764127.4118.9012126.4920.425124.6620.731355.4713.726129.5426.2114161.5431.397112.7817.3715186.2331.398154.2824.9916144.1725.6 試檢驗(yàn), .四，實(shí)驗(yàn)過(guò)程原始記錄（數(shù)據(jù)，圖表，計(jì)算等）data shiyan3;input no x1 x2 ;long=x1-145.29;haigh=

4、x2-22.69;cards;1158.1930.499124.0521.95268.5822.861089.9226.82391.4415.851183.2116.764127.4118.9012126.4920.425124.6620.731355.4713.726129.5426.2114161.5431.397112.7817.3715186.2331.398154.2824.9916144.1725.6; proc means data=shiyan3 probt; /* probt 用于檢驗(yàn)正態(tài)總體均值是否為零*/ var long haigh ; run;proc univari

5、ate data=shiyan3 normal ;var long haigh;run;方法1. 從上表可以看出長(zhǎng)度的檢驗(yàn)結(jié)果p值為0.017，小于0.05.因此，拒絕原假設(shè)，即認(rèn)為羅馬式教堂平均長(zhǎng)度不是145.29m。高度的檢驗(yàn)結(jié)果的p值為0.3157，大于0.05，因此，接受原假設(shè)，即認(rèn)為中殿的平均高度為22.69m。方法2：1.長(zhǎng)度圖2.1.1圖2.1.2由于樣本量0.05,接受原假設(shè)，即羅馬式教堂平均長(zhǎng)度服從正態(tài)分布。并由Students t檢驗(yàn)的p值為0.0170.05，所以y的均值不為0,即即羅馬式教堂平均長(zhǎng)度不是145.29m。2. 高度圖2.2.1 正態(tài)性檢驗(yàn)圖2.2.2 位置

6、檢驗(yàn)由于樣本量0.05，所以y的均值為0,即即羅馬式教堂中殿的平均高度為22.69m。方法3：data shiyan3;input no x1 x2 ;cards;1158.1930.499124.0521.95268.5822.861089.9226.82391.4415.851183.2116.764127.4118.9012126.4920.425124.6620.731355.4713.726129.5426.2114161.5431.397112.7817.3715186.2331.398154.2824.9916144.1725.6; proc iml data=shiyan3;n

7、=16; p=2; /*定義樣本容量和總體的維數(shù)*/u0=145.29 22.69; /*miu 0 均值向量*/use shiyan3; /*打開(kāi)SAS數(shù)據(jù)集*/xn=x1 x2; /*變量名矩陣（行向量）*/read all var xn into x; /*讀取xa矩陣中所有的變量，將其值存放到X矩陣（即數(shù)據(jù)矩陣）中*/print x;in=16 1; /*每個(gè)元素都是 1 的 1行20列矩陣（即行向量）*/print in;xb=(in*x)/n; /*計(jì)算三個(gè)分量的均值*/print xb;xc=xb-u0;print xc;a=x*x-n*xb*xb; /*計(jì)算樣本離差矩陣*/S=a/(n-1); /*計(jì)算樣本協(xié)方差矩陣*/print a s;is=inv(s); /*計(jì)算樣本協(xié)方差矩陣的逆*/print is;t2=n*xc*is*xc; /*計(jì)算霍特林統(tǒng)計(jì)量值*/f=(n-p)*t2/(n-1)*p); /*計(jì)算F值*/print t2 f;p0=1-probf(f,p,n-p); /*計(jì)算P值*/fa=finv(0.95,p,n-p); /*計(jì)算水平為0.05的F臨