Constraining the Evolutionary Stage of Class I Protostars: Multiwavelength Observations and Modeling