input str5 trt week perc lo hi end
P 0 0 0 0 .
P 3 2 1.25 2.75 .
P 6 -.5 -1.2 .5 .
P 9 1 -.5 1.75 .
P 12 2.25 .5 3.25 2.25
P 15 3 -.25 4.25 3.85
S 0 0 0 0 .
S 3 8 7.25 8.75 .
S 6 7.5 6.5 8.75 .
S 9 7.5 6.5 8.75 .
S 12 7.75 6.5 8.5 7.75
S 15 7.5 5.0 8.25 3.7
B 0 0 0 0 .
B 3 12 11.25 12.75 .
B 6 13 12.25 13.75 .
B 9 13.25 12.5 14.0 .
B 12 13 12.25 13.75 13
B 15 12.85 12.1 13.6 3.6
end
replace end=perc if end==.
twoway (rcap hi lo week) ///
(scatter perc week if trt=="P", clcolor(blue ) connect(l) clwidth(thick)) ///
(scatter perc week if trt=="S", clcolor(green) connect(l) clwidth(thick)) ///
(scatter perc week if trt=="B", clcolor(red ) connect(l) clwidth(thick)) ///
(dot perc week if trt=="P", mcolor(blue ) msize(large)) ///
(dot perc week if trt=="S", mcolor(green) msize(large)) ///
(dot perc week if trt=="B", mcolor(red ) msize(large)) ///
(scatter end week if trt=="P" & week>11, clcolor(blue ) connect(l) clwidth(medthick) clpattern(dash) ) ///
(scatter end week if trt=="S" & week>11, clcolor(green) connect(l) clwidth(medthick) clpattern(dash) ) ///
(scatter end week if trt=="B" & week>11, clcolor(red ) connect(l) clwidth(medthick) clpattern(dash) ) ///
, xlabel( 0 (3) 15 ) legend(on) ///
legend(order(5 "placebo" 6 "Treatment 10mg" 7 "Treatment 20mg") col(3))
附件列表