例題：

1.在某工廠生產的型號的圓軸中任取20個，測得其直徑數據如下：

15.63 15.13 15.46 15.40 15.56 15.35 15.56 15.38 15.21
15.58 15.57 15.36 15.48 15.46 15.52 15.29 15.42 15.69
求上序數據的樣本均值，中位數，四分位數，樣本方差，極差，變異系數，二階、三階和四階中心距；求偏度，峰度，并把數據中心距和標準差。

2.下列列出了84個伊特拉斯人男子頭顱的最大寬度（單位：mm）

148 132 138 154 142 150 146 155 158 150 140 147 148 144

150 149 145 149 158 143 141 144 144 126 140 144 142 141 140

145 135 147 146 141 140 136 146 142 137 148 154 137 139 143

140 131 141 143 149 148 135 148 152 143 144 141 143 147 146

150 132 142 142 143 153 149 146 149 138 142 149 142 137 134

144 146 147 140 142 137 152 140 145

在顯著性水平a=0.01下，試檢驗上述頭顱的最大寬度數據是否來自正態總體？

實驗代碼：
data YZTJ;
input d @@;
cards;
15.28 15.63 15.13 15.46 15.40 15.56 15.35 15.56 15.38 15.21
15.48 15.58 15.57 15.36 15.48 15.46 15.52 15.29 15.42 15.69
;
run;
proc means noprint data=YZTJ mean median qrange var range cv skewness kurtosis std;
var d;
output out=YZTJResult1 mean=mean median=median qrange=qrange var=var range=range cv=cv skewness=skewness kurtosis=kurtosis std=std
run;
data YZTJResult2;
input x1-x20 @;
mean=sum(of x1-x20)/20;
B2=sum((x1-mean)**2,(x2-mean)**2,(x3-mean)**2,(x4-mean)**2,
(x5-mean)**2,(x6-mean)**2,(x7-mean)**2,(x8-mean)**2,(x9-mean)**2,(x10-mean)**2,(x11-mean)**2,(x12-mean)**2,
(x13-mean)**2,(x14-mean)**2,(x15-mean)**2,(x16-mean)**2,(x17-mean)**2,(x18-mean)**2,(x19-mean)**2,(x20-mean)**2)/20;
B3=sum((x1-mean)**3,(x2-mean)**3,(x3-mean)**3,(x4-mean)**3,
(x5-mean)**3,(x6-mean)**3,(x7-mean)**3,(x8-mean)**3,(x9-mean)**3,(x10-mean)**3,(x11-mean)**3,(x12-mean)**3,
(x13-mean)**3,(x14-mean)**3,(x15-mean)**3,(x16-mean)**3,(x17-mean)**3,(x18-mean)**3,(x19-mean)**3,(x20-mean)**3)/20;
B4=sum((x1-mean)**4,(x2-mean)**4,(x3-mean)**4,(x4-mean)**4,
(x5-mean)**4,(x6-mean)**4,(x7-mean)**4,(x8-mean)**4,(x9-mean)**4,(x10-mean)**4,(x11-mean)**4,(x12-mean)**4,
(x13-mean)**4,(x14-mean)**4,(x15-mean)**4,(x16-mean)**4,(x17-mean)**4,(x18-mean)**4,(x19-mean)**4,(x20-mean)**4)/20;
run;
cards;
15.28 15.63 15.13 15.46 15.40 15.56 15.35 15.56 15.38 15.21 15.48 15.58 15.57 15.36 15.48 15.46 15.52 15.29 15.42 15.69
;
run;
data YZTJResult;
merge study.YZTJResult1 study.YZTJResult2;
drop x1-x20;
run;
proc print data=study.YZTJResult;
run;
實現結果：

2、
data TLTJ;
input d @@;
cards;
141 148 132 138 154 142 150 146 155 158 150 140 147 148 144
150 149 145 149 158 143 141 144 144 126 140 144 142 141 140
145 135 147 146 141 140 136 146 142 137 148 154 137 139 143
140 131 141 143 149 148 135 148 152 143 144 141 143 147 146
150 132 142 142 143 153 149 146 149 138 142 149 142 137 134
144 146 147 140 142 137 152 140 145
;
run;
proc univariate data=TLTJ normal;
var d;
run;
實現結果：

分析實驗結果：
（1）代碼1中，先建立study.YZTJ數據集把圓軸直徑數據導入，用proc means方法把計算的統計關鍵詞output到study.YZTJResult1中，接著創建臨時數據集temp1和temp2，分別把mean和x1-x20的數據導入其中，創建study.YZTJResult2并接temp1和temp2，算出2到4階的中心矩B2-B4，最后把study.YZTJResult1和study.YZTJResult2并接，得到最終結果數據集study.YZTJResult。

（2）創建數據集study.TLTJ，用proc univariate方法和normal選項進行正態性檢驗，由結果可以看出，在顯著性水平a=0.01下，數據量小于2000，算出的p值0.7>0.01，故接受H0，拒絕H1，即頭顱的最大寬度數據是來自正態總體。

總結

以上是生活随笔為你收集整理的SAS用proc means和proc univariate求数据的样本均值、中位数、四分位数、样本方差、极差、变异系数、二阶、三阶和四阶中心矩、偏度、峰度、标准差和对数据进行正态性检验的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。