Monday, December 7, 2009

paired t-test


proc import datafile="F:\publish\data analysis using R\data\CBT.csv" out=cbt dbms=csv replace;
run;

proc sql;
select t(diff) as t, prt(diff) as pvalue
from cbt
;quit;

proc sql;
create table cbt_long as
select monotonic() as id, 1 as period, prewt as weight from cbt union
select monotonic() as id, 2 as period, postwt as weigt from cbt
;quit;

proc mixed data=cbt_long;
class id period;
model weight = period;
random id;
run;

No comments: