130 likes | 280 Views
Comparison of Repeated Measures and Covariance Analysis for Pretest-Posttest Data. -By Chunmei Zhou. Introduction. We develop a comparison of repeated measures and covariance analysis for pretest-posttest data.
E N D
Comparison of Repeated Measures and Covariance Analysis for Pretest-Posttest Data -By Chunmei Zhou
Introduction • We develop a comparison of repeated measures and covariance analysis for pretest-posttest data. • We consider a study design for which subjects are randomized to a drug or placebo group and measured for systolic blood pressure before and after receiving the treatment. • To develop models to assess the effect of treatment on SBP, we first consider baseline blood pressure as a covariate. As an alternative, the baseline and final blood pressure can be considered to be repeated measures.
Data Description • We simulate data for a pretest-posttest study for comparison of covariance analysis and repeated measures. • We assume the response variable to be normally distributed, the sample size is 30, the average response for baseline is 175 and the treatment effect is –5.
data sbp; %let basemean=175; *overall average response for baseline; %let nsub=30; *Number of subjects; %let err_v=100; *Residual error variance at a time point; %let tsubv=100; *Treatment by subject variance; %let subv=200; *Variance of the subject effects; %let nrep=1000; *Number of replications of the simulation; %let effp1=0; *Placebo effect; %let effp2=-5; *Treatment effect; do trial=1 to &nrep; do subj=1 to ⊄ m=&basemean; sub=sqrt(&subv)*rannor(3345); do time=1 to 2; v=sqrt(&err_v)*rannor(3345); if subj<=&nsub/2 then do; trt=1; if time=1 then y=m+sub+v; if time=2 then y=m+sub+&effp1+v; end; if subj>&nsub/2 then do; trt=2; if time=1 then y=m+sub+v; if time=2 then y=m+sub+&effp2+v+sqrt(&tsubv)*rannor(3345); end; output; end; end; end;
Covariance Model Fitting • We use placebo as reference group to fit the covariance model using the following statement: procmixed data=a; by trial; class subj tref; model sbp2=sbp1 tref/s; make 'solutionf' out=est1; run;
Covariance Model Fitting Covariance Parameter Estimates Cov Parm Estimate Residual 207.11 Solution for Fixed Effects Standard Effect tref Estimate Error DF t Value Pr > |t| Intercept 66.6827 31.0513 27 2.15 0.0409 sbp1 0.6280 0.1834 27 3.42 0.0020 tref 0 -5.9770 5.3231 27 -1.12 0.2714 tref 1 0 . . . .
Repeated-Measures Model Fitting • we fit repeated measures model using placebo and pretest as reference groups. We also create a group variable to have value of 1 for placebo group at post test, 2 for treatment group at post test, and 3 for pretest. procmixed data=b; by trial; class subj tref timref grp; model y=tref timref tref*timref/s; repeated /group=grp; random subj; make 'solutionf' out=est2; run;
Repeated-Measures Model Fitting Covariance Parameter Estimates Cov Parm Group Estimate subj 170.26 Residual grp 1 72.3522 Residual grp 2 259.29 Residual grp 3 49.7515
Repeated-Measures Model Fitting Solution for Fixed Effects Standard Effect tref timref Estimate Error DF t Value Pr >|t| Intercept 168.13 3.8298 28 43.90 <.0001 tref 0 4.6322 5.4162 28 0.86 0.3997 tref 1 0 . . . . timref 0 4.1336 2.8531 28 1.45 0.1585 timref 1 0 . . . . tref*timref 0 0 -7.7003 5.3612 28 -1.44 0.1620 tref*timref 0 1 0 . . . . tref*timref 1 0 0 . . . . tref*timref 1 1 0 . . . .
Comparison of Two Approaches • To compare the two approaches, we set different parameters for three variances (subject variance, residual error variance and treatment by subject variance) to generate data, and then run the simulation 1000 times and get estimates of treatment effect each time from each methods. • However, we could not get model fitting results for the total 1000 simulations from repeated-measures model because of infinite likelihood or convergence failure for some trials.
Method Sub. var. Res. Var. Trt.*Sub. Var. No. of Trials Mean of Est. of trt effect Std. Dev. of Estimate Covariance 200 100 1000 1000 -4.66 9.20 Repeat measures 877 -4.66 9.27 Covariance 200 100 100 1000 -4.97 5.30 Repeat measures 911 -4.99 5.70 Covariance 100 50 100 1000 -4.94 4.14 Repeat measures 962 -4.92 4.39 Covariance 100 100 100 1000 -4.95 5.07 Repeat measures 962 -4.98 5.69 Covariance 200 200 100 1000 -5.00 6.74 Repeat measures 965 -5.05 7.65 Covariance 400 400 100 1000 -5.06 9.21 Repeat measures 970 -5.13 10.51 Table1. Comparison of Estimates for treatment effect from two models with 1000 trials with sample size 30, treatment effect –5 and different variance compontents
Figure1. Box Plots of Estimates for treatment effect Figure2. Box Plots of Estimates for treatment effect from two models with 1000 trials with treatment effect from two models with 1000 trials with treatment –5, subject variance 200, residual error variance 200 effect –5, subject variance 400, residual error var and treatment by subject variance 100 400 and treatment by subject variance 100
Discussion • The repeated-measures analysis of variance and analysis of covariance are two common approaches for analysis of pretest-posttest data. They are closely related but assumptions for the analysis and variance estimates for the parameters are differ. • Our analysis results suggest the covariance analysis method may improve the precision of the estimates of the treatment effects compared to repeated measures analysis. • But we should know the assumption that there is a linear relationship between pre- and post-treatment values may not be true. If this were the case, fitting a baseline covariate could lead to less precise results.