PROC TITLE 'Fidelity and Migration rates of Black Ducks'; PROC MODEL NPAR=20; INLINE DECLARE REAL*8 F1N(2), F2N(2), F1S(2), F2S(2), PI1(2); INLINE DECLARE REAL*8 PI2(2), R1(2), R2(2), B1(2), B2(2); INLINE GRPS=2; INLINE WRITE(*,*); INLINE DO i=0, GRPS-1; INLINE F1N(i+1)=S(i*10+1); INLINE F2N(i+1)=S(i*10+2); INLINE F1S(i+1)=S(i*10+3); INLINE F2S(i+1)=S(i*10+4); INLINE PI1(i+1)=S(i*10+5); INLINE PI2(i+1)=S(i*10+6); INLINE R1(i+1)=S(i*10+7); INLINE R2(i+1)=S(i*10+8); INLINE B1(i+1)=S(i*10+9); INLINE B2(i+1)=S(i*10+10); INLINE END DO; inline write(*,'(10f15.10)')(s(i),i=1,20); INLINE DO i=0, GRPS-1; INLINE j=i+1; INLINE P(1,i*2+1)= (1-B1(j))*F1N(j); INLINE P(2,i*2+1)= B1(j)*F2N(j); INLINE P(3,i*2+1)= (1.-F1N(j))*(1.-PI1(j))*F1S(j); INLINE P(4,i*2+1)= (1.-F1N(j))*PI1(j)*F2S(j); INLINE P(5,i*2+1)= 1-(P(1,i*2+1)+P(2,i*2+1)+P(3,i*2+1)+P(4,i*2+1)); INLINE P(1,i*2+2)= (1-B2(j))*F2N(j); INLINE P(2,i*2+2)= B2(j)*F1N(j); INLINE P(3,i*2+2)= (1.-F2N(j))*PI2(j)*F1S(j); INLINE P(4,i*2+2)= (1.-F2N(j))*(1.-PI2(j))*F2S(j); INLINE P(5,i*2+2)= 1-(P(1,i*2+2)+P(2,i*2+2)+P(3,i*2+2)+P(4,i*2+2)); INLINE WRITE(*,'(10F15.10)')(P(K,i*2+1),K=1,5),(p(k,i*2+2),k=1,5); INLINE END DO; COHORT=306; 8:; 0:; 7:; 0:; COHORT=225; 1:; 15:; 6:; 8:; COHORT=253; 7:; 1:; 8:; 0:; COHORT=342; 0:; 7:; 5:; 12:; LABELS; S(1)=F1N_AF(66); S(2)=F2N_AF(66); S(3)=F1S_AF(66); S(4)=F2S_AF(66); S(5)=PI1_AF(66); S(6)=PI2_AF(66); S(7)=R1_AF(66); S(8)=R2_AF(66); S(9)=B1_AF(66); S(10)=B2_AF(66); S(11)=F1N_AM(66); S(12)=F2N_AM(66); S(13)=F1S_AM(66); S(14)=F2S_AM(66); S(15)=PI1_AM(66); S(16)=PI2_AM(66); S(17)=R1_AM(66); S(18)=R2_AM(66); S(19)=B1_AM(66); S(20)=B2_AM(66); PROC ESTIMATE NOVAR MAXFN=50000 NAME=GENERAL; initial; s(1)=.975; s(2)=.002; s(3)=.913323; s(4)=.037; s(5)=.01; s(6)=.03; s(9)=.97; s(10)=.07; CONSTRAINTS; S(7)= 306; S(8)= 225; S(17)= 253; S(18)= 342; PROC STOP;