无原始数据情况下如何进行方差分析

晚风

在没有原始数据的情况下，进行方差分析有两条途径：

1、  需提供均值、SD和样本数三个统计量。如果能提供完整的组内和组间统计信息，可以进行单因素和多因素方差分析，但这种方法不能进行正态性检验和方差齐性检验。

2、  利用提供的均值、SD和样本数，生成符合N(mean,sd)分布的原始数据，可以进行单因素或多因素方差分析。

以下仅提供第一种方法的R语言代码。

mydata <- read.table(textConnection(
“year   n  FVC    s
wenxian 20 9225.3 433.2701
lankao  20 9184.71 334.2981″
), header=TRUE,row.names=NULL)

mydata$year <- factor(mydata$year, levels=mydata$year)
row.names(mydata) <- mydata$year
library(HH)
mydata.aov <- aov.sufficient(FVC ~ year,data=mydata)
summary(mydata.aov)