This paper numerically investigates radiative magnetohydrodynamic mixed convection boundary layer flow of nanofluids over a nonlinear inclined stretching/shrinking sheet in the presence of heat source/sink and viscous dissipation. The Rosseland approximation is adopted for thermal radiation effects and the Maxwell-Garnetts and Brinkman models are used for the effective thermal conductivity and dynamic viscosity of the nanofluids respectively. The governing coupled nonlinear momentum and thermal boundary layer equations are rendered into a system of ordinary differential equations via local similarity transformations with appropriate boundary conditions. The non-dimensional, nonlinear, well-posed boundary value problem is then solved with the Keller box implicit finite difference scheme. The emerging thermo-physical dimensionless parameters governing the flow are the magnetic field parameter, volume fraction parameter, power-law stretching parameter, Richardson number, suction/injection parameter, Eckert number and heat source/sink parameter. A detailed study of the influence of these parameters on velocity and temperature distributions is conducted. Additionally the evolution of skin friction coefficient and Nusselt number values with selected parameters is presented. Verification of numerical solutions is achieved via benchmarking with some limiting cases documented in previously reported results, and generally very good correlation is demonstrated. This investigation is relevant to fabrication of magnetic nanomaterials and high temperature treatment of magnetic nano-polymers.