exNicolaBullockRuntimeProgressionBootstrapCorrlDiff.docVIP

*ex Nicola Bullock Runtime Progression Bootstrap Corrl Diff SAS.doc; *and ex Nick West probiotics; *ex GIH; *ex Andrea OPTIONS FORMCHAR=||+|+=|-/\*; options ls=90 ps=52 pageno=1; options nonotes nodate nonumber nostimer; *options notes date number stimer; options number; title; %macro assign; *generate data; data dat1; do sim=1 to NoSims; do ID=1 to ss; if IDss/2+1 then do; RandomizedGroup=control; PairRand=ID; end; else do; RandomizedGroup=treat; PairRand=ID-ss/2; end; SubjTrueValue=TrueMean+rannor(0)*TrueSD; Pre=SubjTrueValue+rannor(0); Post=SubjTrueValue+rannor(0); Change=Post-Pre; output; end; end; proc sort; by sim Pre; data dat2; set dat1; select (mod(_n_,4)); when (0,1) RankedGroup=control; when (2,3) RankedGroup=treat; otherwise; end; PairRank=round(mod(_n_-1,ss)+1,2)/2; proc print data=dat2; where sim=1 and ss=20; run; *random assign, no adjustment for baseline; ods listing close; proc mixed data=dat2 covtest cl alpha=0.1; class RandomizedGroup; model Change=RandomizedGroup/ddfm=sat outp=pred alpha=0.1 alphap=0.1 s residual; repeated/group=RandomizedGroup; estimate treat-control RandomizedGroup -1 1/cl alpha=0.1; lsmeans RandomizedGroup/diff cl alpha=0.1; by sim; ods output estimates=est; ods output classlevels=clev; ods output solutionr=solr; ods output solutionf=solf; ods output covparms=cov; ods output lsmeans=lsm; ods output diffs=lsmdiffrandnoadj; run; ods listing; *random assign, adjustment for baseline; ods listing close; proc mixed data=dat2 covtest cl alpha=0.1; class RandomizedGroup; model Change=RandomizedGroup RandomizedGroup*Pre/ddfm=sat outp=pred alpha=0.1 alphap=0.1 s residual; repeated/group=RandomizedGroup; estimate treat-control RandomizedGroup -1 1/cl alpha=0.1; lsmeans RandomizedGroup/diff cl alpha=0.1; by sim; ods output estimates=est; ods output classlevels=clev; ods output solutionr=solr; ods output solutionf=solf; ods output covparms=cov; ods output l