total = d->posterior[i].probofmaxposterior;
DEBUGPRINT1("start total:%f\n",total);
for(j=1;total<=percentage;j++){
- iofmin = d->posterior[i].iofmaxposterior-j;
- iofmax = d->posterior[i].iofmaxposterior+j;
- total += d->posterior[i].prob[iofmax];
- total += d->posterior[i].prob[iofmin];
+ if(d->posterior[i].iofmaxposterior-j >= 0){
+ iofmin = d->posterior[i].iofmaxposterior-j;
+ total += d->posterior[i].prob[iofmin];
+ }
+ if(d->posterior[i].iofmaxposterior+j <= d->posterior[i].n){
+ iofmax = d->posterior[i].iofmaxposterior+j;
+ total += d->posterior[i].prob[iofmax];
+ }
DEBUGPRINT1("end total:%f\n",total);
}