alpha=1 − (ci/100);
sediff=sqrt((p1(1−p1))/n1 +
((p2(1−p2))/n2));
y=probit(alpha/2);
y=abs(y);
LOWERCI = ROUND (((p1 − p2) −
(ysediff)),.01);
UPPERCI = ROUND (((p1 − p2) +
(ysediff)),.01);
proc print;
run;
Figure 7: Program for
proportions test (with
confidence interval for
difference)
*****;
** Filename: COCHRAN.JOB
;
Purpose: Cochran’s Q
test
;
Tests whether
there is a difference between
two or ;
more treatments
in a repeated measures design
when ;
the response
variable is binary, 1 (success)
0 ;
(fail).
;
;
Created: 16 April
1995
;
**;
Appendix 365