Dear Xiaochao, I attach data files with the calculations for the kinematic that you asked for. The resonance region calculations are organized as follows: W (GeV) - Q2 (GeV^2) - APV_LT / Q2 (ppm GeV^-2) - APV/Q2 (ppm GeV^-2)- (APV + delta)/Q2 (ppm GeV^-2)- (APV-delta)/Q2 (ppm GeV^-2) the notation above: APV_LT refers to the Cahn-Gilman formula that comes from parton model and is independent of PDF's - this is to have a benchmark; the last three numbers correspond to the model for F_1,2^{\gamma Z} from fit to e-m DIS data plus isospin rotation, the central value - upper bound - lower bound. These latter calculations do NOT contain the contribution of F_3. The uncertainty is roughly 15% and should compensate for this shortcoming. I believe the uncertainty is conservative but is reasonable. As far as the data in the resonance region are concerned, I am not sure whether I understand the kinematics correctly. I calculated APV/Q2 at fixed Q2 and Eb while varying W2 (and Ep). Is that what is measured? In any case, the comparison with the data seems fine except for the Delta region where, quite surprisingly, the data tends to significantly lower values than the model. Since we believe we understand the Delta quite well this is quite puzzling. For the DIS kinematics I give here just the requested points (W2 - Q2 - APV_LT/Q2 - APV/Q2 - (APV+delta)/Q2 - (APV-delta)/Q2) DIS1: 4.297000e+00 1.085000e+00 -8.438344e+01 -8.508855e+01 -9.725292e+01 -7.245306e+01 DIS2L: 5.430000e+00 1.901000e+00 -8.705350e+01 -8.416709e+01 -9.698139e+01 -7.130367e+01 RES3: 1.595000e+00 9.500000e-01 -8.114092e+01 -9.362097e+01 -1.028383e+02 -8.310811e+01 RES4: 2.530000e+00 8.310000e-01 -8.245911e+01 -8.458297e+01 -9.792383e+01 -6.676043e+01 RES5: 3.450000e+00 7.570000e-01 -8.277284e+01 -8.598995e+01 -1.004504e+02 -7.032773e+01 RES7: 3.923000e+00 1.472000e+00 -8.440149e+01 -8.473894e+01 -9.798336e+01 -7.107142e+01 RES7b 4.122000e+00 1.278000e+00 -8.440843e+01 -8.494511e+01 -9.764983e+01 -7.183164e+01 Best, Misha