In this paper, we extend an alternative, phenomenological approach to inflation by means of an equation of state and a sound speed, both of them functions of the number of $e$-folds and three phenomenological parameters. This approach captures a number of possible inflationary models, including those with non-canonical kinetic terms or scale-dependent non-gaussianities. We perform Markov Chain Monte Carlo analyses using the latest cosmological publicly available measurements, which include Cosmic Microwave Background (CMB) data from the Planck satellite. In the context of the phenomenological description studied here, the sound speed of the primordial curvature perturbation is found to be $c_s> 0.2$ at $95\%$~CL for a number of $e$-folds $N=60$. Within this parametrization, we discard scale invariance with a significance of about $10\sigma$, and the running of the spectral index is constrained as $\alpha_s=-0.62\,^{+0.08}_{-0.09} \times 10^{-3}$. The limit on the tensor-to-scalar ratio is $r<0.003$ at $95\%$~CL from CMB data alone. The maximum amplitude of the equilateral non-gaussianity that we obtain, $|f^{\text{equil}}_{\text{NL}}|\sim 2$, is much smaller than the current Planck mission errors. Future high-redshift, all-sky surveys could reach the required accuracy on equilateral non-gaussianities to allow for additional and independent tests of the parametrization explored here.

We revisit the cosmological viability of the Hu $\&$ Sawicki modified gravity scenario. The impact of such a modification on the different cosmological observables, including gravitational waves, is carefully described. The most recent cosmological data, as well as constraints on the relationship between the clustering parameter $\sigma_8$ and the current matter mass-energy density $\Omega_m$ from cluster number counts and weak lensing tomography, are considered in our numerical calculations. The strongest bound we find is $|f_{R0}| < 3.7 \times 10^{-6}$ at $95\%$~CL. Forthcoming cluster surveys covering 10,000 deg$^2$ in the sky, with galaxy surface densities of $\mathcal{O}(10)$~arcmin$^{-2}$ could improve the precision in the $\sigma_8$-$\Omega_m$ relationship, tightening the above constraint.