用R进行单因素方差分析

脚本与注释

datafilename="http://personality-project.org/r/datasets/R.appendix1.data"  # 数据文件所在位置
mydata=read.table(datafilename,header=T)   #读取数据文件
names(mydata) <-c("group", "Alertness")    # 给变量命名，一个命名为group，另一个为Alertness
aovresult <- aov(Alertness~group,data=mydata)   # 进行单因素方差分析，结果保存到avoresult
summary(aovresult)                              # 显示单因素方程分析总体结果
print(model.tables(aovresult,"means"),digits=3) # 显示平均值、样本数量等
boxplot(Alertness~group,data=mydata)            # 画图