I had add it. But it dosnt work. Perhaps I made an other mistake. My hole SAS code:
data a;
INPUT obs i;
CARDS;
1 1
2 0
3 1
4 1
5 0
6 0
7 0
8 1
9 0
10 1
run;
data a;
set a;
i_sum+i;
run;
data a;
set b;
%let a=0.001;
%let b=0.2;
%let p1=0.1;
%let p2=0.2;
%let q=log(&p2*(1-&p1)/&p1*(1-&p2));
%let h1=log(&b/(1-&a))/&q;
%let h2=log((1-&b)/&a)/&q;
q=&q;
h1=&h1;
h2=&h2;
s=log((1-&p1)/(1-&p2))/&q;
d1=&h1+s*obs;
d2=&h2+s*obs;
run;
proc gplot data=a;
title 'test';
symbol1 i=stepjl repeat=1;
symbol2 i=join;
plot i_sum*obs=1 (d1 d2)*obs=2 /overlay;
run;