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.
``{r} #R 程序
n <-c(201,285,252) #例数-样本数
means <-c(85.26,88.29,90.66) #各组均数
s <-c(6.87,5.32,9.02) #各组标准差
anova<-function(n,means,s){
ss1 <-s^2*(n-1)
ss1 <- sum(ss1) #组内平方和
ms1 <- ss1/(sum(n)-length(n)) #组内均方
X <- sum(means*n)/sum(n) #总均数
ss2 <- sum(n*(means-X)^2) #组间平方和
ms2 <- ss2/(length(n)-1) #组间均方
ms2
F <-ms2/ms1
P <-1- pf(F,length(n)-1,sum(n),lower.tail=F)
results<-list("组间均方：",ms1,"组内均方：",ms2,"F-value:",F,"P-value:",P)
return(results)
}
anova(n,means,s)举例：
原始数据的统计描述：
组别 例数 平均数 标准差 标准误
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
========================================
参考：http://blog.sina.com.cn/s/blog_4e9e6d470100b5ow.html
http://www.360doc.com/content/19/0124/20/52645714_811083449.shtml