Accurate data on electric-dipole transition matrix elements (EDTMs) for bound-bound Rydbergatom transitions become increasingly important in science and technology. Here we compute radial EDTMs of rubidium using the length, velocity and acceleration gauges for electric-dipole-allowed transitions between states with principal and angular-momentum quantum numbers ranging from 15 to 100. Wave-functions are computed based upon model potentials from Marinescu et al., Phys. Rev. A 49, 982 (1994). Length-gauge EDTMs, often used for low-ℓ transitions, are found to deviate from the fundamentally more accurate velocity-gauge EDTMs by relative amounts of up to ∼ 10 −3 . We discuss the physical reasons for the observed gauge differences, explain the conditions for applicability of the velocity and length gauges for different transition series, and present a decision tree of how to choose EDTMs. Implications for contemporary Rydberg-atom applications are discussed.