The authors have quality controlled six global datasets of drifting buoy data, made comparisons of 15-m drogued and undrogued buoy observations, and developed a 2D linear regression model of the difference between drogued and undrogued drifter velocity as a function of wind. The data were acquired from 2334 Surface Velocity Program (SVP) drifters, including 1845 SVP drifters after they lost their drogues; 704 AN/WSQ-6 Navy drifter buoys; and 503 First Global GARP Experiment (FGGE) drifter buoys. Meridional and zonal surface wind velocity components from the global synoptic FNMOC model, the global synoptic ECMWF model, and the global synoptic NCEP model were interpolated to naval AN/WSQ-6, WOCE-TOGA buoy, or FGGE buoy positions and date/times in the datasets. Two-day mean buoy drift velocities and positions were computed: 122 101 SVP drifter mean velocities before they lost their drogues and 58 201 SVP drifter mean velocities after they lost their drogues, 21 799 Navy drifter mean velocities, and 42 338 FGGE drifter mean velocities. A regression analysis was made on selected data in selected 2(G) lat x 8 degrees long bins: U-undrogued = A(undrogued) + B-undrogued W-undrogued, U-drogued = A(drogued) + B-drogued W-drogued, where U-undrogued, U-drogued was ensemble mean buoy velocity and W was wind velocity, and the real and imaginary parts of these quantities were the zonal and meridional components, respectively. The difference in these complex valued regression coefficients, B-difference = B-undrogued = B-drogued measured the linear response to the wind. Navy and SVPL buoy response to the wind was identical and FGGE buoy response was generally the same; the global weighted mean value of |B-difference| was 0.0088 +/- 0.002. The difference in the complex valued y intercept, A(difference) = A(undrogued) - A(drogued), was nearly always zero within error. The buoy response to wind was also estimated by b = ((U) over bar (undrogued) - (U) over bar (drogued))/(W) over bar (undrogued), where the velocity difference (U) over bar (undrogued) - (U) over bar (drogued) was calculated from ensemble mean buoy velocities and (W) over bar was ensemble mean wind velocity. The global-weighted mean value of |b| was 0.0097 +/- 0.005. The analysis also found that the phase angle of either B-difference or b, which measures the velocity difference with respect to the wind, was zero within error and not a function of the surface wind or the Coriolis parameter.