#6.2 Wine Example
x<-rbind(
c(2,3.5,3.5,2,2),
c(1.5,4,1.5,4,4),
c(2,3.5,2,3.5,2),
c(2,3.5,3.5,2,2),
c(4,4,4,1.5,1.5),
c(2,3.5,3.5,2,2),
c(2,3.5,3.5,2,2),
c(2.5,5,2.5,2.5,2.5),
c(2,3.5,3.5,2,2),
c(1.5,4,4,4,1.5))
colnames(x)<-c("A","B","C","D","E")

#Have to get ranks from x*10 as xtotab assumes integer
xtotab(x*10)
U(x*10)
BlockAnalysis(x*10)

#Note that 4 degrees of freedom have been lost - the total degrees of freedom
#is usually (number of products - 1)^2, which is (5-1)^2=16 here

#This is because U(x*10) is not of full rank (see the discussion at the end of chapter 3.6)
eigen(U(x*10))$values
det(U(x*10))