sscp 5 5 5 format labels ) Morrison data p. 225 ) Five 5 by 5 MANOVA hypothesis and error matrices are ) combined into a 3-dimensional array sscp. ) ) You extract matrix 4, say, by h4 <- matrix(sscp[4,,]) To ) extract the part of matrix 4 having to do only with the ) responses and not U, use h4 <- matrix(sscp[4,-1,-1]) ) ) There are 4 hypothesis matrices for a 2x2x2 completely ) randomized factorial experiment with 1 covariate U and 4 ) responses Y1 - Y4 in that order. ) ) Var. 1: U = initial wt of rats ) Var. 2: Y1 = weight loss 1 day after irradiation ) Var. 3: Y2 = weight loss 3 day after irradiation ) Var. 4: Y3 = weight loss 6 day after irradiation ) Var. 5: Y4 = weight loss 7 day after irradiation ) ) Factor A: Radiation level -- 500 r. or 600 r. ) Factor B: Age -- young or adult ) Factor C: Sex -- female or male ) ) There were 3 rats for each combination of factors ) ) The hypothesis matrices in the order they appear are ) HA = sscp[1,,] ) HB = sscp[2,,] ) HC = sscp[3,,] ) H4 = sscp[4,,] = HAB + HAC + HBC + HABC (4 df) ) The error matrix is ) E = sscp[5,,] ) )"6x%lf %lf %lf %lf %lf" (6X,5F8.2) HA 70.04 10.25 59.79 -52.96 -71.75 HA 10.25 1.50 8.75 -7.75 -10.50 HA 59.79 8.75 51.04 -45.21 -61.25 HA -52.96 -7.75 -45.21 40.04 54.25 HA -71.75 -10.50 -61.25 54.25 73.50 HB 11926.04 847.08 66.88 1493.54 1605.00 HB 847.08 60.17 4.75 106.08 114.00 HB 66.88 4.75 .38 8.38 9.00 HB 1493.54 106.08 8.38 187.04 201.00 HB 1605.00 114.00 9.00 201.00 216.00 HC 51987.04 1210.08 4421.46 6469.29 3816.42 HC 1210.08 28.17 102.92 150.58 88.83 HC 4421.46 102.92 376.04 550.21 324.58 HC 6469.29 150.58 550.21 805.04 474.92 HC 3816.42 88.83 324.58 474.92 280.17 H4 3669.50 -204.66 124.84 32.34 172.75 H4 -204.66 86.01 -52.33 -54.33 -61.16 H4 124.84 -52.33 161.16 98.16 71.76 H4 32.34 -54.33 98.16 128.16 89.76 H4 172.75 -61.16 71.76 89.76 73.51 E 1326.00 36.00 284.65 66.67 -90.67 E 36.00 95.98 59.33 53.34 49.33 E 284.65 59.33 411.34 117.67 57.66 E 66.67 53.34 117.67 342.68 276.32 E -90.67 49.33 57.66 276.32 385.32 sscp$labels 15 quoted COLUMNS ) Labels for matrices )"%s %s %s %s %s" "HA" "HB" "HC" "H4" "E" "InitialWt" "Day 1" "Day 3" "Day 6" "Day 7" "InitialWt" "Day 1" "Day 3" "Day 6" "Day 7"