pdf("ex0604.pdf")

ex0604dat <- read.fwf("http://www.stat.ufl.edu/~winner/data/biostat/ex0604.dat", width=c(8,8,8),
col.names=c("subj", "formu", "tmax"))

# Create qualitative factor variable for formu, and assign names to levels

fformu <- factor(ex0604dat$formu, levels=1:3)
levels(fformu) <- c("cap_f", "cap_nf", "entct_f")

# We have to assign subj (Subject id) as a factor level or the linear model will treat
#   it as a numeric (continuous) variable and fit a regression

ex0604 <- data.frame(tmax=ex0604dat$tmax, fformu, subj=factor(ex0604dat$subj))

attach(ex0604)


# Run Friedman's test for the RBD with subject as block and formulation as treatment

friedman.test(tmax ~ fformu|subj)



dev.off()