#
# 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.
#
#-----

source('utils.R')

# Ex. 11.38 (EPage 440):
#
data = c( 200.0, 204.2, 215.5, 219.5, 272.7, 276.9, 299.5, 302.7, 166.6, 172.6, 186.4, 192.0, 232.6, 240.8, 253.4, 261.6 )
D = matrix( data, nrow=2^3, ncol=2, byrow=T )

res = YATES_METHOD( D )


# Ex. 11.39
#
D = matrix( c( 106, 93, 116, 198, 200, 214, 197, 202, 185, 329, 331, 307, 149, 169, 135, 243, 247, 220, 255, 230, 252, 383, 360, 364 ), nrow=2^3, ncol=3, byrow=T )

res = YATES_METHOD( D )


# Ex. 11.40:
#
D = matrix( c( 9.72, 13.50, 9.80, 14.04, 10.13, 11.27, 11.80, 11.30, 12.70, 11.37, 11.96, 12.05, 
               11.38, 9.92, 11.80, 11.10, 13.15, 13.00, 10.60, 12.37, 10.37, 12.00, 11.30, 11.64, 
               13.05, 14.55, 11.15, 15.00, 12.70, 14.10, 13.20, 16.12 ), 
             nrow=2^4, ncol=2, byrow=T ) 

res = YATES_METHOD( D )


# Ex. 11.41:
#
D = matrix( c( 144, 154, 158, 239, 227, 244, 232, 242, 246, 364, 362, 346, 194, 162, 203, 284, 295, 291, 291, 287, 297, 411, 406, 395 ), nrow=2^3, ncol=3, byrow=T ) 

res = YATES_METHOD( D )


# Ex. 11.42:
#
D = matrix( c( 866, 862, 800, 946, 800, 840, 774, 834, 746, 709, 789, 646, 1017, 990, 954, 1028, 906, 977, 817, 783, 771, 829, 806, 691, 988, 808, 650, 966, 976, 876, 702, 658, 650, 784, 700, 596, 922, 808, 868, 1056, 870, 908, 798, 726, 700, 752, 714, 714 ), 
            nrow=2^4, ncol=3, byrow=T ) 

res = YATES_METHOD( D )


# Ex. 11.44:
#
D = matrix( c( 91, 85, 92, 94, 86, 83, 85, 90, 72, 78, 68, 79, 69, 75, 72, 71 ), nrow=2^4, ncol=1 )


# Ex. 11.45:
#
D = matrix( c( 618, 598,
	       583, 560,
	       477, 525,
	       421, 462,
	       601, 595,
	       550, 589,
	       505, 484,
	       452, 451,
	       598, 585,
	       587, 541,
	       480, 508,
	       462, 449,
	       603, 577,
	       571, 552,
	       502, 508,
	       449, 455 ), nrow=2^4, ncol=2, byrow=T )

res = YATES_METHOD( D )


# Ex. 11.48:
#
DF = read.csv( "../../Data/CH11/ex11-48.txt", header=TRUE, quote="'" )


# Ex. 11.49:
#
D = matrix( data = c( 70.4, 72.1, 70.4, 67.4, 68.0, 73.8, 67.0, 67.8, 66.6, 67.5, 64.0, 66.8, 70.3, 67.9, 65.9, 68.0 ), nrow=2^4, ncol=1 )