The results of extensive computations are presented to accurately characterize transitions to chaos for the Kuramoto-Sivashinsky equation. In particular we follow the oscillatory dynamics in a window that supports a complete sequence of period doubling bifurcations preceding chaos. As many as 13 period doublings are followed and used to compute the Feigenbaum number for the cascade and so enable an accurate numerical evaluation of the theory of universal behavior of nonlinear systems, for an infinite dimensional dynamical system. Furthermore, the dynamics at the threshold of chaos exhibit a self-similar behavior that is demonstrated and used to compute a universal scaling factor, which arises also from the theory of nonlinear maps and can enable continuation of the solution into a chaotic regime. Aperiodic solutions alternate with periodic ones after chaos sets in, and we show the existence of a period six solution separated by chaotic regions.