SPSS相关编程：

data list free/x1 x2 x3 s1 s2 s3 .
begin data
17 13 16 7.4 8.8 3.2
end data.
compute ss1=s1**2*49+s2**2*57+s3**2*51.
compute MS1=ss1/157.
compute X=(x1*50+x2*58+x3*52)/160.
compute ss2=50*(x1-X)**2+58*(x2-X)**2+52*(x3-X)**2.
compute MS2=ss2/2.

compute F=MS2/MS1.
compute p=1-CDF.F(F,2,157).
execute.

举例：

原始数据的统计描述：
组别　例数　平均数　　　标准差　　　标准误
1　　　50 　 17.0000　  7.4000 　　 1.0465
2　　　58　　13.0000 　 8.8000　　　1.1555
3 　　 52　　16.0000　  3.2000　　　0.4438

参数统计应用条件检查：
1) 正态性检验(矩法):
样本含量足够大,不必进行正态性检验!

2) 方差齐性检验:
卡方值=46.4058, P=0.0000
按α=0.0500水准，可认为该资料方差不齐!

参数统计结果：

方差分析:
Ho:各个总体均数相等
H1:各个总体均数不相等或不全相等
α=0.0500

方差分析结果
=========================================
变异来源　　SS　　　 ν 　　　MS　　　　　　F 　　　P
------------------------------------------------------------
总 　　　　8095.4600　159
组间　　　475.9000　　　2　　237.9500 　　4.90　0.0086
组内　　　7619.5600　　157　　48.5322
============================================================

结论:经过方差分析,得P=0.0086,按α=0.0500水准拒绝Ho,接受H1,故可认为各组总体均数不相等!

3个样本均数两两比较的q检验(Newman-Keuls法)
=============================================
组别　　　　两均数之差　　　组数　　Q值　　　   P值
------------------------------------------------------------------
第 1与第 2　4.0000　　　　　3　　　 4.2077　　　<0.05
第 1与第 3　1.0000 　　　　 2　　 　1.0249　　　>0.05
第 2与第 3　3.0000　　　　　2 　　　3.1889 　　 <0.05
=============================================

平均住院日（天）

原始数据的统计描述：
组别　例数　　平均数 　　标准差　　　标准误
1　　50 　　　28.0000 　　6.1000　　 0.8627
2　　58　　　 22.0000　　 5.7000　　 0.7484
3 　 52 　　  25.0000 　　1.7000 　　0.2357

参数统计应用条件检查：
1) 正态性检验(矩法):
样本含量足够大,不必进行正态性检验!

2) 方差齐性检验:
卡方值=70.8595, P=0.0000
按α=0.0500水准，可认为该资料方差不齐!

参数统计结果：

方差分析:
Ho:各个总体均数相等
H1:各个总体均数不相等或不全相等
α=0.0500

方差分析结果
=============================================
变异来源　　　　SS　　   ν 　　 MS　　　　　F 　　　P
------------------------------------------------------------
总 　　　　　4791.0100　 159
组间　　　　 968.4000　　2　　 484.2000　  19.89　0.0000
组内　　　　 3822.6100　 157 　24.3478
============================================================

结论:经过方差分析,得P=0.0000,按α=0.0500水准拒绝Ho,接受H1,故可认为各组总体均数不相等!

3个样本均数两两比较的q检验(Newman-Keuls法)
=======================================
组　别 　　　两均数之差　　组数　  Q值　　　P值
------------------------------------------------------------------
第 1与第 2　　6.0000 　　　 3　　 8.9109 　 <0.01
第 1与第 3 　 3.0000 　　　 2　　 4.3410　　<0.01
第 2与第 3  　3.0000　　　  2　　 4.5022　　<0.01
========================================