Monday, December 7, 2009

independent t-test


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

proc ttest data=dental;
class treatment;
var resp;
run;

proc mixed data=dental;
class treatment;
model resp=treatment;
run;

proc mixed data=dental;
class treatment;
model resp=treatment / ddfm=Satterthwaite;
repeated / group=treatment;
run;

/* ML: df is not right */
proc nlmixed data=dental;
parms mu=100, tau=10, s2e=2;
if (treatment=1) then eta=mu;
else eta=mu+tau;
model resp ~ normal(eta, s2e);
run;


No comments: