load d[270,11]=data0.prn; output file=1out0 reset; library optmum; #include optmum.ext; declare k[1]; declare m[1]; declare n[1]; declare mistake[270]; declare casdec[270]; declare pridec[270]; declare stupid[270]; declare pra[270]; declare x[270]; declare xmat1[4,3]; declare xmat2[8,4]; declare xmat3[16,5]; declare xmat4[32,6]; declare xmat5[64,7]; sess=d[1:270,1]; pd=d[1:270,2]; rd=d[1:270,3]; sbj=d[1:270,4]; dec=d[1:270,6]; sig=d[1:270,5]; d1=d[1:270,7]; d2=d[1:270,8]; d3=d[1:270,9]; d4=d[1:270,10]; d5=d[1:270,11]; dc1=dec[1:45]; @ This sets up a matrix xmat1 which has a row corresponding to each @ @ possibility for a round two person. The form of the matrix is as follows: @ @ d1 s @ @ 0 0 @ @ 0 1 @ @ 1 0 @ @ 1 1 @ @ Later a third column is added which contains the probability of observing @ @ the particular combination. @ dec1 = 1; do while dec1 <= 2; signal = 1; do while signal <= 2; index=2*(dec1-1) + signal; xmat1[index,1] = dec1 - 1; xmat1[index,2] = signal - 1; signal = signal + 1; endo; dec1 = dec1 + 1; endo; @ The following lines do the same thing for all of the possible @ @ combinations of d1, d2, and s which may be observed by a round three @ @ person. There are eight possible combinations. @ dec1 = 1; do while dec1 <= 2; dec2 = 1; do while dec2 <= 2; signal = 1; do while signal <= 2; index=4*(dec1-1)+2*(dec2-1)+signal; xmat2[index,1] = dec1 - 1; xmat2[index,2] = dec2 - 1; xmat2[index,3] = signal - 1; signal = signal + 1; endo; dec2 = dec2 + 1; endo; dec1 = dec1 + 1; endo; @ The following lines do the same thing for all of the possible @ @ combinations of d1, d2, d3, and s which may be observed by a round three @ @ person. There are sixteen combinations. @ dec1 = 1; do while dec1 <= 2; dec2 = 1; do while dec2 <= 2; dec3 = 1; do while dec3 <= 2; signal = 1; do while signal <= 2; index=8*(dec1-1)+4*(dec2-1)+2*(dec3-1)+signal; xmat3[index,1] = dec1 - 1; xmat3[index,2] = dec2 - 1; xmat3[index,3] = dec3 - 1; xmat3[index,4] = signal - 1; signal = signal + 1; endo; dec3 = dec3 + 1; endo; dec2 = dec2 + 1; endo; dec1 = dec1 + 1; endo; @ This is for round four and there are 32 @ @ combinations. @ dec1 = 1; do while dec1 <= 2; dec2 = 1; do while dec2 <= 2; dec3 = 1; do while dec3 <= 2; dec4 = 1; do while dec4 <= 2; signal = 1; do while signal <= 2; index=16*(dec1-1)+8*(dec2-1)+4*(dec3-1)+2*(dec4-1)+signal; xmat4[index,1] = dec1 - 1; xmat4[index,2] = dec2 - 1; xmat4[index,3] = dec3 - 1; xmat4[index,4] = dec4 - 1; xmat4[index,5] = signal - 1; signal = signal + 1; endo; dec4 = dec4 + 1; endo; dec3 = dec3 + 1; endo; dec2 = dec2 + 1; endo; dec1 = dec1 + 1; endo; @ More stuff for the next round. sixty four choices. @ @ @ dec1 = 1; do while dec1 <= 2; dec2 = 1; do while dec2 <= 2; dec3 = 1; do while dec3 <= 2; dec4 = 1; do while dec4 <= 2; dec5 = 1; do while dec5 <= 2; signal = 1; do while signal <= 2; index=(32*(dec1-1))+(16*(dec2-1))+(8*(dec3-1))+ (4*(dec4-1))+(2*(dec5-1))+(signal); xmat5[index,1] = dec1 - 1; xmat5[index,2] = dec2 - 1; xmat5[index,3] = dec3 - 1; xmat5[index,4] = dec4 - 1; xmat5[index,5] = dec5 - 1; xmat5[index,6] = signal - 1; signal = signal + 1; endo; dec5 = dec5 + 1; endo; dec4 = dec4 + 1; endo; dec3 = dec3 + 1; endo; dec2 = dec2 + 1; endo; dec1 = dec1 + 1; endo; @ The next loop sets pra for the first round people. The value depends on @ @ the individual's signal. @ i=1; do while i<46; if sig[i] == 1; pra[i] = 2/3; endif; if sig[i] == 0; pra[i] = 1/3; endif; x[i] = 2*pra[i] - 1; i=i+1; endo; x = 2*pra-1; @ This section sets x to be equal to 2*pra - 1 and runs OLS for the first @ @ round people to get a starting value for ML. @ let th0[1,1] = 10; proc like(th); local logl, beta, j, i; beta=th; logl=0; @ This fills in the values of the third column of xmat1. @ xmat1[1,3]=exp((2*2/3-1)*beta)/(1+exp((2*2/3-1)*beta)); xmat1[4,3]=exp((2*2/3-1)*beta)/(1+exp((2*2/3-1)*beta)); xmat1[2,3]= 1 - xmat1[1,3]; xmat1[3,3]= 1 -xmat1[4,3]; @ This section computes pra for round 2 people and computes the fourth column @ @ for xmat2. @ i=46; do while i < 91; if d1[i] == 1; if sig[i] == 1; pra[i] = ((xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3)); x[i] = 2*pra[i] - 1; xmat2[8,4]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat2[6,4] = 1 - xmat2[8,4]; endif; if sig[i] == 0; pra[i] = (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat2[7,4]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat2[5,4] = 1 - xmat2[7,4]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat2[4,4]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat2[2,4] = 1 - xmat2[4,4]; endif; if sig[i] == 0; pra[i] = (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat2[3,4]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat2[1,4] = 1 - xmat2[3,4]; endif; endif; i=i+1; endo; @ This section computes pra and x for the persons who exist in the round 3 @ i=91; do while i < 136; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat3[16,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[14,5] = 1 - xmat3[16,5]; endif; if sig[i] == 0; pra[i] = (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat3[15,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[13,5] = 1 - xmat3[15,5]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat3[8,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[6,5] = 1 - xmat3[8,5]; endif; if sig[i] == 0; pra[i] = (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat3[7,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[5,5] = 1 - xmat3[7,5]; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat3[12,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[10,5] = 1 - xmat3[12,5]; endif; if sig[i] == 0; pra[i] = (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat3[11,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[9,5] = 1 - xmat3[11,5]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat3[4,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[2,5] = 1 - xmat3[4,5]; endif; if sig[i] == 0; pra[i] = (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat3[3,5]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat3[1,5] = 1 - xmat3[3,5]; endif; endif; endif; i=i+1; endo; @ This section computes pra and x for the persons who exist in the round 4 @ i=136; do while i < 181; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[32,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[30,6] = 1 - xmat4[32,6]; endif; if sig[i] == 0; pra[i] = (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[31,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[29,6] = 1 - xmat4[31,6]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[16,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[14,6] = 1 - xmat4[16,6]; endif; if sig[i] == 0; pra[i] = (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[15,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[13,6] = 1 - xmat4[15,6]; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[24,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[22,6] = 1 - xmat4[24,6]; endif; if sig[i] == 0; pra[i] = (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[23,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[21,6] = 1 - xmat4[23,6]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[8,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[6,6] = 1 -xmat4[8,6]; endif; if sig[i] == 0; pra[i] = (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[7,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[5,6] = 1 -xmat4[7,6]; endif; endif; endif; endif; if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[28,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[26,6] = 1 - xmat4[28,6]; endif; if sig[i] == 0; pra[i] = (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[27,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[25,6] = 1 - xmat4[27,6]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[12,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[10,6] = 1 - xmat4[12,6]; endif; if sig[i] == 0; pra[i] = (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[11,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[9,6] = 1 - xmat4[11,6]; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[20,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[18,6] = 1 - xmat4[20,6]; endif; if sig[i] == 0; pra[i] = (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[19,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[17,6] = 1 - xmat4[19,6]; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat4[4,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[2,6] = 1 -xmat4[4,6]; endif; if sig[i] == 0; pra[i] = (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat4[3,6]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat4[1,6] = 1 -xmat4[3,6]; endif; endif; endif; endif; i=i+1; endo; @ This is the section for round 5 @ i=181; do while i <= 225; @if sess[i] < 3;@ if d4[i] == 1; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i] = (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[64,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[62,7] = 1 - xmat5[64,7]; endif; if sig[i] == 0; pra[i] = (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[63,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[61,7] = 1 - xmat5[63,7]; endif; endif; @ d1 =1 @ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[32,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[30,7] = 1 - xmat5[32,7]; endif; if sig[i] == 0; pra[i] = (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[31,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[29,7] = 1 - xmat5[31,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 1 @ if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i] = (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[48,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[46,7] = 1 - xmat5[48,7]; endif; if sig[i] == 0; pra[i] = (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[47,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[45,7] = 1 - xmat5[47,7]; endif; endif; @ d1 = 1 @ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[16,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[14,7] = 1 - xmat5[16,7]; endif; if sig[i] == 0; pra[i] = (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[15,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[13,7] = 1 - xmat5[15,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 0 @ endif; @ d3 = 1 @ if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[56,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[54,7] = 1 - xmat5[56,7]; endif; if sig[i] == 0; pra[i] = (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[55,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[53,7] = 1 - xmat5[55,7]; endif; endif; @ d1 = 1@ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[24,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[22,7] = 1 - xmat5[24,7]; endif; if sig[i] == 0; pra[i] = (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[23,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[21,7] = 1 - xmat5[23,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 1 @ if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[40,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[38,7] = 1 - xmat5[40,7]; endif; if sig[i] == 0; pra[i] = (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[39,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[37,7] = 1 - xmat5[39,7]; endif; endif; @ d1 =1 @ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[8,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[6,7] = 1 - xmat5[6,7]; endif; if sig[i] == 0; pra[i] = (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[7,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[5,7] = 1 - xmat5[7,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 0 @ endif; @ d3 = 0 @ endif; @ d4 = 1 @ if d4[i] == 0; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[60,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[58,7] = 1 - xmat5[60,7]; endif; if sig[i] == 0; pra[i] = (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[59,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[57,7] = 1 - xmat5[59,7]; endif; endif; @ d1 = 1 @ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[28,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[26,7] = 1 - xmat5[28,7]; endif; if sig[i] == 0; pra[i] = (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[27,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[25,7] = 1 - xmat5[27,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 1 @ if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[44,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[42,7] = 1 - xmat5[44,7]; endif; if sig[i] == 0; pra[i] = (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[43,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[41,7] = 1 - xmat5[43,7]; endif; endif; @ d1 = 1 @ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[12,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[10,7] = 1 - xmat5[12,7]; endif; if sig[i] == 0; pra[i] = (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[11,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[9,7] = 1 - xmat5[11,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 0 @ endif; @ d3 = 1 @ if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[52,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[50,7] = 1 - xmat5[52,7]; endif; if sig[i] == 0; pra[i] = (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[51,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[49,7] = 1 - xmat5[51,7]; endif; endif; @ d1 = 1@ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[20,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[18,7] = 1 - xmat5[20,7]; endif; if sig[i] == 0; pra[i] = (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[19,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[17,7] = 1 - xmat5[19,7]; endif; endif; @ d1 = 0@ endif; @ d2 = 1 @ if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[36,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[34,7] = 1 - xmat5[36,7]; endif; if sig[i] == 0; pra[i] = (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[35,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[33,7] = 1 - xmat5[35,7]; endif; endif; @ d1 = 1@ if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; xmat5[4,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[2,7] = 1 - xmat5[4,7]; endif; if sig[i] == 0; pra[i] = (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; xmat5[3,7]=exp(x[i]*beta)/(1+exp(x[i]*beta)); xmat5[1,7] = 1 - xmat5[3,7]; endif; endif; @ d1 = 0 @ endif; @ d2 = 0 @ endif; @ d3 = 0 @ endif; @ d4 = 0 @ @ endif; session < 3 @ i=i+1; endo; @ This is the section for round 6 @ i=226; do while i <= 270; @if sess[i] < 3;@ if d5[i] == 1; if d4[i] == 1; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[64,7]*2/3+xmat5[63,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[64,7]*2/3+xmat5[63,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[64,7]*1/3+xmat5[63,7]*2/3)* (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[64,7]*2/3+xmat5[63,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[64,7]*2/3+xmat5[63,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[64,7]*1/3+xmat5[63,7]*2/3)* (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[32,7]*2/3+xmat5[31,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[32,7]*2/3+xmat5[31,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[32,7]*1/3+xmat5[31,7]*2/3)* (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[32,7]*2/3+xmat5[31,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[32,7]*2/3+xmat5[31,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[32,7]*1/3+xmat5[31,7]*2/3)* (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[48,7]*2/3+xmat5[47,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[48,7]*2/3+xmat5[47,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[48,7]*1/3+xmat5[47,7]*2/3)* (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[48,7]*2/3+xmat5[47,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[48,7]*2/3+xmat5[47,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[48,7]*1/3+xmat5[47,7]*2/3)* (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[16,7]*2/3+xmat5[15,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[16,7]*2/3+xmat5[15,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[16,7]*1/3+xmat5[15,7]*2/3)* (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[16,7]*2/3+xmat5[15,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[16,7]*2/3+xmat5[15,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[16,7]*1/3+xmat5[15,7]*2/3)* (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[56,7]*2/3+xmat5[55,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[56,7]*2/3+xmat5[55,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[56,7]*1/3+xmat5[55,7]*2/3)* (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[56,7]*2/3+xmat5[55,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[56,7]*2/3+xmat5[55,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[56,7]*1/3+xmat5[55,7]*2/3)* (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[24,7]*2/3+xmat5[23,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[24,7]*2/3+xmat5[23,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[24,7]*1/3+xmat5[23,7]*2/3)* (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[24,7]*2/3+xmat5[23,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[24,7]*2/3+xmat5[23,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[24,7]*1/3+xmat5[23,7]*2/3)* (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[40,7]*2/3+xmat5[39,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[40,7]*2/3+xmat5[39,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[40,7]*1/3+xmat5[39,7]*2/3)* (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[40,7]*2/3+xmat5[39,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[40,7]*2/3+xmat5[39,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[40,7]*1/3+xmat5[39,7]*2/3)* (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[8,7]*2/3+xmat5[7,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[8,7]*2/3+xmat5[7,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[8,7]*1/3+xmat5[7,7]*2/3)* (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[8,7]*2/3+xmat5[7,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[8,7]*2/3+xmat5[7,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[8,7]*1/3+xmat5[7,7]*2/3)* (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; endif; @ d4 = 1 @ if d4[i] == 0; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[60,7]*2/3+xmat5[59,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[60,7]*2/3+xmat5[59,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[60,7]*1/3+xmat5[59,7]*2/3)* (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[60,7]*2/3+xmat5[59,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[60,7]*2/3+xmat5[59,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[60,7]*1/3+xmat5[59,7]*2/3)* (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[28,7]*2/3+xmat5[27,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[28,7]*2/3+xmat5[27,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[28,7]*1/3+xmat5[27,7]*2/3)* (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[28,7]*2/3+xmat5[27,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[28,7]*2/3+xmat5[27,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[28,7]*1/3+xmat5[27,7]*2/3)* (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[44,7]*2/3+xmat5[43,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[44,7]*2/3+xmat5[43,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[44,7]*1/3+xmat5[43,7]*2/3)* (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[44,7]*2/3+xmat5[43,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[44,7]*2/3+xmat5[43,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[44,7]*1/3+xmat5[43,7]*2/3)* (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[12,7]*2/3+xmat5[11,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[12,7]*2/3+xmat5[11,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[12,7]*1/3+xmat5[11,7]*2/3)* (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[12,7]*2/3+xmat5[11,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[12,7]*2/3+xmat5[11,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[12,7]*1/3+xmat5[11,7]*2/3)* (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[52,7]*2/3+xmat5[51,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[52,7]*2/3+xmat5[51,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[52,7]*1/3+xmat5[51,7]*2/3)* (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[52,7]*2/3+xmat5[51,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[52,7]*2/3+xmat5[51,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[52,7]*1/3+xmat5[51,7]*2/3)* (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[20,7]*2/3+xmat5[19,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[20,7]*2/3+xmat5[19,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[20,7]*1/3+xmat5[19,7]*2/3)* (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[20,7]*2/3+xmat5[19,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[20,7]*2/3+xmat5[19,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[20,7]*1/3+xmat5[19,7]*2/3)* (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[36,7]*2/3+xmat5[35,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[36,7]*2/3+xmat5[35,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[36,7]*1/3+xmat5[35,7]*2/3)* (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[36,7]*2/3+xmat5[35,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[36,7]*2/3+xmat5[35,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[36,7]*1/3+xmat5[35,7]*2/3)* (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[4,7]*2/3+xmat5[3,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[4,7]*2/3+xmat5[3,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[4,7]*1/3+xmat5[3,7]*2/3)* (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[4,7]*2/3+xmat5[3,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[4,7]*2/3+xmat5[3,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[4,7]*1/3+xmat5[3,7]*2/3)* (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; endif; @ d4 = 0 @ endif; @ d5 = 1 @ if d5[i] == 0; if d4[i] == 1; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[62,7]*2/3+xmat5[61,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[62,7]*2/3+xmat5[61,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[62,7]*1/3+xmat5[61,7]*2/3)* (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[62,7]*2/3+xmat5[61,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[62,7]*2/3+xmat5[61,7]*1/3)* (xmat4[32,6]*2/3+xmat4[31,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[62,7]*1/3+xmat5[61,7]*2/3)* (xmat4[32,6]*1/3+xmat4[31,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[30,7]*2/3+xmat5[29,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[30,7]*2/3+xmat5[29,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[30,7]*1/3+xmat5[29,7]*2/3)* (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[30,7]*2/3+xmat5[29,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[30,7]*2/3+xmat5[29,7]*1/3)* (xmat4[16,6]*2/3+xmat4[15,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[30,7]*1/3+xmat5[29,7]*2/3)* (xmat4[16,6]*1/3+xmat4[15,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[46,7]*2/3+xmat5[45,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[46,7]*2/3+xmat5[45,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[46,7]*1/3+xmat5[45,7]*2/3)* (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[46,7]*2/3+xmat5[45,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[46,7]*2/3+xmat5[45,7]*1/3)* (xmat4[24,6]*2/3+xmat4[23,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[46,7]*1/3+xmat5[45,7]*2/3)* (xmat4[24,6]*1/3+xmat4[23,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[14,7]*2/3+xmat5[13,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[14,7]*2/3+xmat5[13,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[14,7]*1/3+xmat5[13,7]*2/3)* (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[14,7]*2/3+xmat5[13,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[14,7]*2/3+xmat5[13,7]*1/3)* (xmat4[8,6]*2/3+xmat4[7,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[14,7]*1/3+xmat5[13,7]*2/3)* (xmat4[8,6]*1/3+xmat4[7,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[54,7]*2/3+xmat5[53,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[54,7]*2/3+xmat5[53,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[54,7]*1/3+xmat5[53,7]*2/3)* (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[54,7]*2/3+xmat5[53,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[54,7]*2/3+xmat5[53,7]*1/3)* (xmat4[28,6]*2/3+xmat4[27,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[54,7]*1/3+xmat5[53,7]*2/3)* (xmat4[28,6]*1/3+xmat4[27,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[22,7]*2/3+xmat5[21,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[22,7]*2/3+xmat5[21,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[22,7]*1/3+xmat5[21,7]*2/3)* (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[22,7]*2/3+xmat5[21,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[22,7]*2/3+xmat5[21,7]*1/3)* (xmat4[12,6]*2/3+xmat4[11,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[22,7]*1/3+xmat5[21,7]*2/3)* (xmat4[12,6]*1/3+xmat4[11,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[38,7]*2/3+xmat5[37,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[38,7]*2/3+xmat5[37,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[38,7]*1/3+xmat5[37,7]*2/3)* (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[38,7]*2/3+xmat5[37,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[38,7]*2/3+xmat5[37,7]*1/3)* (xmat4[20,6]*2/3+xmat4[19,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[38,7]*1/3+xmat5[37,7]*2/3)* (xmat4[20,6]*1/3+xmat4[19,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[6,7]*2/3+xmat5[5,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[6,7]*2/3+xmat5[5,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[6,7]*1/3+xmat5[5,7]*2/3)* (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[6,7]*2/3+xmat5[5,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[6,7]*2/3+xmat5[5,7]*1/3)* (xmat4[4,6]*2/3+xmat4[3,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[6,7]*1/3+xmat5[5,7]*2/3)* (xmat4[4,6]*1/3+xmat4[3,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; endif; @ d4 = 1 @ if d4[i] == 0; if d3[i] == 1; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[58,7]*2/3+xmat5[57,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[58,7]*2/3+xmat5[57,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[58,7]*1/3+xmat5[57,7]*2/3)* (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[58,7]*2/3+xmat5[57,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[58,7]*2/3+xmat5[57,7]*1/3)* (xmat4[30,6]*2/3+xmat4[29,6]*1/3)* (xmat3[16,5]*2/3+xmat3[15,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[58,7]*1/3+xmat5[57,7]*2/3)* (xmat4[30,6]*1/3+xmat4[29,6]*2/3)* (xmat3[16,5]*1/3+xmat3[15,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[26,7]*2/3+xmat5[25,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[26,7]*2/3+xmat5[25,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[26,7]*1/3+xmat5[25,7]*2/3)* (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[26,7]*2/3+xmat5[25,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[26,7]*2/3+xmat5[25,7]*1/3)* (xmat4[14,6]*2/3+xmat4[13,6]*1/3)* (xmat3[8,5]*2/3+xmat3[7,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[26,7]*1/3+xmat5[25,7]*2/3)* (xmat4[14,6]*1/3+xmat4[13,6]*2/3)* (xmat3[8,5]*1/3+xmat3[7,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[42,7]*2/3+xmat5[41,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[42,7]*2/3+xmat5[41,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[42,7]*1/3+xmat5[41,7]*2/3)* (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[42,7]*2/3+xmat5[41,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[42,7]*2/3+xmat5[41,7]*1/3)* (xmat4[22,6]*2/3+xmat4[21,6]*1/3)* (xmat3[12,5]*2/3+xmat3[11,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[42,7]*1/3+xmat5[41,7]*2/3)* (xmat4[22,6]*1/3+xmat4[21,6]*2/3)* (xmat3[12,5]*1/3+xmat3[11,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[10,7]*2/3+xmat5[9,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[10,7]*2/3+xmat5[9,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[10,7]*1/3+xmat5[9,7]*2/3)* (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[10,7]*2/3+xmat5[9,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[10,7]*2/3+xmat5[9,7]*1/3)* (xmat4[6,6]*2/3+xmat4[5,6]*1/3)* (xmat3[4,5]*2/3+xmat3[3,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[10,7]*1/3+xmat5[9,7]*2/3)* (xmat4[6,6]*1/3+xmat4[5,6]*2/3)* (xmat3[4,5]*1/3+xmat3[3,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; if d3[i] == 0; if d2[i] == 1; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[50,7]*2/3+xmat5[49,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[50,7]*2/3+xmat5[49,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[50,7]*1/3+xmat5[49,7]*2/3)* (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[50,7]*2/3+xmat5[49,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[50,7]*2/3+xmat5[49,7]*1/3)* (xmat4[26,6]*2/3+xmat4[25,6]*1/3)* (xmat3[14,5]*2/3+xmat3[13,5]*1/3)* (xmat2[8,4]*2/3+xmat2[7,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[50,7]*1/3+xmat5[49,7]*2/3)* (xmat4[26,6]*1/3+xmat4[25,6]*2/3)* (xmat3[14,5]*1/3+xmat3[13,5]*2/3)* (xmat2[8,4]*1/3+xmat2[7,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[18,7]*2/3+xmat5[17,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[18,7]*2/3+xmat5[17,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[18,7]*1/3+xmat5[17,7]*2/3)* (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[18,7]*2/3+xmat5[17,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[18,7]*2/3+xmat5[17,7]*1/3)* (xmat4[10,6]*2/3+xmat4[9,6]*1/3)* (xmat3[6,5]*2/3+xmat3[5,5]*1/3)* (xmat2[4,4]*2/3+xmat2[3,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[18,7]*1/3+xmat5[17,7]*2/3)* (xmat4[10,6]*1/3+xmat4[9,6]*2/3)* (xmat3[6,5]*1/3+xmat3[5,5]*2/3)* (xmat2[4,4]*1/3+xmat2[3,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; if d2[i] == 0; if d1[i] == 1; if sig[i] == 1; pra[i]= (xmat5[34,7]*2/3+xmat5[33,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3/ ((xmat5[34,7]*2/3+xmat5[33,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*2/3+ (xmat5[34,7]*1/3+xmat5[33,7]*2/3)* (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[34,7]*2/3+xmat5[33,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3/ ((xmat5[34,7]*2/3+xmat5[33,7]*1/3)* (xmat4[18,6]*2/3+xmat4[17,6]*1/3)* (xmat3[10,5]*2/3+xmat3[9,5]*1/3)* (xmat2[6,4]*2/3+xmat2[5,4]*1/3)* (xmat1[4,3]*2/3+xmat1[3,3]*1/3)*1/3+ (xmat5[34,7]*1/3+xmat5[33,7]*2/3)* (xmat4[18,6]*1/3+xmat4[17,6]*2/3)* (xmat3[10,5]*1/3+xmat3[9,5]*2/3)* (xmat2[6,4]*1/3+xmat2[5,4]*2/3)* (xmat1[4,3]*1/3+xmat1[3,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; if d1[i] == 0; if sig[i] == 1; pra[i] = (xmat5[2,7]*2/3+xmat5[1,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3/ ((xmat5[2,7]*2/3+xmat5[1,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*2/3+ (xmat5[2,7]*1/3+xmat5[1,7]*2/3)* (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*1/3); x[i] = 2*pra[i] - 1; endif; if sig[i] == 0; pra[i] = (xmat5[2,7]*2/3+xmat5[1,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3/ ((xmat5[2,7]*2/3+xmat5[1,7]*1/3)* (xmat4[2,6]*2/3+xmat4[1,6]*1/3)* (xmat3[2,5]*2/3+xmat3[1,5]*1/3)* (xmat2[2,4]*2/3+xmat2[1,4]*1/3)* (xmat1[2,3]*2/3+xmat1[1,3]*1/3)*1/3+ (xmat5[2,7]*1/3+xmat5[1,7]*2/3)* (xmat4[2,6]*1/3+xmat4[1,6]*2/3)* (xmat3[2,5]*1/3+xmat3[1,5]*2/3)* (xmat2[2,4]*1/3+xmat2[1,4]*2/3)* (xmat1[2,3]*1/3+xmat1[1,3]*2/3)*2/3); x[i] = 2*pra[i] - 1; endif; endif; endif; endif; endif; @ d4 = 0 @ endif; @ d5 = 0 @ @endif; session < 3@ i=i+1; endo; j=1; do while j<=270; logl = logl - dec[j]*ln(exp(x[j]*beta)/(1+exp(x[j]*beta))) - (1- dec[j])*ln(1/(1+exp(x[j]*beta))); j=j+1; endo; retp(logl); endp; @ This section sets the options for the optimization subroutine and calls it. @ _opshess=1; _opalgr=5; _opbtol=1.e-6; _opgtol=1.e-6; {th,f,g,retcode}=optmum(&like,th0); @ This sections prints interesting and useful results. @ i=1; do while i<=270; mistake[i] = 0; if pra[i] < .5; if dec[i] == 1; mistake[i] = 1; endif; endif; if pra[i] > .5; if dec[i] == 0; mistake[i] = 1; endif; endif; i=i+1; endo; i=1; do while i<=270; casdec[i] = 0; if sig[i] ne dec[i]; if dec[i] == 1; if pra[i] > .5; casdec[i] = 1; endif; endif; if dec[i] == 0; if pra[i] < .5; casdec[i] = 1; endif; endif; endif; i=i+1; endo; i=1; do while i <= 270; pridec[i] = 0; if dec[i] == 1; if pra[i] < .5; if sig[i] == 1; pridec[i] = 1; endif; endif; endif; if dec[i] == 0; if pra[i] > .5; if sig[i] == 0; pridec[i] = 1; endif; endif; endif; i=i+1; endo; i=1; do while i <=270; stupid[i] = 0; if sig[i] == 0; if pra[i] < .5; if dec[i] == 1; stupid[i] = 1; endif; endif; endif; if sig[i] == 1; if pra[i] > .5; if dec[i] == 0; stupid[i] = 1; endif; endif; endif; i=i+1; endo; c=sess~pd~rd~sbj~d1~d2~d3~d4~d5~sig~dec~pra~mistake~casdec ~pridec~stupid; c; vg=inv(hessp(&like,th)); t=th./sqrt(diag(vg)); m = cdftc(t,i-2); if m < .00001; m = 0; endif; output on; "Output file for logit estimation of information cascade data with one beta."; "";"";" Final Estimates "; "";"";"Sample size: " round(i-1); "";"Starting value: " th0; "";"Value of the likelihood function:" f; "";"Beta: " th; "";"Standard Error: " sqrt(diag(vg)); "";"t-score: " t; "";"Significance level: " m; output off; output on; "";"stupid decisions";""; i=1; do while i < 271; if stupid[i] == 1; c[i,.]; endif; i=i+1; endo; "";"private decisions";""; i=1; do while i < 271; if pridec[i] == 1; c[i,.]; endif; i=i+1; endo; "";"cascade decisions";""; i=1; do while i < 271; if casdec[i] == 1; c[i,.]; endif; i=i+1; endo; output off; output on; "";"session period round subject d1 d2 d3 d4 d5 sig dec pra mis casdec pridec stupid"; ""; c; output off; end; @ output on; "";xmat1; "";xmat2; "";xmat3; "";xmat4; "";xmat5; output off; @