- double q1 = rs[i].auc / ( 2 - rs[i].auc);
- double q2 = 2 * pow2 (rs[i].auc) / ( 1 + rs[i].auc);
-
- se = rs[i].auc * (1 - rs[i].auc) + (rs[i].n1 - 1) * (q1 - pow2 (rs[i].auc)) +
- (rs[i].n2 - 1) * (q2 - pow2 (rs[i].auc));
+ se = rs[i].auc * (1 - rs[i].auc) + (rs[i].n1 - 1) * (rs[i].q1hat - pow2 (rs[i].auc)) +
+ (rs[i].n2 - 1) * (rs[i].q2hat - pow2 (rs[i].auc));