Modeling, optimization, process monitoring, and product development in a toner process using multiway principal component analysis and multiway partial least square method is described. Process measurements and product quality values of past successful batches were collected in a data matrix and preprocessed through time alignment, centering, and scaling. Following the identification of latent variables, an empirical model was built through a fourfold cross validation that can represent the operation of a successful batch. The prepared model provided a realistic prediction of process behavior, realistically represented the operation of the industrial unit, and is mathematically simple enough to be used in online optimization and for automatic control strategies of selected abnormal batches.