# 
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

library(MASS)

fm = formula( Rabbit )

cbind( 
Treatment = unstack( Rabbit, Treatment ~ Animal )[,1], 
Dose = unstack( Rabbit, Dose ~ Animal )[,1], 
unstack( Rabbit, BPchange ~ Animal )
)

# Another way of doing this: 
require(reshape)
recast( Rabbit, Treatment + Dose ~ Animal, measure.var="BPchange" )