Displaying 1 result from an estimated 1 matches for "ncparm".
Did you mean:
charm
2007 Oct 31
0
set initial parameter values for GLMM estimation
...to
predicted scores (no error)
proc mixed data = dataSet;
class treat group;
model y = treat*time ;
random intercept / subject=group(treat);
PARMS (2.1) (1.2) / NOITER;
contrast 'slopes' treat*time 1 -1 0,treat*time 1 0 -1;
ods output contrasts=c;
run;
data dataSet;
set c;
alpha=0.05;
ncparm=numdf*fvalue;
fc=finv(1-alpha,numdf,dendf,0);
power=1-probf(fc,numdf,dendf,ncparm);
run;
proc print;run;
cheers,
Wilfried
Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm