OPTIONS NOCENTER LINESIZE=81 PAGESIZE=66 FORMCHAR = "|----|+|---+="; FILENAME IN1 'C:\DATA\PROZAC\WTDELTbM.DAT'; TITLE1 'Prozac Weight Data - data from visits 2 through 5'; DATA A; INFILE IN1; INPUT PATIENT 5. (WTKGCH2-WTKGCH5)(4*8.3) compdata 3. COND 3. SEXM 3. BMIV1 8.3 (STAT2-STAT5)(4*3.); PROC FORMAT; VALUE COND 1 = 'Placebo' 2 = 'LowDose' 3 = 'HighDose'; VALUE SEXM 0 = 'Female' 1 = 'Male'; VALUE STAT 0 = 'Abstainer' 1 = 'Chipper' 2 = 'Smoker'; VALUE COMPDATA 0 = 'No' 1 = 'Yes'; PROC MEANS; RUN; /* means across time by condition for all subjects */ PROC SORT; BY COND ; PROC MEANS; VAR WTKGCH2-WTKGCH5; BY COND ; FORMAT COND COND. ; RUN; /* means across time by condition for completers */ DATA C; SET A; IF COMPDATA=1; PROC SORT; BY COND ; PROC MEANS; VAR WTKGCH2-WTKGCH5; BY COND ; FORMAT COND COND. ; RUN; /* set up the data for plotting the group means across time */; /* all data */; PROC SUMMARY NWAY DATA=A; CLASS COND; VAR WTKGCH2-WTKGCH5; OUTPUT OUT=B MEAN = MWTCH2-MWTCH5; DATA PLOT; SET B; ARRAY MWTCH(4) MWTCH2-MWTCH5; DO TIME = 1 TO 4; MEANWTCH = MWTCH(TIME); OUTPUT; END; DROP MWTCH2-MWTCH5; PROC PLOT; PLOT MEANWTCH*TIME=COND; FORMAT COND COND.; RUN; /* set up the data for plotting the group means across time */; /* subjects with complete data across time */; PROC SUMMARY NWAY DATA=C; CLASS COND; VAR WTKGCH2-WTKGCH5; OUTPUT OUT=D MEAN = MWTCH2-MWTCH5; DATA PLOT; SET D; ARRAY MWTCH(4) MWTCH2-MWTCH5; DO TIME = 1 TO 4; MEANWTCH = MWTCH(TIME); OUTPUT; END; DROP MWTCH2-MWTCH5; PROC PLOT; PLOT MEANWTCH*TIME=COND; FORMAT COND COND.; RUN; /* perform multivariate repeated measures analysis */ TITLE2 'multivariate repeated measures analyses'; PROC GLM DATA=A; CLASS COND; MODEL WTKGCH2-WTKGCH5 = COND / NOUNI; CONTRAST 'low vs placebo' COND -1 1 0 ; CONTRAST 'high vs placebo' COND -1 0 1 ; REPEATED TIME POLYNOMIAL / SUMMARY PRINTM PRINTE; LSMEANS COND; RUN; Prozac Weight Data - data from visits 2 through 5 The MEANS Procedure Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- PATIENT 989 797.0687563 462.0553200 1.0000000 1583.00 WTKGCH2 894 0.2972472 1.3872107 -4.5450000 5.9090000 WTKGCH3 852 0.4676068 1.5940949 -5.4550000 8.6360000 WTKGCH4 785 0.6933529 1.8282857 -7.0450000 7.7270000 WTKGCH5 724 1.0754144 2.1833625 -7.0450000 10.9090000 compdata 989 0.7229525 0.4477666 0 1.0000000 COND 989 1.9949444 0.8179259 1.0000000 3.0000000 SEXM 989 0.3953488 0.4891729 0 1.0000000 BMIV1 984 25.3690274 4.4910708 16.6370000 48.5510000 STAT2 898 0.5645880 0.7634005 0 2.0000000 STAT3 858 0.6247086 0.7677546 0 2.0000000 STAT4 790 0.7936709 0.8034845 0 2.0000000 STAT5 727 0.8363136 0.8344521 0 2.0000000 ------------------------------------------------------------------------------- COND=Placebo The MEANS Procedure Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 310 0.7738290 1.2972127 -3.8640000 4.5450000 WTKGCH3 285 1.2215263 1.4740143 -2.5000000 8.6360000 WTKGCH4 263 1.6471445 1.6507896 -3.8640000 7.7270000 WTKGCH5 244 1.9340410 2.0173265 -6.0910000 8.8640000 ------------------------------------------------------------------------------- COND=LowDose Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 296 0.1353176 1.2591884 -4.0910000 5.9090000 WTKGCH3 286 0.1848986 1.4311115 -4.3180000 6.5000000 WTKGCH4 263 0.3457871 1.6566482 -7.0450000 5.2730000 WTKGCH5 238 0.8194790 1.9473696 -7.0450000 5.0000000 ------------------------------------------------------------------------------- COND=HighDose Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 288 -0.0493125 1.4673853 -4.5450000 4.1820000 WTKGCH3 281 -0.0093060 1.5957434 -5.4550000 5.6820000 WTKGCH4 259 0.0777645 1.7810518 -5.0000000 6.1360000 WTKGCH5 242 0.4613967 2.2975734 -6.7270000 10.9090000 ------------------------------------------------------------------------------- COND=Placebo The MEANS Procedure Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 240 0.7376917 1.3046634 -3.8640000 4.5450000 WTKGCH3 240 1.2200208 1.4966333 -2.5000000 8.6360000 WTKGCH4 240 1.6433917 1.6732090 -3.8640000 7.7270000 WTKGCH5 240 1.9458625 2.0153420 -6.0910000 8.8640000 ------------------------------------------------------------------------------- COND=LowDose Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 236 0.1003729 1.2689539 -4.0910000 5.9090000 WTKGCH3 236 0.1401610 1.4197710 -4.3180000 6.5000000 WTKGCH4 236 0.3152288 1.6352823 -7.0450000 5.2730000 WTKGCH5 236 0.8414449 1.9216063 -7.0450000 5.0000000 ------------------------------------------------------------------------------- COND=HighDose Variable N Mean Std Dev Minimum Maximum ------------------------------------------------------------------------------- WTKGCH2 239 -0.0215774 1.4953592 -4.5450000 4.0910000 WTKGCH3 239 -0.0079707 1.6128224 -5.4550000 5.6820000 WTKGCH4 239 0.0916987 1.7929479 -5.0000000 6.1360000 WTKGCH5 239 0.4637657 2.3071275 -6.7270000 10.9090000 ------------------------------------------------------------------------------- Plot of MEANWTCH*TIME. Symbol is value of COND. MEANWTCH | | 2.0 + | P | | 1.8 + | | | P 1.6 + | | | 1.4 + | | | 1.2 + P | | | 1.0 + | | | 0.8 + L | P | | 0.6 + | | | H 0.4 + | L | | 0.2 + L | L | H | 0.0 + H | H | | -0.2 + | ---+-----------------+-----------------+-----------------+-- 1 2 3 4 TIME Plot of MEANWTCH*TIME. Symbol is value of COND. 2.00 + | P | | | | 1.75 + | | | P | | 1.50 + | | | | | 1.25 + | P | MEANWTCH | | | 1.00 + | | | | L | 0.75 + P | | | | | 0.50 + | H | | | L | 0.25 + | | | L | L H | 0.00 + H | H | | | | -0.25 + ---+-----------------+-----------------+-----------------+-- 1 2 3 4 TIME Prozac Weight Data - data from visits 2 through 5 multivariate repeated measures analyses The GLM Procedure Class Level Information Class Levels Values COND 3 1 2 3 Number of observations 989 NOTE: Observations with missing values will not be included in this analysis. Thus, only 715 observations can be used in this analysis. Repeated Measures Level Information Dependent Variable WTKGCH2 WTKGCH3 WTKGCH4 WTKGCH5 Level of TIME 1 2 3 4 Partial Correlation Coefficients from the Error SSCP Matrix / Prob > |r| DF = 712 WTKGCH2 WTKGCH3 WTKGCH4 WTKGCH5 WTKGCH2 1.000000 0.723842 0.701899 0.617306 <.0001 <.0001 <.0001 WTKGCH3 0.723842 1.000000 0.770861 0.689480 <.0001 <.0001 <.0001 WTKGCH4 0.701899 0.770861 1.000000 0.802360 <.0001 <.0001 <.0001 WTKGCH5 0.617306 0.689480 0.802360 1.000000 <.0001 <.0001 <.0001 TIME_N represents the nth degree polynomial contrast for TIME M Matrix Describing Transformed Variables WTKGCH2 WTKGCH3 WTKGCH4 WTKGCH5 TIME_1 -.6708203932 -.2236067977 0.2236067977 0.6708203932 TIME_2 0.5000000000 -.5000000000 -.5000000000 0.5000000000 TIME_3 -.2236067977 0.6708203932 -.6708203932 0.2236067977 E = Error SSCP Matrix TIME_N represents the nth degree polynomial contrast for TIME TIME_1 TIME_2 TIME_3 TIME_1 1024.80338 158.20555 5.65509 TIME_2 158.20555 459.93302 -0.46597 TIME_3 5.65509 -0.46597 370.72444 Partial Correlation Coefficients from the Error SSCP Matrix of the Variables Defined by the Specified Transformation / Prob > |r| DF = 712 TIME_1 TIME_2 TIME_3 TIME_1 1.000000 0.230438 0.009175 <.0001 0.8068 TIME_2 0.230438 1.000000 -0.001128 <.0001 0.9760 TIME_3 0.009175 -0.001128 1.000000 0.8068 0.9760 Sphericity Tests Mauchly's Variables DF Criterion Chi-Square Pr > ChiSq Transformed Variates 5 0.699289 254.21905 <.0001 Orthogonal Components 5 0.699289 254.21905 <.0001 Manova Test Criteria and Exact F Statistics for the Hypothesis of no TIME Effect H = Type III SSCP Matrix for TIME E = Error SSCP Matrix S=1 M=0.5 N=354 Statistic Value F Value Num DF Den DF Pr > F Wilks' Lambda 0.80019538 59.09 3 710 <.0001 Pillai's Trace 0.19980462 59.09 3 710 <.0001 Hotelling-Lawley Trace 0.24969479 59.09 3 710 <.0001 Roy's Greatest Root 0.24969479 59.09 3 710 <.0001 Manova Test Criteria and F Approximations for the Hypothesis of no TIME*COND Effect H = Type III SSCP Matrix for TIME*COND E = Error SSCP Matrix S=2 M=0 N=354 Statistic Value F Value Num DF Den DF Pr > F Wilks' Lambda 0.91682904 10.50 6 1420 <.0001 Pillai's Trace 0.08347056 10.32 6 1422 <.0001 Hotelling-Lawley Trace 0.09038912 10.69 6 944.89 <.0001 Roy's Greatest Root 0.08661644 20.53 3 711 <.0001 NOTE: F Statistic for Roy's Greatest Root is an upper bound. NOTE: F Statistic for Wilks' Lambda is exact. Manova Test Criteria and Exact F Statistics for the Hypothesis of no TIME*low vs placebo Effect H = Contrast SSCP Matrix for TIME*low vs placebo E = Error SSCP Matrix S=1 M=0.5 N=354 Statistic Value F Value Num DF Den DF Pr > F Wilks' Lambda 0.94451457 13.90 3 710 <.0001 Pillai's Trace 0.05548543 13.90 3 710 <.0001 Hotelling-Lawley Trace 0.05874492 13.90 3 710 <.0001 Roy's Greatest Root 0.05874492 13.90 3 710 <.0001 Manova Test Criteria and Exact F Statistics for the Hypothesis of no TIME*high vs placebo Effect H = Contrast SSCP Matrix for TIME*high vs placebo E = Error SSCP Matrix S=1 M=0.5 N=354 Statistic Value F Value Num DF Den DF Pr > F Wilks' Lambda 0.93272452 17.07 3 710 <.0001 Pillai's Trace 0.06727548 17.07 3 710 <.0001 Hotelling-Lawley Trace 0.07212792 17.07 3 710 <.0001 Roy's Greatest Root 0.07212792 17.07 3 710 <.0001 Tests of Hypotheses for Between Subjects Effects Source DF Type III SS Mean Square F Value Pr > F COND 2 861.639174 430.819587 49.02 <.0001 Error 712 6258.014062 8.789346 Contrast DF Contrast SS Mean Square F Value Pr > F low vs placebo 1 512.2737000 512.2737000 58.28 <.0001 high vs placebo 1 754.7494244 754.7494244 85.87 <.0001 Univariate Tests of Hypotheses for Within Subject Effects Source DF Type III SS Mean Square F Value Pr > F TIME 3 263.571138 87.857046 101.14 <.0001 TIME*COND 6 53.974977 8.995829 10.36 <.0001 Error(TIME) 2136 1855.460841 0.868661 Adj Pr > F Source G - G H - F TIME <.0001 <.0001 TIME*COND <.0001 <.0001 Error(TIME) Greenhouse-Geisser Epsilon 0.7918 Huynh-Feldt Epsilon 0.7969 Contrast DF Contrast SS Mean Square F Value Pr > F TIME*low vs placebo 3 29.85525355 9.95175118 11.46 <.0001 TIME*high vs placebo 3 46.23263710 15.41087903 17.74 <.0001 Adj Pr > F Contrast G - G H - F TIME*low vs placebo <.0001 <.0001 TIME*high vs placebo <.0001 <.0001 Analysis of Variance of Contrast Variables TIME_N represents the nth degree polynomial contrast for TIME Contrast Variable: TIME_1 Source DF Type III SS Mean Square F Value Pr > F Mean 1 254.327845 254.327845 176.70 <.0001 COND 2 38.491859 19.245930 13.37 <.0001 Error 712 1024.803381 1.439331 Contrast DF Contrast SS Mean Square F Value Pr > F low vs placebo 1 16.18986751 16.18986751 11.25 0.0008 high vs placebo 1 37.18810476 37.18810476 25.84 <.0001 Contrast Variable: TIME_2 Source DF Type III SS Mean Square F Value Pr > F Mean 1 8.7834208 8.7834208 13.60 0.0002 COND 2 14.9267307 7.4633653 11.55 <.0001 Error 712 459.9330191 0.6459733 Contrast DF Contrast SS Mean Square F Value Pr > F low vs placebo 1 13.20620598 13.20620598 20.44 <.0001 high vs placebo 1 8.67545746 8.67545746 13.43 0.0003 Contrast Variable: TIME_3 Source DF Type III SS Mean Square F Value Pr > F Mean 1 0.4598726 0.4598726 0.88 0.3476 COND 2 0.5563867 0.2781934 0.53 0.5863 Error 712 370.7244404 0.5206804 Contrast DF Contrast SS Mean Square F Value Pr > F low vs placebo 1 0.45918006 0.45918006 0.88 0.3480 high vs placebo 1 0.36907489 0.36907489 0.71 0.4001 Least Squares Means WTKGCH2 WTKGCH3 WTKGCH4 WTKGCH5 COND LSMEAN LSMEAN LSMEAN LSMEAN 1 0.73769167 1.22002083 1.64339167 1.94586250 2 0.10037288 0.14016102 0.31522881 0.84144492 3 -0.02157741 -0.00797071 0.09169874 0.46376569