仅有各组样本例数,均数、标准差如何进行方差分析?提供R语言代码

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