wbw <- read.table("http://www.stat.ufl.edu/~winner/data/whole_breast_weight.dat", header=F,col.names=c("trt","repnum","base","meth","wtg")) attach(wbw); names(wbw) # Treatment Ordering: BS, BSM, BC, BCM trt.f <- factor(trt) oneway.test(wtg ~ trt.f, var.equal=T) oneway.test(wtg ~ trt.f, var.equal=F) bartlett.test(wtg ~ trt.f)