##    Create data
set.seed(1)
m <- matrix(rep(1,3*24),ncol=3)  
m[9:16,3] <- 3 ; m[17:24,] <- 3    #create 3 groups
m <- m+rnorm(24*3,0,0.5)           #add noise
m <- floor(10*m)/10                #just one digits
r2xcluster(m)
## And once you have Xcluster program:
## Not run: 
#   system('Xcluster -f xcluster.txt -e 0 -p 0 -s 0 -l 0')
#   h <- xcluster2r('xcluster.gtr')
#   plot(h,hang=-1)
#   ## End(Not run)
Run the code above in your browser using DataLab