- for( i=0; i<lgth1; i++ ) for( j=0; j<lgth2; j++ )
- {
- tbppmtx[i][j] = 1.0;
- impmtx2[i][j] = 0.0;
- }
- for( i=0; i<lgth1; i++ ) for( pairpt1=pairprob1[i]; pairpt1->bestpos!=-1; pairpt1++ )
- {
- for( j=0; j<lgth2; j++ ) for( pairpt2=pairprob2[j]; pairpt2->bestpos!=-1; pairpt2++ )
- {
- uido = pairpt1->bestpos;
- ujdo = pairpt2->bestpos;
- prob = pairpt1->bestscore * pairpt2->bestscore;
- if( uido > -1 && ujdo > -1 )
- {
- if( uido > i && j > ujdo )
- {
- impmtx2[i][j] += prob * pairedribosumscore53( nseq1, nseq2, oseq1, oseq2, eff1, eff2, i, j, uido, ujdo ) * consweight_multi;
- tbppmtx[i][j] -= prob;
- }
- else if( i < uido && j < ujdo )
- {
- impmtx2[i][j] += prob * pairedribosumscore35( nseq1, nseq2, oseq1, oseq2, eff1, eff2, i, j, uido, ujdo ) * consweight_multi;
- tbppmtx[i][j] -= prob;
- }
- }
- }
- }
-
-
- for( i=0; i<lgth1; i++ )
- {
- for( j=0; j<lgth2; j++ )
- {
- impmtx2[i][j] += tbppmtx[i][j] * singleribosumscore( nseq1, nseq2, oseq1, oseq2, eff1, eff2, i, j ) * consweight_multi;
- }
- }