The partial slip effects for radiative convective nanofluid flow over a stretching sheet in porous medium are analytically explored in this work. The Navier–Stokes equations, the momentum and the energy equations are converted into a set of non-linear ODEs by the similarity transformation. Using the modified optimal homotopy asymptotic method (OHAM), the resulting non-linear ODEs are analytically approximately solved. The impact of various parameters, such as: the velocity exponential factor n, the wall thickness parameter γ, the dimensionless velocity slip parameter δ1, the Prandtl number Pr, the radiation parameter R, and the dimensionless temperature jump parameter δ2, on the behaviour of the mass and heat transfer is presented. The influence of these parameters is tabular and graphically presented. An excellent agreement between the approximate analytical solution and the corresponding numerical solution is highlighted. The results obtained confirm that modified OHAM is a useful and competitive mathematical tool to explore a large class of non-linear problems with applications in various fields of science and engineering.