多元統(tǒng)計實驗三方差分析(共5頁)

1、桂林電子科技大學(dxu) 數(shù)學與計算(j sun)科學學院實驗報告院系數(shù)學與計算科學學院學號姓名成績課程名稱應用多元統(tǒng)計分析實驗實驗項目名 稱實驗三 方差分析一 ，實驗目的用SAS軟件進行單個總體均值的推斷二，實驗原理1means過程 功能：對數(shù)值型變量給出簡單描述性統(tǒng)計值。 (16種統(tǒng)計量，P51,默認的只有5種) 也可以檢驗正態(tài)總體均值是否為零。 格式： proc means probt; /* probt 用于檢驗正態(tài)總體均值是否為零*/ var x ; run; 2. Univariate過程 功能：UNIVARIATE過程除了可以完成與MEANS過程相同的基本統(tǒng)計量外，還可以計算變

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

3、16個哥特式教堂的樣本數(shù)據(jù)如下： 序號長度/m中殿高度/m序號長度/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 試檢驗, .四，實驗過程原始記錄（數(shù)據(jù)，圖表，計算等）data shiyan3;input no x1 x2 ;long=x1-1

4、45.29;haigh=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 用于檢驗正態(tài)總體均值是否為零*/ var long haigh ; run;

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

6、。2. 高度圖2.2.1 正態(tài)性檢驗圖2.2.2 位置檢驗由于樣本量0.05,接受原假設，即羅馬式教堂中殿的平均高度服從正態(tài)分布。并由Students t檢驗的p值為0.91570.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

7、.5431.397112.7817.3715186.2331.398154.2824.9916144.1725.6; proc iml data=shiyan3;n=16; p=2; /*定義樣本容量和總體的維數(shù)*/u0=145.29 22.69; /*miu 0 均值向量*/use shiyan3; /*打開SAS數(shù)據(jù)集*/xn=x1 x2; /*變量名矩陣（行向量）*/read all var xn into x; /*讀取xa矩陣中所有的變量，將其值存放到X矩陣（即數(shù)據(jù)矩陣）中*/print x;in=16 1; /*每個元素都是 1 的 1行20列矩陣（即行向量）*/print in;

8、xb=(in*x)/n; /*計算三個分量的均值*/print xb;xc=xb-u0;print xc;a=x*x-n*xb*xb; /*計算樣本離差矩陣*/S=a/(n-1); /*計算樣本協(xié)方差矩陣*/print a s;is=inv(s); /*計算樣本協(xié)方差矩陣的逆*/print is;t2=n*xc*is*xc; /*計算霍特林統(tǒng)計量值*/f=(n-p)*t2/(n-1)*p); /*計算F值*/print t2 f;p0=1-probf(f,p,n-p); /*計算P值*/fa=finv(0.95,p,n-p); /*計算水平為0.05的F臨界值*/print p0 fa ;run;qui